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RAPIER - A FORTRAN IV PROGRAM FOR MULTIPLE LINEAR REGRESSION 
ANALYSIS PROVIDING INTERNALLY EVALUATED REMODELING 
by Steven M. Sidik and Bert Henry 
Lewis Research Center 

SUMMARY 

RAPIER is a digital computer program which can be used with ease to perform ex- 
tensive regression analyses or a simple least-squares curve fit, and it includes a back- 
ward term rejection option. The program is written in FORTRAN IV, version 13, for 
the IBM 7094/7044 DCS. The major value of the program is its comprehensiveness of 
calculations and options. 

RAPIER computes the variance-covariance matrix of the independent variables, re- 
gression coefficients, t- statistics for individual tests, and analysis of variance tables 
for overall testing of regression. There is a provision for a choice of three strategies 
for the variance estimate to be used in computing t- statistics. 

Also, more than one set of response or dependent variables can by analyzed for the 
same set of independent variables. 

A backward rejection option method based on the first dependent variable may be 
used. In this case, a critical significance level is supplied as input. The least signifi- 
cant independent variable is deleted and the regression recomputed. This process is 
repeated until all remaining variables have significantly nonzero coefficients. 

The algorithm uses the triangular form of symmetric matrices throughout. It also 
allows for the use of weighted regression, computation of predicted values at any com- 
bination of independent variables, a table of residuals, and a chi-square test for the 
normality of the distribution of residuals. 


INTRODUCTION 

RAPIER is an almost entirely new multiple regression computer program. It is the 
result of 5 years of development in meeting the needs of several statistical investigators 
posing a variety of problems. The problems included analysis of nuclear reactor com- 
ponents, determining predictive models from corrosion and fracture data of both metals 
and alloys, investigating the behavior of processing variables in the manufacture of 



solar cells, optimizing fuel-cell experiment procedures, and predicting personnel per- 
formance from academic histories. 

The nucleus of the program is based on a program written by Kunin (ref. 1). How- 
ever, in its present expanded form, it allows the user to choose from a number of sets 
of options which include options of input, of methods for calculation, and of output, 
thereby providing great flexibility. 

With the aid of a few control cards, the program can be used readily for a wide 
range of applications which can vary from a simple least- squares curve-fitting problem 
to a complete regression analysis. It can provide the variance-covariance matrix of 
independent variables, regression coefficients, the variance-covariance matrix of the 
regression coefficients, individual t- statistics with their significance levels, analysis 
of variance tables for significance of regression, special usage of replicated data to 
estimate the error due' to lack of fit, any one of three pooling procedures which may be 
used to estimate the error variance, tests for normality of distribution of the residuals, 
weighted regression, and the use of more than one dependent variable. 

The mathematical analysis of the computations and their reliability is aided further 
by the option of obtaining an eigenvector decomposition of both the variance- covariance 
matrix and the correlation matrix of the independent variables. 

The program also provides an option to perform a backward rejection regression at 
any given level of significance. 

Despite its sophistication, RAPIER is relatively easy to use, but it presupposes 
that the user has at least a basic knowledge and/or experience in the application of sta- 
tistical techniques. 

To provide a framework for the discussion of the calculations and statistical options 
available in RAPIER, a brief description of multiple linear regression is presented, 
with no attempt to make the discussion thorough or rigorous. Notable presentations 
of applied regression analysis are those by Draper and Smith (ref. 2) and Graybill 
(ref. 3). Reference 4 by Kendall and Stuart is a useful guide to both applications and 
theoretical justifications. Rao (ref. 5) presents a more mathematically sophisticated 
treatment of the subject of linear statistical models. 

After discussion of the calculations and options available, the card input necessary 
is described in detail and illustrated by an example which uses almost all of the options. 

* SYMBOLS 


A 

matrix 

A’ 

transpose of A 

A' 1 

inverse of A 
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I 


B 

b 

b i 

b i 

b 0 

b l> 

c 

c ij 

D 



E(x) 


f 

fjfcp • • • > ^j^) 

H 0 

H 1 

J 

K 

k 

LOF 

M 

MS (source) 

m 

N 

N( M ,a 2 ) 

n i 


matrix 

vector (column) 

true regression coefficient 

estimated regression coefficient 

constant term 

unknown parameters 

correlation matrix 

elements of C 

indicator variable, equal to 0 if no bp coefficient is estimated and 
equal to 1 if bp is estimated 

expected value of x (i. e. , average of x over all possible values 
of x) 

vector of observation errors 

statistic distributed as variance ratio with a and b degrees of 
freedom 

expected number of observations in each cell of a partitioned range 
of studentized residuals 

term of regression equation 

statistical hypothesis to be tested 

alternate hypothesis to be accepted if Hq is judged to be false 
number of coefficients estimated, excluding bg 
number of independent variables observed 

number of segments or cells in range of possible studentized resid- 
uals 

lack of fit 

total number of independent and dependent variables 
mean square due to source, where source is REG, RES, etc. 
moment about origin 
number of observations 

o 

normal distribution with mean ^ and variance a 
number of studentized residuals in i n cell 
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Ill 


Pi 

R 

REG 

REP 

RES 

r 

S 


SSQ (source) 
TOT 

*n 

V(x) 

w,x 

X 


probability in i** 1 cell 

number of sets of replicates 

regression 

replication 

residual 

number of replicates in set 
diagonal matrix 

sum of squares correction if D = 1, and 0 if D = 0 
elements of diagonal matrix 

sum of squares due to source, where source is REG, RES, etc. 
total 

statistic distributed as Student’s t with n degrees of freedom 
variance of x, expected value of (x - E(x)) 2 
matrices 
vector (column) 


x(J) 


y 

Z; 


z 


1’ • 


[l x 


M 

; 2 




<v 


J K 


X j 


1 N 

Sl> 

W i=l 


ij 


vector (column) 
studentized residual 
variables 

mean of x defined as E(x) 

estimate of ju based on observation of random sample 

variance of x defined as V(x) 

2 

estimate of a based on observation of random sample 
statistic distributed as chi-square with n degrees of freedom 
is distributed as 


Superscript: 

’ transpose 
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ESTIMATION OF BASIC LINEAR MODEL 


Basic Linear Model 

In multiple linear regression, a dependent or response variable Y (such as tem- 
perature or pressure) measured on an object or experiment is assumed to be correlated 
with a function of one or more other variables (z^, . . z^) measured on the same ob- 

ject or experiment. This function includes a number of unknown parameters 
(bj, . . . , bj) and can be represented as 

y = h(b r . . . , bj, z v . . . , z R ) + e (1) 

The only restriction imposed on this function is that it be linear in the parameters 

J 

y = £ b j f j (z l’ ' • *» Z K } (2) 

j=l 

where fj(Zj, . . z R ) is a TERM of the regression equation. (A TERM is a quantity 

which may be a variable or a function of a variable, e. g. , T is a TERM and Z, after 

it is defined as Z = log T, is also a TERM. ) 

Suppose that there are N observations of the dependent variable. Let the subscript 

i indicate that the values are associated with the i tb observation; in particular, the 

value of the response variable y. would depend on the observed values of the variables 

(z-j, . . ., z 1K ). Also, let the subscript j denote the j term in the regression model 

so that x.. = f.(z n , . . . , z.^) describes the transformations of the (z. 1? . . . , z.^) to 
13 j ±± liv th th 11 

produce the value of x.j for the j term at the i observation. 

The regression model can now be rewritten as 

y i = b l x il + b 2 x i2 + • • • + ViJ + e i i = 1, ■ • • , N (3) 

where e^ denotes the difference between the observed value and the expected value of 
y^. For the N observations, it is convenient to write this regression model in matrix 
notation as y = Xb + e where 
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More often than not, the analyst feels the model 

y { = b 0 + b i x u + • • • + V^j + e i i = 1, . . . , N (5) 

is more appropriate. Let a Q = b Q + bjX ^ + . . . + bjX j. Then, as a result of adding 
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this equation to, and subtracting it from, equation (5) and rearranging terms 

y i =( b 0 + b ix > i + . . . +bjx j) 

+ b l(Xii - x a ) + . . . + bj(Xij- x j) + e i i = l, . . ., N (6) 

If then, a dummy variable x Q is introduced such that for all values of i, x iQ = 1. 0, 
equation (6) may be written as 

y i = a 0 x i0 + b l^ x il " x . 1^ + • • • + b j( x ij ‘ x . j) + e i i = 1, . . . , N (6a) 

Equation (6a) now resembles equation (3) and may be written in matrix notation, similar 
to equation (4), as y = Xb + e where now 




e s 



*1J 




X 2J ' x . J 


' x . jy 


( 7 ) 
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Estimating b 


Equations (4) and (7) are similar in form and for N > J are an overdetermined set 

a 

of linear equations. There will be some vector b which is a "best'' vector to use. If 

O 

the vector e is composed of random variables ej such that E^) = 0, V(e i ) = a < +°°, 
and the e. are uncorrelated, then as is well known, the method of least squares gives 
the linear minimum variance estimators b for b. And b is given by 


b = (X’X)" 1 X'y 


( 8 ) 


The matrix X’X divided by N - 1 is called the variance-covariance matrix of the 

A 

independent variables. The variance -covariance matrix of b is given by 

V(b) = o 2 (X’X) _1 (9) 


It is important to note that when the form of equation (7) is used, X'X is 




X'X = 



N 

2>ii- 

i 


|)(x 


iJ 


This is seen to be symmetric and of the form 


t 

X'X = I 
\° 



0 

N 

i: *u - x .i )<x u 

l 


N / 
Z( x u- x .j ) 2 J 

( 10 ) 
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Hence, 



RAPIER uses this relation to advantage by storing only the upper triangular part of 
B and computing only the coefficients bj, . . bj by matrix manipulations. Then 
bp is given by the simple equation 


b o = y 


b l x . 1 ” b 2 x . 2 


b J x .J 


( 11 ) 


where y = = ag. It can also be shown that 

V(b Q ) = V(y) + V(b'x) = jl + x’(X*X ) _1 xja 2 


COV(b Q ,b) = -(X’X )" 1 xcr 2 


When there is no bp term in the regression model, 

’ N « N 

E x u E x n x i2 • • 

1 1 


N 

E x n x tJ 


1 


X'X = 


( 12 ) 


N N N „ 

E X U X U E x 12 x U ' ' • E X U 

1 1 1 


Comparing this to equation (10) shows this form of X'X to be similar to the lower right 
submatrix in equation (10). This similarity is used to simplify notation by assuming 
that X'X represents either the form of equation (12) or the lower right portion of equa- 
tion (10) and considering the calculation of bg as a special case. Thus, further refer- 
ence to b implies 
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I I III 


I 1 1 II II ■■■>-! 



Correlation Matrix 


Another matrix of interest both computationally and statistically is the correlation 
matrix C. The elements of C, which are denoted C^, are the sample correlation co- 
efficients between the terms X. and Xj. These are 


N 



and all these numbers are between 1. 0 and -1.0. 

The calculation of C can be expressed in matrix notation conveniently by defining a 
diagonal matrix S = diagfej, s 2 , . . . , Sj) with elements 



(14) 


Then 


C = S(X’X)S 


(15) 


and 


(X’X) 1 = sjs^cx’x) 1 S -3 JS = SC _1 S (16) 

The algorithm of RAPIER performs the following operations: (1) constructs the 
X’X matrix, (2) computes C, (3) inverts C, (4) computes (X’X)" 1 from C" 1 by equa- 
tion (16), and (5) computes the b estimates. Because C is a normalized matrix, the 
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inversion of C is likely to be more accurate than direct inversion of X'X. Examination 
of the structure of X’X and/or C is of assistance in evaluating the possible numerical 
problems. 

It may also be that the independent variables are random variables. Then X'X 
divided by N - 1 represents the variance-covariance matrix and C the sample corre- 
lation matrix. If the independent variables are considered to be from a multivariate 
distribution, it is useful in some cases to consider the eigenvalues and eigenvectors of 
X'X and/or C. 

For these reasons, RAPIER includes options to compute and print these quantities. 
As a partial check on the accuracy of the inversion process, it is also possible to have 
C- C * computed and printed. This should be the identity matrix. 


Estimating o 2 

For any regression model y = Xb + e, there are possibly two methods of estimat- 
ing a . First, if the assumed regression model is, in reality, the true model, it is 
well known that an unbiased estimator is given by 

J, _ y'y - b'x'y 

ct RES(J) - ^TF 

= SSQ(RES) 

N - J - D 


= MS(RES(J)) (17) 

o 

Second, where there are replicated data points, another estimator of a , depending only 

on V(e.) = a for all i and not on the correctness of the assumed model, is the pooled 

mean squares computed from the replicated data points. 

Assume the observations are grouped into replicate sets in sequence. Let R be 

f h 

the number of sets of replicates and r i be the number of replicates in the i repli- 
cate set. Let 


r*+r^ 

SSQ(1) = Y (y n - y,) 2 
n=r*+l 


( 18 ) 
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where 


i-1. 

r * = E r i 

j=i 


It is assumed y is from the i th replicate set and y. is calculated only from those 
iu n. 1 

y n in the i replicate set. Then define the pooled sum of squares due to replication as 

R R 

SSQ(REP) = 2 SSQ(i) and the pooled degrees of freedom as NPDEG = £ (r, - 1). The 
i=l i=l 


second estimator of cr becomes 


2 _ SSQ(REP) 

REP NPDEG 

= MS (REP) (19) 

It can be shown (ref. 2) that the sums of squares due to residuals can be partitioned 
into a component due to replication and a component due to lack of fit; that is, 

SSQ(RES) = SSQ(LOF) + SSQ(REP) (20) 

o 

This partitioning is used later to determine the estimate of a to use in tests of hy- 
potheses. 


HYPOTHESIS TESTING 
Test NE - Normality of e 

As stated before, the only assumption necessary for b to be a linear minimum 

o 

variance estimator is that E(e.) = 0. 0, V(e.) = a < +°°, and e. be uncorrelated. If it 

1 9 1 * 

can further be assumed that e^ ~ N(0, & ), a number of standard tests become available. 
RAPIER computes a chi-square statistic and the sample skewness and kurtosis for test- 
ing this hypothesis. 

2 

Under the hypothesis e. ~ N(0,a ), the studentized resiudals defined by 
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will be distributed as Student’s t with the degrees of freedom associated with the esti- 
mate a. If the degree of freedom is 30 or more, the t distribution is very close to the 
normal. 

The range of possible studentized residuals is (-°°, +°°) and may be divided into k 
segments or cells each with probability p., so that each segment will have Np^ as the 
expected number of observations falling into it. Let n. denote the number of student- 
ized residuals in the 1 cell. Then a chi-square goodness-of-fit statistic may be cal- 
culated as 



i=l 


RAPIER computes this statistic by using an even number of cells greater than or equal to 
four and less than or equal to 20, such that the expected numbers of observations per cell 
is five or more. The bounding values for the i in cell are Z where F(Z.) = (i- k)/N 
and F(Z) is the cumulative normal distribution function. Then each cell has the same 
expected number of observations, say f = N/k. Then 



i=l 


*2 

This statistic is not computed for less than 30 degrees of freedom for the estimate a . 

Two other statistics which may be used to test the normality of an empirical distri- 
bution are skewness and kurtosis. Define the moments about the origin as 


m 2 


m 4 


1 _ 

N 

1 _ 

N 

1 _ 

N 


E z ? 
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th 9 / 9 

where Z. is the i studentized deviate. Then skewness is RELSKW = mT./m'z, which 
1 / 2 o' ^ 

should be nearly zero. Kurtosis is RELKUR = which should be nearly 3. 

Probability points for these are tabulated in reference 6. 

If these statistics indicate nonnormality, there are three possible courses of action. 

First, perhaps a transformation of the response variable or the independent variables 

can be found which will bring the distribution of residuals closer to normal. RAPIER 

makes this task quite easy. Second, a different candidate model might be used (see 

ref. 2, ’’Analysis of Residuals"). As the last choice, it is possible to do nothing and 

simply rely on the robustness of the tests involved. See reference 4 for definition and 

discussion of robustness. 


Also note that the individual observations may be weighted to perform a weighted 
regression analysis. RAPIER permits the use of weights (ref. 2). In this case, the 
X’X and X’y matrices take the form 
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Analysis of Variance Table 


For most hypothesis testing of the regression model, it is convenient to summarize 
the available information in an Analysis of Variance (ANOVA) table, as follows: 


Source 

Sums of squares 

Degrees of 
freedom 

Mean squares 

Regression 

(REG) 

SSQfREG) = b'X'y - S c a 

j 

MS(REG) = SSQ(REG)/J 

Residual 

(RES) 

SSQ(RES) = y'y - b’X'y 

N - J - D b 

MS(RES) = SSQ(RES)/(N - J - D) 

Total 

SSQ(TOT) = y'y - S Q 

N - D 



j 0 if no bg coefficient is estimated. 

if a bg coefficient is estimated. 

0 if no bg coefficient is estimated. 

1 if bg is estimated. 

If there are replicated data points, another ANOVA table can be constructed to show the 
separation of the residual sums of squares into components from lack of fit and replica- 
tion, as in the following table: 



Source 

Sums of squares 

Degrees of 
freedom 

Mean squares 

Lack of fit 
(LOF) 
Replication 
(REP) 

SSQ(LOF) = SSQ(RES) - SSQ(REP) 
SSQ(REP) 

N - J - D - NPDEG 
NPDEG 

MS (LOF) = SSQ(LOF)/(N - J - D - NPDEG) 
MS (REP) = SSQ(REP) /NPDEG 

Residual 

(RES) 

y'y - b'X'y 

N - J - D 
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o 

Choice of Estimator for o 

o 

As mentioned previously, there are two possible methods of estimating a depend- 
ing on whether there are replicated data points. This is true for any given model equa- 
tion. When the backward rejection option of RAPIER is used, there is no longer one hy- 

o 

pothetical model but a series of models. Thus, there is the choice of estimator for a 
to be made after each rejection of a term in the previous model. 

As an example, consider the model 

y = b Q + bjXj + bgXg + bgXg + e (22) 


with replicated data points. The first step is to estimate bg, bp bg, and bg. There 
will then be the estimators 0 res(J) anc * ^REP‘ ^ the mo< ^ e ^ * n equation (22) has not 
left out any significant terms, both estimators are valid. 

The ratio F = MS(LOF)/MS(REP) can be used to test the hypothesis that there is no 
lack of fit, where F ~ F a b with a = N- J- D- NPDEG and b = NPDEG degrees of 

freedom. If the test accepts the hypothesis of no lack of fit, MS(RES) is a pooled esti- 

2 *2 
mate of a with more degrees of freedom and will usually make tests using 

more sensitive than those using a pEp . But there is the possibility that the hypothesis 

was accepted as a result of random fluctuation when there really is some lack of fit; 

that is, there is the possibility that is a biased estimator. If lack of fit is not 

concluded to be significant, the decision to pool or not is usually made on the basis of 

the number of degrees of freedom for replication. If this is "large" (no definition of 

large is given herein), a pEp is used. If "small, " the pooled estimate is 

used. 

In testing equation (22), should it be decided that bg is not significantly different 
from zero (see section Test TT - t-Tests), the coefficients of the following model would 
be estimated: 


y = bg + bjXj + bgXg + e 

From this model there is an estimate This estimate could also be biased 

since bg may be small but nonzero and the decision of bg = 0 may have been due to 
random fluctuation. 

At the first step, the lack of fit can be considered a random sample of an infinite 
possibility of biases. But the biases due to pooling mean squares after rejecting terms 
can be considered to be systematic biasing and hence less desirable. 
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RAPIER provides three strategies of pooling estimates for use in the decision pro- 
cedure: 

(1) Never pool. This is appropriate only when there are replicated data points. The 
estimator used in all t-tests is <* RE p- 

(2) Always pool initial residual. This will always pool the lack of fit and replication 
from the first model only. Additional mean squares due to rejected terms will be 
ignored. 

(3) Always pool. This strategy will always use ^ EE g(j_j) for the model with i re- 
jected terms. 

A rule for pooling lack of fit and replication mean squares is discussed by Draper 
and Smith (ref. 2). Related work as applied to factorial designs is presented by Holms 
(ref. 7) and Bozivich, Bancroft, and Hartley (ref. 8). 


Test OR - Overall Regression 

One of the first tests usually applied to a regression model is the test of the overall 
significance of the model. In the notation of hypothesis testing this is stated Hq! b = 0; 
Hj: b * 0 where 


b = 



a o 

The statistic for this test is F = MS(REG)/ct . Then F ~ F„ . with a = J - D, and 

b equals the degrees of freedom associated with cr . 

Another useful statistic for judging the significance of overall regression is 

R^ = SSQ(REG)/SSQ(TOT). The sampling distribution of R does not lend itself to very 

2 

simple tests except in the case of H rt : R = 0. The main value of R is that it must be 

u 2 

a number in the range 0 to 1 and 100 R is a measure of the percentage of variation in 
the y values that is accounted for by the regression model. 


Test SF - Sequential F-Test 


There is often reason to consider a partitioned form of 
hypotheses 


b’ = 



for testing the 
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Hq : w 2 “ ^pfl’ ' ' •’ b j) _ ® 

Hji w£ * 0 

Partition the matrix X corresponding to the partitioning of b and denote it as 
X = (WpW 2 ) where Wj is Nxp and W g is N x (J - p). Then the test statistic is 
F = (SSQ(REG)^/p)CT REg ^, where SSQ(REG)^ = wJWJy and Wj = (WJW^ -1 WJy. Then 
F ~ F & b with a = p and b = N - p - D. Sometimes this test is performed with p = 1, 
p = 2, . . . , p = J. This is then referred to (ref. 2) as the sequential F-test. RAPIER 
computes regressions for p = 1, . . . , p = J upon request. 


Test TT - t-Tests 

In many cases, the regression model contains terms whose estimated coefficients 
are ’’small. ” This may be an indication that the term does not have a real effect on the 
dependent variable and that the coefficient is nonzero due to random sampling variation. 
If this is true, it is desirable to delete the term from the regression model. A test 
statistic for deciding this is 


b i 

t = i (23) 

5 2 (x'x)j 1 1 

— 1 f ft __ i 

where (X'X)T denotes the i diagonal element of the (X’X) - matrix. The statistic 
t ~ tjj_j jy An equivalent test statistic is 


b a 

F = t = i (24) 

(X’X)^ 

where F ~ Fj This is often referred to (ref. 2) as the partial F-test. The 

quantity b i/[^ ^)i/] is callec ^ the sum squares due to b^ if x. were last to enter 
the equation. RAPIER computes and prints the t-statistics and the probability associated 
with the interval (-t,t). 

This particular test is the basis for the rejection option of RAPIER. The analyst 

A 2 

has chosen which a estimator to use by the choice of strategy. Then the analyst may 
choose a significance level which all coefficients must meet. For example, suppose a 
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significance level of 0. 900 is chosen. The t- statistic is then computed for each coeffi- 
cient, and the coefficient with minimum |t| is identified. If min|t| > q all 

terms are concluded to be significant at the 0. 900 (or 90. 0 percent) level of significance. 
If min|t| < t N _j_ D q 950 , the term corresponding to the minimum |t| is dropped from 
the hypothetical model, and the regression is recomputed. This process is repeated 
until all remaining coefficients are significant at the specified level of probability. 


PREDICTING VALUES FROM ESTIMATED REGRESSION EQUATION 


Regression equations are often used to predict an estimated response at some condi- 
tion of the independent variables. Useful estimates of parameters to know are the vari- 
ance of the regression equation and the variance of a single further observation at the 
desired combination of the independent variables. 

Let x’ = (xj, . . . , Xj) denote the vector of independent variables at which a pre- 


diction is desired. Let x* = x - x. Let a „ denote the estimated variance of the 

-2 M' x 

regression equation at x. Let <x x denote the estimated variance of a single further 
observation at x. Then, 


; m- x = S2 [n + x *' (X,X> ’ 1x *] 

<7y >x = a 2 |l.O + ^ + x^X’X)" 1 x*j 


(25) 

(26) 


where, as before, D = 1 if a bp coefficient is estimated and D = 0 if a bp coeffi- 
cient is not estimated. The quantity s = * s called the standard error of esti- 


mate and often is used as a simple approximation to cr, 
if N is very large and x = x, in which case, 


y-x 


This approximation is close 


WK 1 - 0 ^) 


» s 


When x ^ x, this may be a poor approximation. RAPIER accepts input vectors x and 
computes y = b Q + bjXj + . . . + bjXj, as well as <r 2 -x , ct^ > x , <Jy. x , £ y . x , and the 
standard error of estimate. 
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USER'S GUIDE TO INPUT 
Sample Regression Problem 


Let 

Xj temperature 

Xg time 

Xg pressure 

yj output, lb 

Yg cost of operation 

The data are coded into standardized units, as is often done in experimental design anal- 
ysis. The y 1 variable is assumed to be of primary interest in this problem. 

Table I contains the x and y data. Table n contains a summary of the type of in- 


TABLE I. - x AND y DATA 


Group 

x i 

x 2 

X 3 


?2 

Group 

x i 

x 2 

x 3 

*1 

9.61 
10.01 
10. 12 
9.95 

*2 

1 

-1 

-1 

-1 

9.17 

38.5 

8 

0 

0 

0 

49.1 
49.3 

50.1 
51.8 

2 

1 

1 

-1 

-1 

-1 

-1 

12.76 

12.97 

43.1 

44.0 







3 

-1 

-1 

1 

1 

-1 

-1 

9.11 

8.96 

58.3 

58.7 




9 

-2 

0 

0 

11.78 

50.1 

10 

2 

2 

0 

0 

0 

0 

23.83 

22.90 

57.6 

58.1 

4 

1 

1 

-1 

17.03 

63.2 

5 

-1 

-1 

-1 

-1 

1 

i 1 

9.05 

8.86 

38.7 

39.6 

11 

0 

-2 

0 

7.99 

28.6 

12 

0 

0 

2 

2 

0 

0 

12.11 

11.70 

71.0 

70.2 

6 

1 

, 1 

-1 

-1 

1 

1 

12.60 

13.21 

44. 1 
43.8 

13 

0 

0 

0 

0 

-2 

-2 

10. 11 
10.01 

48.7 

49.9 

7 

1 

1 

1 

1 

1 

1 

17.20 

17.04 

62. 1 
62.8 

14 

0 

0 

2 

10.02 

50.8 


TABLE n. - FUNCTIONS OF INPUT CARDS 


Type of input 

Function 

1 

Identification 

2 

Definition of problem size 

3 

Definition of problem logic 

4 

Terms, transformations, and constants 

5 

Control rejection option 

6 

Provide replication information 

7 

Data input unit and format 

8 

Data 

9 

Prediction information 
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Figure L - SwmIi Input form. 







put cards and their basic functions. Figure 1 shows a sample set of data for a complete 
regression as it would be written on a FORTRAN coding sheet. 

Use the model equations 

= b 0 + b l x l + b 2 x 2 + b 3 x 3 + b 12 x l x 2 + b 13 x X x 3 + b 23 x 2 x 3 

2 2 2 

+ bjjXj + b 22 x 2 + b 33 x 3 + e i = 1, 2 (27) 

Test whether the interaction terms as a group are significant, given that the linear terms 
are in the model. Predict the response at the point (x^, x 2 , x^) = (-1, +1, +1). 


Types of Input Cards 

Nine types of data card may be used to define a regression analysis for RAPIER. A 
summary of the types and their functions is presented in table n. 

Type 1. - In type 1 input, as many as 100 cards with Hollerith information may be 
read to identify and describe the problem. At least one card is read. The first two 
columns of the first card used to specify the additional number of identification cards to 
be read, and columns 3 to 80 are used for Hollerith information. Each following card 
uses columns 1 to 78. (See lines 1 to 16 in fig. 1.) 

Type 2. - In type 2 input, one card with three four-column fields followed by a five- 
column field specifies 

(1) Number of independent variables to read 

(2) Number of dependent variables to read 

(3) Number of terms in the model equation (not counting bg) 

(4) Number of observations 
(See line 17 in fig. 1. ) 

Type 3, - In type 3 input, one card with 10 one-column fields specifies 

(1) The bg term in the model equation (T or F) 

(2) Computation of t-statistics and their confidence levels (T or F) 

(3) Weighting factor either of 1. 0 (T) or supplied with each data point (F) 

(4) Computation of residuals and chi-square test (T or F) 

(5) Computation of eigenvalues and eigenvectors of correlation matrix (T or F) 

(6) Computation of eigenvalues and eigenvectors of X'X (T or F) 

(7) Computation of product of correlation matrix and its inverse (T or F) 

(8) Use of bordering inversion technique for computation of sequential regression 

(T or F); see Test SF - Sequential F-Test 


22 



I 


(9) Use of an economy version of output which does not print the matrices X'X, 

(X’X)”*, x’y, C, or C - * (T or F) (If item 7 of this set is set T, then C * C - * 

is printed. ) 

(10) The pooling strategy: 

(1) Never pool. Always use replication error. (If there is no replication, the 

program sets this to 3.) 

(2) Pool initial residual. 

(3) Pool all residuals. 

(See line 18 in fig. 1.) 

Type 4, typ e 4A, type 4B, and type 4C. - Type 4 card has two four-column fields 
specifying 

(1) Number of transformations 

(2) Number of constants 
(See line 19 of fig. 1.) 

If the number of transformations is zero , and t herefore the number of constants is 
zero, the type 4A, type 4B, and type 4C cards are not expected by the program . In this 
case the program assumes the independent and the dependent variables are arranged on 
the input cards as 


x l’ x 2’ • • •’ X J’ y l’ • ' •’ y NODEP 

where NODEP is the number of dependent variables. 

If a weighting factor other than 1.0 is to be used (i. e. , if item 3 of the type 3 card 
contains an F), the value of the weighting factor for each observation must appear as the 
last item in the list, so that in this case the data for one observation is entered on the 
cards as 


x l’ x 2’ - • ' ’ X J’ y i> • • • > y NODEP ,WT 

For each observation, RAPIER reads a total of M numbers, where M is the sum 
of the number of independent and dependent variables. These numbers are stored con- 
secutively in an array called VAR, beginning with location 01 and ending with location M. 
If the weighting factor is not identically 1. 0, then M + 1 numbers are read, but the last 
number, being the weighting factor, is treated and stored separately. The data in VAR 
are used with appropriate weighting factors to cumulatively create X’X and X’y as 
shown in equations (21). 

When a more complex model is desired, information must be supplied instructing 
the program as to (1) where to find the values for the TERMS of the equation, (2) how to 


23 



create the TERMS from the variables and constants, and (3) what the values are for the 
constant terms. This can be achieved easily by use of the type 4A(TERMS), type 
4B(TRANSFORMATIONS), and type 4C(CONSTANTS) control cards (see p. 29). These 
three types and their functions can best be described by considering the sample model 
given by equation (27) as an example which illustrates their application. 

An array called CON has a twofold purpose. First, if the number of constants 
designated in the second field of the type 4 card is nonzero, that number of constants 
will be read from the type 4C card and stored consecutively in this array beginning with 
location 01. If the number of constants i s zero, the type 4C card is not expected by the 
program . Second, all the intermediate and final results of transformations are also 
stored in the CON array as the program obeys the instructions of the type 4B cards. 

The type 4A card must identify the relative location in the CON array where the value 
for each TERM is to be found for constructing the X'X and X'y matrices. 

The VAR and CON arrays for this example are illustrated in figure 2. Five numbers 
are read for each observation: Xj, Xg, x^, y^, and yg. These numbers automatically 
enter the VAR array beginning with location 01. Using transformation codes packed in 
fields of eight columns each on the type 4B cards, the program stores the result of each 
transformation into the appropriate relative location in the CON array as designated by 


Location VAR 


Location CON 


01 

02 

03 

04 

05 

06 


56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 


95 

96 

97 

98 

99 



Input constants 
followed by terms 
for regression 
equation 


Data stored are 
stored in same 
relative location 
in both arrays 
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Figure 2. - Map of VAR and CON arrays. Data transferred into any location of CON 
array beyond location 60 are immediately duplicated in same relative location in 
VAR array. 






the last two digits of the field. Each transformation code is made up of four subfields 
of two card columns each, with the following interpretation: 


Subfield 

Interpretation 

1 

Operand VI 

Relative location in VAR 

2 

Operation OP 

Arithmetic operation 

3 

Operator Cl 

Relative location in CON 

4 

Result CS 

Relative location in CON 


Thus, subfield 1 always references the VAR array and subfields 3 and 4 reference 
the CON array. The result of every transformation is a term which is stored in the 
designated location of the CON array, with the added feature that if the term is stored 
in relative location 61 or beyond, it is also stored in the VAR array. This is illustrated 
by the arrows in figure 2. This feature allows successive transformations to be per- 
formed more easily. The OP (operation codes) are tabulated in table III. 


TABLE in. - OPERATIONS* 1 AND CODE NUMBERS 


Operation 

code 

(OP) 

Resulting operation 

Operation 

code 

(OP) 

Resulting operation 

00 

No operation 

16 

1. 0/SQRT(VAR) 

01 

VAR + CONST 

17 

CONST **VAR 

02 

VAR*CONST 

18 

10. 0**VAR 

03 

CONST /VAR 

19 

SINH(VAR) 

04 

EXP(VAR) 

20 

COSH (VAR) 

05 

VAR**CONST 

21 

(1. 0-COS(VAR))/2. 0 

06 

A LOG (VAR) 

22 

ATAN(VAR) 

07 

ALOG10(VAR) 

23 

AT AN2 (VAR/C ONST) 

08 

SIN (VAR) 

24 

VAR**2 

09 

COS(VAR) 

25 

VAR**3 

10 

SIN (jt *CONST *VAR) 

26 

ARCSIN (SQRT (VAR)) 

11 

COS(ix*CONST*VAR) 

27 

2. 0*17 *VAR 

12 

1. 0/VAR 

28 

1. 0/(2. 0*!r*VAR) 

13 

EXP(CONST/VAR) 

29 

ERF(VAR) 

14 

15 

EXP(CONST/VAR**2) 
SQRT (VAR) 

30 

GAMMA (VAR) 


a All function names and operations are consistent with 
FORTRAN IV mathematical subroutines. 


25 



TABLE IV. - SEQUENCE OF TRANSFORMATIONS 


Transformation 

VI 

OP 

Cl 

CS 

Interpretation 

number 






1 

01 

00 

00 

11 

Xj - CON(ll) 

2 

01 

00 

00 

61 

- VAR(61), CON (61) 

3 

02 

00 

00 

12 

x 2 - CON(12) 

4 

02 

00 

00 

62 

x 2 - VAR(62), CON (62) 

5 

03 

00 

00 

13 

x 3 - CON(13) 

6 

03 

00 

00 

63 

Xg - VAR(63), CON (6 3) 

7 

61 

02 

61 

17 

- CON(17) 

8 

62 

02 

62 

18 

x 2 - CON(18) 

9 

63 

02 

63 

19 

Xg - CON(19) 

10 

61 

02 

62 

14 

XjX 2 - CON(14) 

11 

61 

02 

63 

15 

XjXg - CON(15) 

12 

62 

02 

63 

16 

x 2 Xg - CON(16) 

13 

04 

00 

00 

20 

y x - CON(20) 

14 

05 

00 

00 

21 

y 2 - CON(21) 


Table IV shows the sequence of transformations used to construct the terms of the 
example in equation (27). (See lines 21 and 22 of fig. 1. ) 

The arrays VAR and CON are shown in figure 3 both before and after one set of 
transformations performed on an observation. Note that CON now contains Xj, Xg, x^, 
x l x 2’ x l x 3’ X 2 X 3’ x l» x 2’ x 3’ y l» y 2 along with unused locations. It may be that not 
all of these quantities are needed to express the model equation. The type 4A(TERMS) 
card must be used to supply the locations of CON which contain the terms of the model 
equation. (Note that this allows the user to somewhat arbitrarily assign terms to loca- 
tions in CON. ) The terms identifying the independent variables must be first, and the 
terms identifying the dependent variables last, just as in the assumed convention when 
no transformations are performed. Thus, in this case the terms needed are, according 
to equation (27), 


X 1 x 2 x 3 x l x 2 x l x 3 X 2 X 3 

4 

x 2 

X 3 

y l 

y 2 

11 12 13 14 15 16 

17 

18 

19 

20 

21 


(See line 20 of fig. 1.) 

After the set of transformations has been performed on an observation, the contents 
of the relative locations of the CON array specified on the terms card are transferred 
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back to VAR in consecutive locations beginning with location 01. Thus, for this example, 
after selection of proper terms, the VAR array contains 


x v x 2 , x 3 , Xl x 2 , Xl x 3 , X2 X 3 , x\, x| y 1 , y 2 


in consecutive locations as required by equation (27) and the convention on independent 
and dependent variables. The bp term is accounted for by setting the proper item of 
type 3 input to T. 

There are three important facts to note concerning types 4A, 4B, and 4C. First, 
the transformation with OP = 00 is an identity transformation which simply transfers 
data from VAR to CON. If transformations are desired at all, a minimum requirement 
is that all variables at least be moved to CON so that the selection of terms will have a 
number to move back to VAR. Second, constants used in the transformations are stored 
in CON in locations beginning with 01 in sequence. Suppose there are NC constants ini- 
tially supplied. Then if a transformation has any of the loc ations 01 through NC refer- 
enced in subfield 4, a term will replace the constant . Third, the use of TRANSFOR- 
MATIONS and TERMS overrides the convention that independent variables must precede 
dependent variables on the input data cards. The convention holds true for the terms 
card data in this case. 

The sequence of input and formats is as follows: 

(1) Type 4A(TERMS): One or more cards, as necessary, with two-column fields 
denoting the relative locations of the CON array containing the final terms to be used in 
regression model (Up to 60 independent and nine dependent terms may be supplied. See 
line 20 of fig. 1.) 

(2) Type 4B (TRANSFORMATIONS): As many cards as necessary, with 10 trans- 
formations per card, each transformation being composed of four two-column subfields 
(A maximum of 100 transformations may be performed. See lines 21 and 22 of fig. 1. ) 

(3) Type 4C (CONSTANTS): As many cards as necessary, containing the required 
number of constants in (5E15. 7) format. As many as 60 constants may be supplied. 

A second example is shown to illustrate the flexibility afforded by the TERMS, 
TRANSFORMATIONS, and CONSTANTS. Consider a model with two independent and 
two dependent variables with four terms given by 

b 0 + b l x l + b 2 x 2 + b 3 x l x 2 + b 4 x l = y l = ~ 

y 2 
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where k = 1.0. Then a sequence of transformations which could be used is 


VI 

OP 

Cl 

CS 

Interpretation 

01 

00 

00 

61 

Xj - CON(61), VAR(61) 

01 

02 

61 

62 

Xj - CON(62), VAR(62) 

61 

02 

62 

63 

xj 1 - CON(63), VAR(63) 

01 

00 

00 

14 

Xj - CON(14) 

02 

00 

00 

02 

x 2 - CON(02) 

02 

02 

61 

12 

XjX 2 - CON(12) 

03 

00 

00 

11 

y t - CON(ll) 

04 

12 

01 

09 

— - CON(09) 
?2 


The TERMS information should then be 


x l x 2 


*1 


k_ 


14 02 


12 63 11 09 


The total process is illustrated by figure 4 and the type 4 input given in figure 5. 

Type 5 . - In type 5 input, one card with one one-column field and one three-column 
field specifies 

(1) Use of t- statistics to reject insignificant terms (T or F) 

(2) Probability level that t- statistic must meet to be considered significant (This is 
written without a decimal point; e. g. , 95-percent significance level is supplied as 950, 

99. 9 percent as 999, etc. See line 23 of fig. 1. ) 

Type 6 . - In type 6 input, one card with one one-column field specifies that the data 
contain replicated points (T or F). If there are replicated data points, as many cards as 
are needed are read, containing 20 four-column fields specifying 

(1) The number of replicate sets 

(2) The number of replicates in each replicate set 

Note that it is not safe for the program to assume that all data points with the same levels 
of the independent variables are true replicates. For this reason, the user must arrange 
replicate sets. RAPIER does check that all independent terms are the same within a 
replicate set. If not, the program stops. A nonreplicated data point is considered to be 
a group of size 1. Note that the data in table I are grouped to clearly indicate the repli- 
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cated data points. There are 14 such groups. Thus, the first field of the second type 6 
card contains a 14, and the remaining fields contain the count of the replicates in each 
group. (See lines 24 and 25 of fig. 1.) 

Type 7. - In type 7 input, one card with one two- column field specifies the input unit 
number the data will be on. The remainder of the card is used to supply the format in 
which the data will be supplied. Note that if a weighting factor other than 1.0 is to be 
used, it will be read with each data point, and the format must allow for this. The cur- 
rent example uses a weighting factor of 1. 0. 

The format is (5F6. 0) since there are three independent variables (x^,Xg, x^) and 
two dependent variables If a weighting factor other than 1. 0 is used, it must 

appear with every data point, and the format could, for example, be (5F6. 0, F10. 3). 

(See line 26 of fig. 1.) 

Type 8. - Type 8 input consists of the input variables. Each observation consisting 
of the given x’s and y’s is read by execution of one READ statement. Thus, there 
will be at least one card for each observation. If the transformation option is not used, 
the program expects the first variables read to be the independent variables and the last 
ones to be the dependent variables; that is, the data must be arranged as 

x i ,x 2’ . . . , Xj, y^yg, • • •> ^nodep 
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Otherwise, if transformations are used, appropriate use of the terms card information 
allows for more flexibility of input. (See lines 27 to 50 of fig. 1. ) 

Type 9. - In type 9 input, one card with one column is used to indicate if predictions 
are desired (T or F). (See line 51 of fig. 1.) If this is false, a new case is started. If 
it is true, the following cards are read: One card with one four-column field specifies 
the number of predictions desired. This is followed by cards with the values of the in- 
dependent variables at which predictions are desired. Only the final regression model 
is used, but the number of independent and dependent variables originally supplied on the 
type 2 data cards are read. All transformations in type 4 input are performed. Then 
the proper terms are chosen by the program to correspond to the final model. The de- 
pendent variables are not needed and, hence, may be left off the data card unless one 
of the transformations of the dependent variables might lead to an impossible operation 
(e.g. , log(y 2 )). (See lines 52 and 53 of fig. 1.) 



SAMPLE OUTPUT 


sample rapier problem 

MODEL EQUATION I- 1.2 

Yin » BO ♦ 8 1* XI *B2*X2 «-B3*X3 

*■ B 12* XI# X2 *B13*X1*X3 +B23*X2*X3 

♦ B11*X1**2 +B22*X2**2 +B33*X3**2 ERR 

XI - TEMP ** 

X2 * TIME ** DATA CODED TO STANDARIZED UNIT 

X 3 - PRESS ** 

VI - POUNDS OUTPUT 
Y 2 » COST OF OPERATION 

THE DATA IS FROM AN INCOMPLETE FACTORIAL OES IGN WITH OME 
REPLICATION IFICTITIOUS DATA) 

3 2 9 25 

T TTTFFFFF 2 

THERE IS A BO TO ESTIMATE 
N T ERMIK )*= 


11 12 13 1* 15 16 17 18 19 20 21 

THE TRANSFORMATIONS ARE 


1 

0 

0 

LI 

1 

0 

0 

61 

2 

0 

0 

12 

2 

0 

0 

62 

3 

0 

0 

13 

3 

0 

0 

63 

61 

2 

61 

17 

62 

2 

62 

ie 

63 

2 

63 

19 

61 

2 

62 

l* 

61 

2 

63 

1 5 

62 

2 

63 

16 

* 

0 

0 

20 

5 

0 

0 
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< 5F6-01 

SAMPLE RAPIER PROBLEM 

OBSERVED VARIABLES, WEIGHT = 1.000000 OBSERVATION = , l 

-l. 000000 -1.000000 -1.000000 9.170330 38.50300 

TERMS OF THE EQUATION, OBSERVATION * l 

-l. 000000 -l. 000000 -l. 000000 1.000330 L. 000300 l. 000000 l. 000000 1.003000 1.033333 

9.170000 38.50000 

*♦ REPLICATE SET 1 ******♦**************+**#********************+*****************+***************#******************** 

* ************** *«** *********** ******** **************************** *********** ******** ** + ******** *** ** * *** *** **♦*♦***<■* * ****** 

OBSERVED VARIABLES, WEIGHT = l. 000000 OBSERVATION * , 2 

1.000000 -l. 000000 -l. 000000 12.76330 *3.10300 

TERMS OF THE EQUATION, OBSERVATION = 2 

L. 000000 -1.000000 -1.000000 -1.000330 -1.000300 1.000000 l. 000000 1.000000 1.333333 

12.76000 *3.10000 

OBSERVED VARIA3LES, WEIGHT = l. 000000 
L. 000000 -l. 000000 -l. 000000 

TERMS OF THE EQUATION, OBSERVATION = 

1. 000000 -1.000000 -l. 000000 

12.97000 **.00000 

** REPLICATE SET 2 **♦**********************************************+************************** ********************* 

DEP. VAR . 1 SSQ = 0. 220* 895 E— 01 SUM= 25.733000 ME AN* 12.865000 

DEP . VAR. 2 SSQ= 0.*050598 SUM= 87.099999 MEAN* *3.550000 

** ********************* **** ****** ******* ********** ** **** * ** * **** * **+ ************************** * ****************************** 

OBSERVEO VARIA3LES, WEIGHT = 1.000000 OBSERVATION * , * 

-l. 000000 l. 000000 -1.000000 9.110300 58.30300 

TERMS OF THE EQUATION, OBSERVATION = * 

-I. 000000 I. 000000 -L. 000000 -1.033330 1.000300 -1.000000 1.000000 1.000000 1.333333 

9.110000 58.30000 

OBSERVEO VARIABLES, HEIGHT = l. 000000 OBSERVATION = , 5 

-l. 000000 I. 000000 -l.OOOOOO 8.963333 58.70300 

TERMS OF THE EQUATION, OBSERVATION = 5 

-1.000000 1.000000 -1.000000 -1.333333 1.003300 -1.000000 1.000000 1.000000 1.33)333 

B. 960000 58.70000 

** REPLICATE SET 3 *************************************** ******************************** **** ************** *********** 

DEP. VAR. 1 S SQ= 0. H251*5E-0l SUH= 18.073300 MEAN* 9.0350000 

DEP. VAR. 2 SSQ= 0. 79956 05E-01 SUM= 117.03000 MEAN* 58.500000 

** ** *** * * ** * •*** * titi** ** * **** ** * * ****** ** * ** * *** * + * * * **** * +** * *** * ** * ****** ** ** * + * * * * * * * * * ** * ***** * * * * **** * * * * ** * ** * * * * **** * * * 


OBSERVED VARIABLES, WEIGHT = 1.000000 OBSERVATION = , 6 

1.000000 1.000000 -1.000000 17.03030 63.20300 

TERMS OF THE EQUATION, OBSERVATION * 6 

1.000000 1.000000 -1.000000 1.003330 -1.003300 -1.000000 1.000000 1.000000 1.033333 

17.03000 63.20000 


OBSERVAT I ON = , 3 

12.97333 **.00300 

-1.033330 -1.003300 1.000000 1.000000 1.000000 1.333333 
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** REPLICATE SET 4 *************** ****** + *♦*+** ** * **** * *** ********** *** *** ******** * ***** * ****** ***♦ ****** *** 4 * ****** *** 

** ************* *********************** *********************** *** ************ ************************************************* 

OBSERVED VARIABLES, HEIGHT = I. OOOOOO OBSERVATION = , 7 

-L.OOOOOO -I. OOOOOO 1.000000 9.050330 38.70000 

TERMS OF THE EQUATION, OBSERVATION = 7 

-1.000000 -l. OOOOOO l. OOOOOO 1.000330 -1.003)00 -1.000000 1.000000 l. OOOOOO 1.033003 

9.050000 38.70000 


OBSERVED VARIA3LES, WEIGHT = l. OOOOOO 
- 1.000000 - 1.000000 1.000000 

TERMS 3F THE EQUATION, OBSERVATION = 8 

-1.000000 -l. OOOOOO l. OOOOOO 

8.860000 39.60000 

** REPLICATE SET 5 **************************************************************************************************** 

DEP. VAT • 1 SSQ= 0. 18051 15E-01 SUM= 17.913003 ME AN = 8.9549999 

DEP. VAR. 2 SSQ= 0.4050293 SUM= 78.299999 ME AN = 39.150000 

***************************************************************************************************************************** 

UBSERVED VARIABLES, WEIGHT = 1.000000 OBSERVATION = . 9 

1.000000 -L.OOOOOO 1 • OOOOOO 12.60330 44.10300 

TERMS OF THE EQUATION, OBSERVATION = 9 

1.000000 -1. OOOOOO l. OOOOOO -1.030333 1.003300 -l. OOOOOO 1.000000 1.000000 1.033333 

12.60000 44.10000 

OBSERVED VARIABLES, WEIGHT = 1.000000 

1.000000 -l. OOOOOO 1.000000 

TERMS OF THE EQUATION, OBSERVATION = 10 

1. 000000 -I. OOOOOO 1.000000 

13.21000 43.80000 

** REPLICATE SET 6 ****************************************************************** * * *t ***************** ************* 

DEP. VAR. i SSQ= 0.1860542 SUM= 25.8 13303 ME AN= 12.905000 

DEP. VAR. 2 SSQ- 0. 4501 343E-01 SUM= 87.903000 ME AN= 43.950000 

************ ********* *********************************** * ************ ******************************************************** 


OBSERVED VARIABLES, WEIGHT = 1.000000 OBSERVATION = , 11 

1.000000 1.000000 1.000000 17.20330 62.10300 

TERMS 3F THE EQUATION, OBSERVATION * 11 

1.000000 1.000000 l. OOOOOO 1.000330 1.003300 l. OOOOOO 1.000000 1.003000 1.003333 

17.20000 62.10000 

OBSERVED VARIABLES, WEIGHT - l. OOOOOO 

1.000000 1.000000 l. OOOOOO 

TERMS OF THE EQUATION, OBSERVATION ■= 12 

L.OOOOOO 1.000000 1.000000 

17.04000 62.80000 

** REPLICATE SET 7 **************************************************************************************************** 

DEP. VAR. 1 SSQ= 0. 1280212E-01 SUM= 34.243300 ME AN= 17.120000 

DEP. VAR. 2 SSQ = 0.2449951 SUM= L24. 93000 ME AN = 62.450000 

********************************************************************************************** ********** ********************* 


UBSERVED VARIABLES, WEIGHT = L.OOOOOO OBSERVATION = , L3 

000 9.610333 49.10300 

TERMS OF THE EQUATION, OBSERVATION - 13 

000000003 
9.610000 49.10000 

OBSERVED VARIABLES, WEIGHT * 1.000000 OBSERVATION = , 14 

000 10.01333 49.30300 

TEFMS 3F THE EQUATION, OBSERVATION = 14 

OOOOOOOOO 
10.01000 49.30000 

OBSERVED VARIABLES, WEIGHT - l. OOOOOO OBSERVATION = . 15 

000 10.12330 50.10300 

TERMS 3F THE EQUATION, OBSERVATION = 15 

OOOOOOOOO 
10.1200 0 50. 10000 

0 BSERVEO VARIABLES, WEIGHT = l. OOOOOO 
0 0 0 

TERMS OF THE EQUATION, OBSERVATION = 16 

0 0 0 

9.950000 51.80000 

** REPLICATE SET 8 ******** ********************************+************************** ******** ***►**+****♦********** *** 

DEP. VAT. I SSQ= 0.1450768 SUM= 39.693000 MEA'I= 9.9224999 

DEP. VAR. 2 SSQ= 4.5275879 SJM= 200.33000 MEAN= 50.074999 

************************************************************* **************************************************************** 

OBSERVED VARIABLES, WEIGHT = l. OOOOOO OBSERVATION = , 17 

-2.000000 0 0 Li. 78000 50.L0300 

TERMS OF THE EQUATION, OBSERVATION * 17 

-2.000000 0 0-0-0 0 4.000000 0 0 

11.78000 50.10000 

** REPLICATE SET 9 * + ************************************************************************************************* * 

**************!< ***************************************************************************************** ********************* 

OBSERVED VARIABLES, WEIGHT = 1.000000 OBSERVATION = , 18 

2. 000000 0 0 2 3.83330 57.60300 

TEFMS 3 F THE EQUATION, OBSERVATION = 18 

2.000000 0 0 0 3 0 4. OOOOOO 0 3 

23.83000 57.60000 


QBSERVAT l ON = , 16 

9.953330 51.80300 

0 3 0 0 0 0 


08SERVAT ION - , 12 

17.34333 62.83300 

1.033333 1.000300 1.000000 1.000000 1.003003 1.33333) 


OBSERVAT ION = , 10 

13.21330 43.83300 

-1.003333 1 .7)00300 -1.000000 1.000000 1.000033 1.33)33) 


OBSERVAT ION = * 8 

8.85333) 39.60300 

1.000330 -i. 003300 -1.000000 1.000000 1.003003 1.033333 
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OBSERVED VARIABLES, WEIGHT = l. OOOOOO 

2.000000 0 0 

TERMS OF THE EQUATION, OBSERVATION = 19 

2.000000 0 0 

22.90000 58.10000 


OBSERVAT ION * , 19 

22.90000 58.10000 

0 0 


0 


*.000000 0 


0 


** REPLICATE SET 10 ♦.*****•****♦*+**♦♦♦*****,.*** * 4 .** *******+*****+*«.****♦♦♦*******«.**•«.*►*#.***►**♦►***** + ** + *♦* ********* 

D6P. VAR. I $SQ= 0. *32*3*1 $UM= *6. 730000 MEAN- 23.365000 

OEP . VAR. 2 SSQ- 0.1250000 SUM= 115.70000 MEAN* 57.850000 

***************************************************************************************************************************** 


OBSERVED VARIA3LES, WEIGHT = 1.000000 OBSERVATION = , 20 


0 

TERMS OF THE 

-2.000000 0 
EQUATION, OBSERVATION « 

7.993330 

20 

28.60000 




0 

7.990000 

-2.000000 0 
28.60000 

-0 

0 

-0 

0 

*.000000 


** REPLICATE SET LL **> ************************************************** ************ *********** * **** ******************** 
** ***************************** ******* *********************** *************************** *********************** ************** 


OBSERVED VARIABLES, WEIGHT = 1.000000 

0 2.000000 0 
TERMS OF THE EQUATION, OBSERVATION = 

21 

12. 

OBSERVAT ION * , 
11330 71.00300 

21 




0 

12. 11000 

2.000000 

71.00000 

0 


0 

0 


0 

0 

*.003300 

OBSERVED VARIABLES, WEIGHT * 
0 2.000000 

l. OOOOOO 

0 


11. 

OBSERVAT ION * , 
70330 70.20300 

22 




TERMS OF THE 

EQUATION, OBSERVATION * 

22 







0 

11.70000 

2.000000 

70.20000 

0 


0 

0 


0 

0 

*. OOOOOO 


** REPLICATE SET 12 ***********************♦**********************+***+***+**+*********+*+********** ************** ****** 

OEP. VAR. 1 SSQ* 0. 8*0530*6-01 SUM= 23.B13303 MEAN* 11.905000 

DEP. VAR. 2 SSQ* 0-320190* SUM* 1*1.20000 MEAN* 70.599999 

********* *** *** *** * ******** ******* *********** ***************************** ****** ************ ************************** ******* 


OBSERVED VARIABLES, WEIGHT = l. OOOOOO 

0 0 - 2.000000 

TERMS 3 F THE EQUATION, OBSERVATION = 

0 0 - 2.000000 

10 . 1 1000 *8.70000 

QBSERVEO VARIABLES, WEIGHT = l. OOOOOO 

0 0 - 2.000000 

TERMS OF THE EQUATION, OBSERVATION = 

0 0 - 2.000000 

10.01000 *9.90000 


OBSERVATION = , 23 

LO. 11330 *8.70000 

23 

0 -0 -0 


OBSERVATION * , 2* 

10.01330 *9.90300 

2 * 

0 -0 -0 


0 0 *.333333 


0 0 *.033333 


** REPLICATE SET 13 **************************************************************************************************** 

OEP. VAR. I SSQ = 0. 50029 756-02 SUM* 20.123000 MEAN* 10.060000 

DEP. VAR. 2 SSQ = 0.7200317 SUM* 98.599999 MEAN* *9.300000 

*************** ************************************************* ****** ******************************************************* 


OBSERVED VARIABLES, WEIGHT = 1.000000 

0 0 2.000000 

TERMS 3 F THE EQUATION, OBSERVATION = 25 

0 0 2.000000 

10.02000 50.80000 


OBSERVAT ION = , 25 

10.02330 50.80300 

0 0 


0 


0 


0 


*.033303 


** REPLICATE SET 1* **************************************************************************************************** 

***************************************************************************************************************************** 


SUMS JF IN DEP AND OEP VARIABLES 


*.0000000 

0 

-2. 0000000 

0 

2.0000300 

-2.0000000 

2*. OOOOOO 

2*. 003000 

2*. 0000 00 

308. 10000 

1282.2000 







X TRANSPOSE 

X MATRIX 








ROW 

1 

2*. 00000 








ROW 

2 

0 

24.00000 







ROW 

3 

2.000000 

-2.000000 

24.00330 






ROW 

* 

-2.000000 

2.000000 

4.000030 

12.03000 





ROW 

5 

0 

*. OOOOOO 

2.000333 

-2.003000 

12.00000 




ROW 

6 

*.000000 

0 

-2.003333 

2.003000 

0 

12.00000 



ROW 

7 

10.00000 

-2. OOOOOO 

0 

3 

2.000000 

-2.000000 

60. 03030 


ROW 

8 

2 .OOOOOO 

6. OOOOOO 

0 

3 

2.000000 

-2.000000 

12.00000 

53. 03333 

ROW 

9 

2.000000 

-2. OOOOOO 

-8.000330 

3 

2.000000 

-2.000000 

12.03330 

12.03330 

ROW 

9 

60.00000 








X TRANSPOSE 

Y MATRIX 








ROW 

1 

127.5600 

260.5000 







ROW 

2 

22.36000 

238.5000 







ROW 

3 

- 12.2*000 

-U 0. 3000 







ROW 

* 

8.740000 

12.90000 







ROW 

5 

26.62000 

13 9. 7000 







ROW 

6 

-9.680000 

-95.90000 







ROW 

7 

382.0000 

1260. LOO 







ROW 

8 

275. 1600 

1276.100 







ROW 

9 

268.5200 

U 94.5 00 







MEANS OF IN OEP AND OEP VARIABLES 







C. 

1600000 

0 

-0. 8000000E 

-01 0 

3 . 8003000E 

-01 - 0. 8000000E- 

01 0.9600000 

0. 9503333 


0. 

9600000 

12.32*000 

51.288000 
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X TRANSPOSE X MATRIX WHERE X IS DEVIATION FROM MEAN 


ROW 

1 

23.36000 







ROW 

2 

-0 

24.00000 






ROW 

3 

2.320000 

-2.000000 

23.84000 





ROW 

4 

-2.000000 

2.000000 

4.000000 

12.03000 




ROW 

5 

-0.320000 

4.000000 

2.160333 

-2.003000 

11.84000 



ROW 

6 

4.320000 

0 

-2.160333 

2.003000 

0.160000 

11. 84000 


ROW 

7 

6.160000 

-2. OOOOOO 

1.923333 

-o 

0. 800000F- 0 1 

-0. 8000005-01 

36. 95330 

ROW 

8 

-1.840000 

6. OOOOOO 

1.920330 

-3 

0 . 0OOQOOE-O 1 

—0. 800000E-0 1 

-11.04000 

ROW 

9 

-1.840000 

-2.000000 

-6.083333 

-3 

0.800000E-01 

-0.8000005-01 

-11.04330 

ROW 

9 

36.96000 








X TRANSPOSE 

Y MATRIX WHERE 

X AND Y ARE DEVIATIONS FROM MEAN 

ROW 

1 

78.26400 

55.34800 

ROW 

2 

22.36000 

238.5000 

ROW 

3 

12.40800 

-7.723999 

ROW 

4 

8.740000 

12.90000 

ROW 

5 

L. 972000 

37.12400 

ROW 

6 

14.96800 

6.676000 

ROW 

7 

86.22400 

29.18800 

ROW 

8 

-20.61600 

45.18800 

ROW 

9 

-27.25600 

-36-41200 


CORRELATION COEFFICIENTS 


ROW 

ROW 

ROW 

ROW 

ROW 

ROW 

ROW 

ROW 

ROW 


l -OOOOOO 
-0 

0 • 9331028 —01 
-0.119455 
-0.192414E-01 
0.259760 
0.209642 
-0. 62620 3E -01 
-0.626203E-01 


l. OOOOOO 
-0. 3361256-01 
0. 117351 
0.237289 
0 

- 0. 671 5 19E-0 L 
0. 201456 
-0. 67I519E-01 


1.030300 
0.236492 
0.128566 
-0.128566 
0.6468 l 0E— 0 l 
0.64681 8E— 0 1 
-0.234826 


1.003000 
-3 .167789 
0.167709 
-0 
-0 
-0 


1.000000 

0. L35135E-01 1.000000 

0.382427E-02 - 0. 3 82 42 7E - 02 1.003330 

0.382427E-02 -0. 3 8242 7E - 02 -0.298731 

0. 38242 7E- 02 -0. 38242 7E-02 -0.298731 


ROW 9 1.000000 


SAMPLE RAPIER PROBLEM 

EACH COLUMN CONTAINS THE COEFFICIENTS FOR ONE DEPENDENT TERM 
CONSTANT TERM (BO) 

9.973995 50.L5123 

REGRESSION COEFFICIENTS (BI.....BK1 

1 2.940374 2.031365 

2 0.987377 10.16032 

3 — 0. 1383 3 2E- 01 0.370808 

4 1.068721 -0.467007 

5 0.820396E-01 -0.397416 

6 0.2037 6 7E- 01 -0. 191 401E-OI 

7 1.909437 0.965708 

8 0.910342E-02 - 0. 571182E-01 

9 0.3302 73E-01 -0.596032E-03 


ANOVA OF REGRESSION ON DEPENDENT VARIABLE l 
********* ****** **************************************************************** 

SOURCE SUMS OF SQUARES DEGREES OF FREEDOM MEAN SQJARES 


REGRESSION 425.390369 9 47.2655964 

RESIOUAL 1.08622742 15 0 . 7241 5 160E-0 l 


TOTAL 426.476597 24 

**************************************************************** *************** 

R SQUARED = SSQ(REG) / SSQITOT) = 0.997453 R = .998726 

STANOARD ERROR OF ESTIMATE 0.269101 

USING POOLING STRATEGY 2 THE ERROR MEAN SQUARE = 0 . 724 1 5 1 6 E-D 1 WITH DEGREES OF FREE OOM = 15 

F*MS(REG)/MS(ERR)= 652.70 COMPARE TO F( 9, 15) 


ANOVA OF LACK OF FIT 

**************************************************************** *************** 

SOURCE SUMS OF SQUARES DEGREES OF FREEDOM MEAN SQJARES 


LACK OF FIT 0.16945267 5 0. 33890533 E-01 

REPLICATION 0.91677475 13 0. 9 1677475 E-01 

RESIOUAL 1.08622742 15 0. 7241 5 160E-01 

F « MS( LOF) /MS(REPS) = 0.370 

******************************************************************************* 


36.96333 
-1 1. 34333 


1.333333 

-3.293731 
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ANOVA 3F REGRESSION ON 3c PENDENT VARIABLE 2 
*********** ******************************************************************** 

SOURCE SUMS OF SQUARES DEGREES OF FREEDOM MEAN SQJARES 


REGRESSION 2539.42480 9 282.150310 

RESIDUAL 10.4015808 15 0.69343872 


TOTAL 2 549. 8263 9 24 

**************************************************** *************************** 

R SQUARED * SSQIREG1 / SSO(TOT) * 0.995921 R = .997958 

STANDARD ERROR OF ESTIMATE 0.832730 

USING POOLING STRATEGY 2 THE ERROR MEAN SOUARE = 0.6934387 WITH DEGREES OF FREEDOM = 15 

F»MS(REG)/NS(ERR1» 406.90 C OH PARE TO FC 9, 15) 


ANGVA OF LACK OF FIT 

*************** **************************************************************** 

SOURCE SUMS OF SQUARES DEGREES OF FREEDOM MEAN SQJARES 

LACK OF FIT 3.52871704 5 0.70574340 

REPLICATION 6.87286377 ID D. 68728638 

RESIDUAL L 0. 401 5808 15 0.69343872 


F = M SI LOF l/MSIREPS) = 1.027 

********* *************************************************************** ******* 


SUMS OF SQUARES DUE TO EACH VARIABLE IF IT HERE LAST TO ENTER REGRESSION 


1 

164.5587 

78.54016 

2 

19.45477 

2 063.2 86 

3 

0.3499 7 2E- 02 

2.514697 

4 

10.44572 

2.001442 

5 

0. 66536 7E- 0 1 

1.561362 

6 

0.40 5068E- 02 

0.357697E-02 

7 

96.39906 

24.65779 

8 

0.22 10 86E- C2 

0. 870364E-01 

9 

0.28829 9E- 01 

0. 938933E-05 

STANOARD DEVIATION 

OF REGRESSION 

0 

0. 128439 

0.397453 

l 

0.6 168 1 8E- 01 

0. 190874 

2 

0.602399E- 01 

0. 186412 

3 

0.629246E- 01 

0. 194720 

4 

0.889836E- 01 

0.275359 

5 

0.8558706- 01 

0.264848 

6 

0.86 1556E-01 

0.2666 08 

7 

0.523340E- 01 

0. 161947 

8 

0.52 10 Q2£- 01 

0.161224 

9 

0. 52 344 0E— 01 

0.161978 


(DERIVED FROM DIAGONAL ELEMENTS DF (X TRANSPOSE XMNVERSE MATRIX ) 


(X TRANSPOSE X) INVERSE MATR 


ROH 

ROM 

ROW 

ROW 

ROW 

ROW 

ROW 

ROW 

ROW 


0 . 525393E- 01 
-0.5 1441 IE -02 
— 0.116754E-01 
0.19112 2E -01 
0 . 892031E - 02 
-0 ■ 247 18 7E - 01 
-0. 8929956-02 
0 .683955E- 03 
-0.2119306-02 


IX 

0. 5011 17E-01 
0. U5598E-01 
-0. 1 79948E-01 
-0. 222912E-01 
0. 73 0194E-02 
0. 149314E-02 
-0. 77161 7E -02 
0. 2 56242E -02 


0. 5467 79 E— 01 
-0. 28552 7E-01 
—0 .1934346-0 l 
0.1 939 19 E-0 1 
0. 19B443 E-02 
-0. 1803 10E-02 
0.915151E-02 


0. 109343 
0 . 307 362E-0 l 
-0.311 087E -0 1 
-D.355785E-02 
D . 262 664E-02 
-3.513 127E-02 


0. 101155 
-0. 133644E-01 
-0.22381 IE-02 
0. 296380E-02 
-0.39 75 2 9E -02 


0. 102503 
0.4717046-02 
-0.6197346-03 
0.38275 BE- 02 


0. 37321 V E -01 
3. 1525096-01 
0. I 503O5E-O1 


0.37VBV3E-31 
0. I 535 02E-31 


ROW 9 0.378359E-01 

SAMPLE RAPIER PROBLEM 


CALCULATED T STATISTICS 

THE I STATISTICS CAN BE USED TO TEST THE NET REGRESSION COEFFICIENTS Bill. 
47.67006 10.64246 

16.39073 54.54760 

0.219838 1.904315 

12.01031 1.698899 

0.958553 1.500541 

0.236510 0. 718214E-01 

36.48561 5.963113 

0.174729 0.354280 

0.630967 0.367971E-02 

UNOER NULL HYPOTHESIS THE INTERVAL I-T.T) 

.999. 


WHERE T IS GIVEN ABOVE, HAS APPROX PROBABILITY LISTED BELOW. 


MINUS 

SIGN INDICATES 

PROB EXCEEDS 

1 

-0.999 

-0.999 

2 

-0.999 

-0.999 

3 

0.1 71 

0. 924 

4 

-0.999 

0.890 

5 

0.647 

0. 845 

6 

0. I 84 

0.056 

7 

-0.999 

-0.999 

8 

0.135 

0.271 

9 

0.463 

0. 002 


THE DESIRED VALUE OF PROBABILITY IS 95.0 PERCENT 
THE TERM X( 8) IS BEING DELETED 

THE NUMBER OF DEPENDENT VARIABLES WAS 2 IT IS BEING SET TO ONE AND THE REJECTION OPTION EXERCISED ON 3EPEN36NT VARIABLE 1 
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CORRELATION COEFFICIENTS 


ROW 

1 

1.000000 






ROW 

2 

-0 

l. 000000 





ROW 

3 

0. 9831026-01 

- 0. 8361 2 5E-01 

1 • 000000 




ROW 

4 

-0.119455 

0.117861 

0.236492 

1.003000 



ROW 

5 

-0.192414E-01 

0.237289 

0. 128566 

-0.167789 

1.000000 


ROW 

6 

0.259760 

0 

-0.128565 

0. 167789 

0.1351356-01 

l. 000000 

ROW 

7 

0.209642 

-0.671519E-01 

0.646818 E-0 1 

-0 

0.382427F-02 

-0.382427E-02 

ROW 

a 

-0.626203E-01 

-0.67151 9E-01 

-0.204326 

-0 

0.382427E-02 

-0.382427E-02 


SAMPLE RAPIER PROBLEM 

EACH COLUMN CONTAINS THE COEFFICIENTS FOR ONE DEPENDENT TERM 
CONSTANT TERM ( BO) 

9.989932 

REGRESSION COEFFICIENTS <&!,«•. *8KJ 

1 2.940208 

2 0.989251 

3 -0. 1 3 374 76- OL 

4 1. 06808 3 

5 Q.813199E-01 

6 0.205272E-01 

7 1.905733 

9 0. 29372 2E- Oi 


1.000300 

-0.298731 


1.333333 


ANOVA OF REGRESSION ON DEPENDENT VARIABLE I 

******* ***** ******************************** ****** ********<.************ 4 .******* 

SOURCE SUMS OF SQUARES OEGREES OF FREEDOM MEAN SQUARES 


REGRESSION 425.388165 8 53.1735206 

RESIDUAL 1.08843231 16 0 . 6802 7023 E-0 1 


TOTAL 426.476597 24 

*********** ************ ***************************************** *************** 

R SQUARED = SSQIRESI / SSQCTOT) = 0.997448 R = .998723 

STANDARD ERROR OF ESTIMATE 0. 260320 

USING POOLING STRATEGY 2 THE ERROR MEAN SJUARE = 0. 724 1 51 6 E-0 1 WITH DEGREES OF FREEDOM = 15 

F-MSIR6G)/MSIERR) = 734.29 COMPARE TO F( 8, 161 


ANOVA OF LACK OF FIT 

*************** **************************************************************** 

SOURCE SUMS OF SQUARES DEGREES OF FREEDOM MEAN SQUARES 


LACK OF FIT 0.17165756 6 0. 28609593 E- 0 1 

REPLICATION 0.91677475 13 0 . 9 1 67 74 75 E- 01 

RESIDUAL 1.08843231 16 0. 6802 7020 £- 01 


F = MSI LOF) /MSlAfcPS) = 0.312 

********************* ********************************************************** 


SUMS 3F SQUARES DUE TO EACH VARIABLE IF IT WERE LAST TO ENTER REGRESSION 

1 164.5792 

2 20.16795 

3 0. 32 772 5£- 02 

4 10.450 84 

5 0. 6552616-01 

6 0.41 U L6E-02 

7 114.8713 

9 0.27 1357E- 01 


STANDARD DEVIATION OF REGRESSION COEFFICIENTS IOERIVEO FROM 01 AGONAL E. EM ENT S OF IX TRANSPOSE XUNVERSE MATRIX ) 

0 0.90430 IE- 01 

1 0.6 1674 5E- 01 

2 0.5927756-01 

3 0.62869SE-01 

4 0.88908 76-01 

5 0.8548786-01 

6 0.86 L5 1 3E~ 01 

7 0.478488E-01 

9 0.4798236-01 


(X TRANSPOSE X) INVERSE MATRIX 


ROW 

1 

0 • 5252686-01 








ROW 

2 

-0.500332E-02 

0. 4852 34E-01 







ROW 

3 

-0.1 16409E - 01 

0. IU7UE-01 

0.5458ZB E-3 1 






ROW 

4 

0.190643E-01 

-0. 1745416-01 

-0.2842046-01 

3.139 159 





ROW 

5 

0.886623E-02 

-0. 2168116-01 

-0. 1919416 -OL 

0 . 335 285E-0 1 

0. 100920 




ROW 

6 

-0 . 2470746-01 

0. 7L7437E— 02 

0.1935376-01 

-0.313653E-01 

-0.1331546-01 

0. 1 02493 



ROW 

7 

-0.9208226-02 

0. 463255E-02 

0.2752636-02 

-0.462653E-02 

-0. 344396E- 02 

0.49691 96-02 

0. 316L646-0L 


ROW 

8 

-0 .23939 IE -02 

0. 5660526-02 

0. 9909 53 E-02 

-0.6185896-02 

-0.516528E-0? 

0.40764 lc -02 

0. 970715F-0? 

0.3179316 
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SAMPLE RAPIER PROBLEM 


CALCULATED T STATISTICS 

THE T STATISTICS CAN BE USED TO TEST THE NET REGRESSION COEFFICIENTS BID. 

47 .67 303 
16.68845 
0.212736 
12.01326 
0.951245 
0.238269 
39.82824 
0.612147 

UNDER NULL HYPOTHESIS THE INTERVAL 1- T * Tl WHERE T IS GIVEN ABOVE ♦ HAS APPROX PROBABILITY LISTED BELOW. 
MINUS SIGN INDICATES PROB EXCEEDS .999. * 

1 -0.999 

2 -0.999 

3 C.165 

4 -0.999 

5 0.644 

6 0.185 

7 -0.999 

9 0.451 


THE DESIRED 

VALUE OF PROBABILITY IS 95.0 

PERCENT 





THE TERM XI 3) IS 

BEING OELETEO 






CORRELATION 

COEFFIC IENTS 







ROW l 

l .000000 







ROW 2 

-0 

1. 000000 






ROW 3 

-0 .1L9455 

0. 1 1 7851 

1.003333 





ROW 4 

-0. 1924 14E- 01 

0.237289 

-0.167739 

1.033000 




ROW 5 

0.259760 

0 

0.167739 

0. 135 135E-01 

L. 000000 



ROW 6 

0.209642 

-0.671519E-01 

-0 

0. 382427E-02 

-0. 38242 7E- 0 2 

l. 000000 


ROW 7 

-0.6262036-01 

-0.671519E-01 

-0 

0 .382427E-02 

-0.382427E-02 

-0.298701 

l , 003330 

SAMPLE RAPIER PROBLEM 







EACH COLUMN 

CONTAINS THE COEFFICIENTS FOR 

ONE DEPENDENT 

TERM 





CONSTANT TERM 1801 
9.989235 

REGRESSION COEFFICIENTS <Bl,...,BK) 
L 2.937356 
2 0.991988 

4 1.061119 

5 0.766L66E-01 

6 0.25 268 8E- 01 

7 1.906408 

9 0. 3 1 80 04E- 01 


ANOVA OF REGRESSION ON DEPENDENT VARIABLE l 

******************************************************************************* 

SOURCE SUMS OF SQUARES DEGREES OF FREEDOM MEAN SQ J ARES 


REGRESSION 425.384892 7 63.7692695 

RESIDUAL 1.09170532 17 0 . 64 2 1 7 960 E-0 l 


TOTAL 426.476597 24 

***************************************** **************** ********************** 

R SQUARED = SSQIRE3) / SSQITQT) = 0.997440 R = .998719 

STANDARO ERROR OF ESTIMATE 0.253413 

USING POOLING STRATEGY 2 THE ERROR MEAN SQUARE = 0 . 7 24 1 5L6 6-0 l WITH DEGREES OF FREEDOM = 15 

F=MSIREG I/MSIERRl* 839. 18 COMPARE TO F( 7, 15) 


ANOVA OF LACK OF FIT 

*************** **************************************************************** 

SOURCE SUMS OF SQUARES DEGREES OF FREEOGM MEAN SQUARES 


LACK OF FIT 0.17493057 7 0. 24990082 E-01 

REPLICATION 0.91677475 10 0.916774756-01 

RESIOUAL 1.09170532 17 3 . 642 l 7960E-01 


F = MSI LOF) /MSIREPS) = 0.273 

************ *** **************************************************************** 


SUMS OF SQUARES DUE TO EACH VARIABLE IF IT WERE LAST TO ENTER REGRESSION 

1 L72.4090 

2 21.28251 

4 11.93265 

5 0.623348E- 01 

6 0.66 76 68E- 02 

7 115.4596 

9 0.33 71 56E- 01 
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STANDARD DEVIATION OF REGRESSION COEFFICIENTS {DERIVED FROM DIAGONAL E_ EH ENT S OF (X TR ANSPOS r X)INV C RSE MATRIX I 

0 0-90 3708E- 01 

1 0.60 1993E- 0 1 

2 0-5 7 8642E- 01 
A 0. 8266 2 9E- 0 I 

5 0- 82 5796E- 01 

6 0.832182E-01 

7 0.477436E-01 

9 0.466049E- 0 1 


(X TRANSPOSE XJ INVERSE MATRIX 


ROW 

1 

0. 50044 IE-01 




ROW 

2 

-0.262083E-02 

0. 462 3 70E-01 



ROW 

3 

0 • 130030E-01 

-0. 1163 74E-01 

0.9436086-01 


ROW 

4 

0.477268E-02 

-0. 177527E-01 

0.2053 45 E-0 1 

0.941 707E-0 l 

ROW 

5 

— 0.2058 04E - 0 1 

0.321397E-02 

-0.209897E-01 

-0.65I070E-02 

ROW 

6 

-0. 862 l 17E-02 

0.406919E-02 

—0. 319328E— 02 

-0.247600E-02 

ROW 

7 

— 0. 2804756 - 03 

0. 363238E-02 

-0. 10261 IE-02 

-0 . 163055E-02 


SAMPLE RAP IER PROBLEM 


0.9 56329E- 01 
0. 399 332E-0 2 
0.563249E-03 


0. 31477 6E -01 
0. 9207405-02 


D. 2999V OE-DI 


CALCULATED T STATISTICS 

THE T STATISTICS CAN BE USEO TO TEST THE NET REGRESSION COEFFICIENTS Bill. 

48.79386 

17.14339 

12.83671 

0.927792 

0.303645 

39.93009 

0.682340 

UNDER NULL HYPOTHESIS THE INTERVAL I- T * T > WHERE T IS GIVEN ABOVE. HAS APPROX PROBABILITY LISTED BELOW. 
MINUS SIGN INDICATES PROB EXCEEDS .999. 

1 -0.999 

2 -0.999 

4 -0.999 

5 0.632 

6 C.233 

7 -0.999 

9 0.495 


THE DESIRED VALUE OF PROBABILITY IS 95.0 PERCENT 
THE TERM XI 6) IS BEING DELETED 


CORR ELAT 

ION 

COEFFICIENTS 





ROW 

1 

l .000000 





ROW 

2 

-0 

l. 000000 




ROW 

3 

-0.119455 

0. I 1 7851 

1.000333 



ROW 

4 

-0. 1924 14E - 01 

0.237289 

-0.167789 1.033000 



ROW 

5 

0.209642 

-0.671519E-01 

-0 0.382427E-02 

1.000000 


ROW 

6 

-0.626203E-01 

-0.671519E-01 

-0 0 . 38242 7E-02 

-0.298701 

1.000000 

SAMPLE 

RAPIER PROBLEM 





EACH COLUMN 

CONTAINS THE COEFFICIENTS FOR 

ONE DEPENDENT TERM 




CONSTANT TERM CBO) 


9.987362 

REGRESSION COEFFICIENTS (BI.....BK) 

1 2.942794 

2 0.991L39 

4 1.066665 

5 0.78 3369E- 01 

7 1.905353 

9 0.3I6516E-0I 


ANOVA OF REGRESSION ON DEPENDENT VARIABLE l 
*************** *************************************+*****+*+*********+*+****** 

SOURCE SUMS DF SQUARES DEGREES OF FREEDOM MEAN SQUARES 


REGRESSION 425.378220 5 70.8963690 

RESIOUAL 1.09837723 18 0 . 6 1020957 E- 0 1 


TOTAL 426.476597 

*********************************************** 

R SQUARED = SSJ1REG1 / SSQIT0T1 = 0.997425 
STANDARD ERROR OF ESTIMATE 0.247024 
USING P30LING STRATEGY 2 THE ERROR MEAN SQUARE 
F=MS(REG)/HSIERR>= 979.03 COMPARE TO F( 


24 

.******************************* 

R = . 99 B 7 U 

= 0.7 24 1516 E-0 1 WITH DEGREES OF FREEDOM = 15 

6. 151 
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ANOVA OF LACK OF F 1 T 

******************************************************************************* 

SOURCE SUMS OF SOUARES DtGREES OF FREEDOM MEAN SjJARES 


LACK OF FIT 0.18160240 8 0. 2 2 700 3 ID E-0 1 

REPLICATION 0. >1677475 ID 0. 9 1 67 74 75 E- 01 

RESIDUAL 1.09837723 13 0 . 6 LU20957 E-OI 


F = MS(LOF) /MSUEPS) = 0.248 

********************* ********************************************************** 


SUMS OF SJUARES DUE TO EACH VARIABLE IF IT WERE LAST TO ENTER REGRESSION 
l 189.8499 
i 21.29583 

4 12.67660 

5 0.654736E-01 

7 II5.946C 

9 0.334045E-01 


STANDARD DEVIATION OF REGRESSION COEFFICIENTS (DERIVED FROM DIAGONAL ELEMENTS OF (X TRANSPOSE XMNVERSE 

0 0. 901600E- 01 

1 0.5747 3 7E- 0 1 

2 0.57796 5E - 01 

4 0.8061 9 7E- 01 

5 0. 82 38 5 OE- 01 

7 0.476 1 7 OE- 01 

9 0. 466024E- 01 


(X TRAN 

SPOSE 

X) INVERSE MATRI 

X 



ROW 

1 

0 . 4561 52E — 01 




ROW 

2 

- 0 . 19291 8E - 02 

0.461290E-01 



ROW 

3 

0.848598E-02 - 

0. 109320E-01 

Q.d97539E-0I 


ROW 

4 

0.33715 6E- 02 - 

0. 175339E-01 

0.191055E-01 

3 .937 275E-0 l 

ROW 

5 

-0. 776179E-02 

0. 3 93498E-02 

-0.231532E-02 

-0.22D4L3E-Q2 

ROW 

6 

-0.159262E-03 

0. 36 1 345E-02 

-0.902^ 3 7 E-3 3 

-3 . 16422IE-02 


SAMPLE RAP IER PROBLEM 


0.3131 09E - 0 L 

0.918388E-02 0.299907E-01 


CALCULATED T STATISTICS 

THE T STATISTICS CAN BE USED TO TEST THE NET REGRESSION COEFFICIENTS Bill. 

51.20241 
17. 14875 
13.23082 
C. 950 864 
40.01411 
0.679185 

UNDER NULL HYPOTHESIS THE INTERVAL 1-T.T) WHERE T IS GIVEN ABOVE. HAS APPROX PROBABILITY LISTED BELOW. 
MINUS SIGN INDICATES PROB EXCEEDS .999. 

1 -0.999 

2 -0.999 

4 -0.999 

5 C.643 

7 -0.999 

9 C.493 


THE DESIRED VALUE OF 
THE TERM X 


PROBA0I LI TY 
9) IS BEING 


S 95.0 PERCENT 
DELETED 


CORRELATION 

COEFFICIENTS 



ROW 

1 

1 .000000 



ROW 

2 

-0 


l. 000000 

ROW 

3 

-0.119455 


0. 11 7051 

ROW 

4 

-0. 192414c- 

01 

0.237289 

ROW 

5 

0 .209 642 


-0.671519E 

sample 

RAPIER PROBLEM 



EACH COLUMN 

CONTAINS THE 

COEFFICIENTS I 

CONS TAN 

T TERM (BO) 



10 

.02609 



REGRESS 

I3N 

COEFFICIENTS 

( B l , 

,. . . . B K) 


1 . 0000 DO 
-0.167739 


l. 003000 
0. 38242 7E-02 


2.942962 
0.907325 
1.067618 
0.80070 1E-0L 
1.895660 


1.000000 


/ 


MATRIX I 
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ANOVA OF REGRESSION ON DEPENDENT VARIABLE l 

**************************** ************************** ************************* 

SOURCE SUMS OF SQUARES DEGREES OF FREEDOM MEAN S0JARES 


REGRESSION 425.344822 5 85.0689640 

RESIDUAL 1.13177490 19 0. 595671OD6-01 


TOTAL 426.476597 2* 

************************************************************************ ******* 

R S3UARED * SSUIREGl / SSQITOT) = 0.997345 R = ,998672 

STANDARO ERROR OF ESTIMATE 0.244054 

USING POOLING STRATEGY 2 THE ERROR MEAN SUUARE = 3 .72% 151 6E-Q l WITH DEGREES OF FREEDOM = 15 

F-MS(REG)/MS(ERR)= 174. 74 COMPARE TO F I 5, 151 


ANOVA OF LACK OF F I T 

******************************* ****** ****** ************************************ 

SOURCE SUMS OF SQUARES DEGREES OF FREEDOM MEAN SQUARES 


LACK OF FIT 0.21500015 9 0. 230809O6E-O1 

REPLICATION 0.9L677475 10 0 . 9 167 74 75 E-01 

RESIDUAL 1.13177490 19 0. 59567 lOOE-Ol 


F * MSI L0F1/MSIREPS1 = 0.261 

************ *** ************ **************************************************** 


SUMS OF SQUARES DUE TO EACH VARIABLE IF IT WERE LAST TO ENTER REGRESSION 

1 189.8751 

2 21.33362 

4 12.70310 

5 0.684685E-01 

7 126.0952 


STANDARD DEVIATION OF REGRESSION COEFFICIENTS IDERIVEO FROM DIAGONAL ELEMENTS OF (X TRANSPOSE XJINVERSE MATRIX ) 

0 0.688638E-01 

1 0.574732E-01 

2 0.57 52 3 2E -01 

4 0.806076E-01 

5 0.82345 5E- 01 

7 0.45428 2E- 0 1 


IX TRANSPOSE X) INVERSE MATRIX 


ROW 

1 

0 . 456 143E-01 




ROW 

2 

-Q. 1909996-02 

0. 456936E-0 1 



ROW 

3 

0 • 848 1 1 9c -02 

-0. 1 082 33 E-Ol 

0.897258E-01 


ROW 

4 

0. 3362846-02 

-0. 1733606-01 

0. 190551E-01 

0,935 375E-0 l 

ROW 

5 

—0.771 3Q2E-02 

0.282846E— 02 

-0. 204045 E-02 

-0.170 125E-02 


SAMPLE RAPIER PROBLEM 


0.284985E-01 


CALCULATED T STATISTICS 

THE I STATISTICS CAN BE USED TO TEST THE NET REGRESSION COEFFICIENTS Bill. 

51.20581 

17.16396 

13.24464 

0.972368 

41.72868 

UNDER NULL HYPOTHESIS THE INTERVAL I-T.T) WHERE T IS GWEN ABOVE. HAS APPROX PROBABILITY LISTE3 BF LOW. 
MINUS SIGN INOICATES PROB EXCEEDS .999. 

1 -0.999 

2 -0.999 

4 -C.999 

5 0.6 54 

7 -0.999 


THE DESIRED VALUE OF PROBABILITY IS 95.0 PERCENT 
THE TERM XI 5) IS BEING DELETED 


CORRELATION 
ROW 1 

ROW 2 

ROW 3 

ROW 4 


OEFF IC IfcNTS 
t. 000000 
-0 

-0.119455 

0.209642 


I. 000000 
0.117851 
-0.671 519E 


1.000300 
-01 -0 


1.000000 


SAMPLE RAPIER PROBLEM 

EACH COLUMN CONTAINS THE COEFFICIENTS FOR ONE DEPENDENT TERM 
CONSTANT TERM IBO) 

10.03236 

REGRESSION COEFFICIENTS !B1*...,BK) 

1 2.940086 

2 1.002149 

4 1.051323 

7 1.897115 
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anova of regression on dependent variable l 

**************************** ************************************ *******«,******„ 

SOURCE SUMS OF SQUARES DEGREES OF FREEDOM MEAN SQUARES 


REGRESSION 425.276363 4 136.319091 

RESIDUAL 1.20023346 23 0. 6001 1 6 72 E- 01 


TOTAL 426.476597 24 

************ ********* ********************************************************** 

R SJUARED = SSJ(REG) / SSQ(TGT) = 0.997US R = .998592 

STANDARD ERROR OF ESTIMATE 0.244973 

USING POOLING STRATEGY 2 THE ERROR MEAN SOUARE = 3 . 7 24 l 5 l 6 E-0 1 WITH DEGREES OF FREEDOM = 15 

F=MS(REG )/M SI ERR I = 468. 19 COMPARE TO Ft A, 151 


ANOVA OF LACK OF FIT 

************ ♦ ** *** * **************** *** ********** *** ****** 4 > * + ** *************** 

SOURCE SUMS OF SOUARE S DEGREES OF FREEDOM MEAN SQ J ARES 


LACK OF FIT 

0.20345871 

13 

3. 28345871E-01 

REPL icat iqn 

0.916774 75 

10 

0.9 16774 75 E-01 

RESIDUAL 

1.20023346 

23 

0. 6001 L672E-01 


F = MSI LOFI/MSIREPS) = 0.309 

*************** ****************************************** ********************** 


SUMS OF SQUARES DUE TO EACH VARIABLE IF IT WERE LAST TO ENTER REGRESSION 

1 190.0073 

2 23.63953 

4 12.87474 

7 126.4260 


STANDARD DEVIATION OF REGRESSION COEFFICIENTS (DERIVED FROM 01 AGONAL ELEMENTS OF 

0 0.66633 7E- 0 l 

1 0. 5 7 39 7 IE- 0 1 

2 0.55466 IE- 01 

4 0. 788463E-01 

7 0. 4540 3 6E- 01 


(X TRANSPOSE 


XIINVFRSF MATRI X I 


( X TRAN 

SPOSE 

X) INVERSE MATRI 

X 



ROW 

1 

0.454936E-01 




ROW 

2 

-0*12874 0£- 02 

0. 42 4 84 IE— 0 1 



ROW 

3 

0. 77968 36-02 - 

0. 729524E-02 

0.8584876-01 


ROW 

4 

-0.76519 2E- 02 

0. 251349 E-02 

-0. 1694246-02 

0.2846766-01 

SAMPLE 

RAP IER PROBLEM 





CALCULATED T STATISTICS 

THE T STATISTICS CAN BE USED TO TEST THE NET REGRESSION COEFFICIENTS B( I ) . 

51.22363 
18 .06778 
13.33382 
41 .78 336 

UNDER NULL HYPOTHESIS THE INTERVAL (-T.T) WHERE T IS GIVEN ABOVE, HAS APPROX PROBABILITY LISTE3 BELOW. 
MINUS SIGN INDICATES PROB EXCEEDS .999. 

1 - C. 999 

2 -0.999 

4 -0.999 

7 -0.999 


THE DESIRED VALUE OF PROBABILITY IS 95.0 PERCENT 
SAMPLE RAPIER PROBLEM 

FOR EACH OEPENDENT TERM AND OBSERVATION IS PRINTED 
OBSERVED RESPONSE (Y OBSERVED! 

CALCULATED RESPONSE IY CALCI 
RESIDUAL {YOBS- YCALC =YDIF ! 

STANDARDIZED RESIDUAL l l ) 


Y OBSERVED 9.1700 

Y CALC 9.0386 

Y DIF 0.1314 

S TUOENT I ZED 0.4884 


Y OBSERVED 

Y CALC 

Y OIF 

STUDENT IZED 


Y OBSERVEO 

Y CALC 

Y OIF 

STUDENT IZED 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZED 


Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZED 


12.760 
12.816 
-0 .5608E-01 
-0 .2084 

12.970 
12.616 
0.1539 
0 . 5720 

9. 1100 
8.9402 
0. 1698 
0.6310 


8.9600 
8.9402 
0. 1979E-01 
0. 7354E-01 


43 


Y OBSERVED 

Y CALC 

Y OIF 

STUDENTIZED 

I 7.030 
16.92 3 
0. 10 70 
0.3975 

Y OBSERVEO 

Y CALC 

Y OIF 

STUDENTIZED 

9.0500 
9.0386 
0. U44E-01 
0.4252E-01 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZED 

8.8600 

9.0386 

-0.1786 

-0.6635 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENTIZED 

12.600 

12.816 

-0.2161 

-0.8030 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZEO 

13.210 
12.816 
0.3939 
1 . 4638 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENTIZED 

17.200 
16.923 
0.2770 
l .0292 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT I ZED 

l 7.040 
L6.92 3 
0.1170 
0.4347 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZED 

9.6100 
10.032 
-0.4224 
-1 . 5695 

Y OBSERVED 

Y CALC 

Y OIF 

STUDENT IZED 

10.010 
10.032 
-0.2236E-01 
-0 • 8308E-0 1 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZED 

10.120 
10.032 
0 • 8 764E-01 
0.3257 

Y OBSERVED 

Y CALC 

Y OIF 

STUDENT IZED 

9.9500 
10.032 
-0. 8236E-01 
-0. 3060 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZED 

11.780 

11.741 

0.393&E-01 

0.1463 

Y OBSERVED 

Y CALC 

Y OIF 

STUDENT IZED 

23.830 

23.501 

0.3290 

1.2226 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZED 

22.900 

23.501 

-0.6010 

-2.2333 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENTIZED 

7.9900 
8. 0281 
-0.3806E-01 
-0.1414 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZED 

1 2. 110 
12.037 
0.7335E-01 
0.2726 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZED 

1 1.700 
12.037 
-0.3367 
-1 .2510 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENT IZED 

10.110 
10.032 
0. 7764E-01 
0.2885 

Y OBSERVED 

Y CALC 

Y DIF 

STUDENTIZED 

10.010 
10.032 
-0 . 2236E-0 l 
-0.8308E-01 

Y OBSERVED 

Y CALC 

Y DIF 

STUOENT IZ ED 

10.020 
10.032 
-0. 1236E-01 
-0.4592E-01 



SAMPLE RAPIER PROBLEM 

SKEWNESS (SHOULD BE NEAR ZEROl 
0. 1928 

KURTQSIS (SHOULD BE NEAR THREE) 

1. 69*9 

CHI-SQUARE IS NOT COMPUTED FOR LESS THAN 30 DEGREES OF FREEDOM FOR ERROR. 


FOR EACH SET OF INDEP VARIABLES THERE IS COMPUTED... 

PREDICTED RESPONSE 

VARIANCE OF REGRESSION LINE 

STANDARD DEVIATION OF REGRESSION 

VARIANCE OF PREDICTED VALUE 

STANOARD DEVIATION OF PREDICTED VALUE 

INPUT DATA FOR THIS PREDICTED RESPONSE 
-l. 000000 I. 000000 I. 000000 -0 -3 

INDEPENDENT TERMS ACCORDING TO FINAL REGRESSION MODEL 
- 1. 000000 1.000000 -1.000000 l .00033 3 

PREDICTED RESPONSE FOR ABOVE INDEP VARIABLES 
8.9*0211 
0. 1928^76-01 
0. 138869 
0.9 16998E— C 1 
0.302820 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, October 16, 1969, 

129-04. 
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APPENDIX A 


PROGRAM DOCUMENTATION AND LISTINGS 

The contents of this appendix include a flow chart of the program, a listing of the 
routines used in RAPIER and their major functions, the call structure of the program, 
a dictionary of the program, and the listing. 


General Mathematical and Logical Flow of Program 


The flow of operation in RAPIER is illustrated in figure 6. 



Figured - Flow chart for RAPIER. 
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Routines and Their Major Functions 


FORTRAN 

name 


Function of routine 


BORD 

CHISQ 

HIST 

INVXTX 

LOC 


MAT IN V 

MFIX 

PREDCT 

RAPIER 


RECT 

RSTATS 

SUMUPS 

TRAN 

TRIANG 

TTEST 

EIGEN 


Inverts symmetric matrix A of order n by adding bordering 
column to already inverted matrix of order n - 1 

Computes residuals at observed points and chi-square statistic to 
test goodness of fit 

Prints histogram of residuals 

Inverts symmetric matrix by Gauss elimination 

When given row and column indices of symmetric matrix element, 
it computes location this element would have if only upper tri- 
angular part were stored as vector. 

Controls inversion process; computes regression coefficients; 
computes eigenvalues and eigenvectors of X'X if requested 

Prints and truncates X’X and computes C 

Computes predicted values, variances, and standard deviations of 
regression line and further observations at specified points 

Executes overall problem control; computes replication error; 
controls deletion of variables when given results of t-tests; 
controls most input and output 

Writes rectangular matrix 

Computes regression statistics and writes regression and lack-of- 
fit analysis of variance tables 

Constructs X'X and X’y matrices one observation at a time, in 
double precision 

Performs transformations 

Writes lower triangular part of symmetric matrix 

Computes t-statistics and their significance levels; determines 
which variable should be deleted 

Computes eigenvalues and eigenvectors of input symmetric matrix 


111 
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Call Structure of Program 

The call structure of the program is illustrated in figure 7. 



Dictionary of Program 


FORTRAN 

name 

Mathematical 

symbol 

Description 

ACTDEV 

e i 

Error in observation i; difference between ob- 
served and predicted response 

ZEAN 

B 

BO 

E(X), m 

Expected or average value (of X) 

Regression coefficients other than the constant 
Constant regression coefficient 

CHI 

CON 

x 2 

Chi- squared statistic 

Constants used in transformations, and results of 
transformations 

CORR 

c 

Correlation matrix 
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FORTRAN 

name 


Mathematical 

symbol 


Description 


DELETE 

ERRMS a 2 

EIG 

FMT 

FMTTRI 

IDENT 

IDOUT 

IFCHI 

IFSSR 

IFTT 

IFWT 

INPUT 
INPUT 5 
INTER 

IOUT 

JCOL 

KONNO 

LENGTH 

LIST 


Lexical variable set to TRUE when deletion of 

terms is desired 
o 

Estimate of a used in hypothesis tests 

Overflow area for storing part of modal matrix 

Variable input format 

Format for printing matrix 

First identification printed at top of each page 

Original sequence number of each term relating 
reduced models to original model 

Logical variable set to TRUE if chi-square option 
is desired 

Logical variable set to TRUE if sequential regres- 
sions are desired 

Logical variable set to TRUE if t- statistics are 
desired 

Logical variable set to TRUE if all weights of ob- 
servations are 1. 0 

Input logical tape unit number for data 

Set equal to 5 to denote input device is card reader 

Tape unit where input data is stored for chi-square 
or prediction routines 

Sequence number of term among those remaining 
which is to be deleted 

Total number of independent and dependent terms 
in regression model 

Number of constants originally supplied for trans- 
formations 

Number of locations in correlation matrix storage 
area currently needed 

Set equal to 6 to denote output device is printer 
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FORTRAN 

name 


Mathematical 

symbol 


Description 


NARAY 


Number of replications per replicate set 

NCON 


Array containing addresses in CON array for use 
in transformations 

NERROR 


Degrees of freedom for error mean square estimate 

NODEP 


Number of dependent variables 

NOOB 

N 

Number of observations 

BZERO 


Logical variable set to TRUE if constant b Q coeffi- 
cient should be in regression model 

NOTERM 

J 

Number of terms in current regression model 

NOVAR 

K 

Number of independent variables to be read 

NPDEG 

NPDEG 

Pooled degrees of freedom for replication error 

NRES 

N - J - D 

Degrees of freedom for estimation of residual vari- 
ance 

NTERM 


Array containing locations of terms in CON 
array that should be in regression model 

NT RAN 


Array containing transformation codes for use in 
performing transformations 

NTRANS 


Number of transformations to perform 

NWHERE 


Location in X array of first dependent variable; 
used in prediction routine to adjust for deleted 
terms 

NXCOD 


Array containing addresses of variables (or terms 
with address > 60 ) for use in transformations 

P 


Probability that the interval (-t, t) must have before 
a term is considered to be significant 

POOLED 

SSQ(REP) 

Array containing pooled sums of squares from rep- 
lications for each dependent term 

PREDCT 


Logical variable set to TRUE if prediction option 
is desired 

RELSKW 

RELSKW 

Skewness of distribution of residuals 
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FORTRAN 

name 

Mathematical 

symbol 

Description 

RELKUR 

RELKUR 

Kurtosis of distribution of residuals 

REPS 


Logical variable set to TRUE if there are replicate 
sets in the data 

REPVAR 


Array containing replication variance of each de- 
pendent term 

RE SMS 


Array containing residual mean square or variance 
of each dependent term 

RWT 


Reciprocal of total weight 

SUMX 

X)x,Hy 

Array containing sums of independent and dependent 
terms 

SUMX2 

Sx 2 ,Zy 2 

Array containing sum of squared independent and 
dependent terms 

SUMXX 

X'X 

Sums of squares and crossproducts matrix, and 
variance- covariance matrix of independent terms 

SUMXXI 

(X'X)" 1 

Variance-covariance matrix of estimated regres- 
sion coefficients 

SUMXY 

X'y 

Array containing sums of crossproducts of inde- 
pendent terms with dependent terms 

TOTWT 

Sw i 

Sum of weight of observations 

X 


Before transformations are performed, this con- 
tains the variables as read in. After the trans- 
formations are performed, appropriate data from 
the CON array are placed here according to in- 
formation on the terms cards. 

NLOF 

N - J - NPDEG - D 

Degrees of freedom for estimating variance due to 
lack of fit 

NREG 

J 

Degrees of freedom for determining variance due 
to regression 

NTOT 

N-D 

Total degrees of freedom 

RNLOF 


Reciprocal of degrees of freedom for lack of fit 

RNREG 


Reciprocal of degrees of freedom for regression 
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FORTRAN 

name 

RNRES 

STORYI 

STORYC 

ST ORYX 

SATRTD 

ECONMY 

XCHK 


Mathematical 

symbol 


Description 


Reciprocal of degrees of freedom for residual 

Logical variable set to TRUE if product C • C" 1 
is to be computed and printed 

Logical variable set to TRUE if eigenvectors and 
eigenvalues of C are to be computed and printed 

Logical variable set TRUE if eigenvectors and 
eigenvalues of X’X are to be computed and 
printed 

Logical variable indicating that there are no de- 
grees of freedom for residual if TRUE 

Logical variable indicating suppress printout of 
x'x, x'x deviations, and C if TRUE 

Array used in checking if all values of independent 
terms are the same within a replicate set 


Program Listing 


* I BFTC BLCV 


6LCCK DATA 


COMM CN/ SHALL / 

BYPASS * BZERC 

1, DELETE, 

FIRST, IFCHI 

X 

I FTT , 

I FRT , 

INPUT , 

INPUTS, 

X 

ISTRAT, 

JCOL, 

KONNO, 

length. 

X 

N ERRCR t 

NOBEP, 

NOOB, 


X 

NOVAP , 

NPCEG , 

NRES , 

NTRANS, 

X 

P, 

PREDCTi 

REPS, 

RHT, 

X 

STORYI | 

STORYC, 

STORYX, 

TOTHTV 

X 

ERRFXC, 

ECONMY 




LOGICAL ECCNKY 

OCUetE PRECISION RWT/TOTWT, WEIGHT 
COMMON /FRMTS/ FMTIHI.FHTTRIIIAI 
DATA INTER/3/, INPUH/5/. LIST/6/ 

DATA IFMTTRII II. I = i.4)/6MI5H RO, 6HW 1 5,2, 6HX,t8Gl, 
ENC 


IFSSR, 

INTER, 

LIST, 

NOTERM, 

NWHERE, 

HEIGHT, 


6H5 . 6 ) ) / 
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SIBFTC RAPIER 


C 

C THIS IS RAP 1 ER»M A IN PROGRAM FOR RE GRESSICN ANALYSIS PROVIDING 

C INTERNAL EVALUATION 0F RESULTS. 

C*** ***************************************************** ********* ****** 

C 

COMMON/ BIG/ SUMXX < 1830), SUMXXII 1830) , E I G ( 1830 ) , SUMXY ( 60 ,9 ) * 

X B(6C,9),C0RR(1830) 

CGMMGN/MiED/ B019), C0NI99), ERRMSI9), 

X I CENT I13)| IDOUT ( 60 ) , NCON ( 200 ) i 

X NTERM < 60 ) » NTRAN(100>| NXCOD(IOO), P0CLED<9), 


XREPVAR( 9U 
XSUMX2 (69 I ». 

CCMMCN /FRMiTS/ 


X< *9), 

FMT { 13) fFMTTRI ( 14 ) 


RESM5I9):, 
ZE AN 1 69 I t 


SUMX ( 1 38 ) • 
SUMY2I 181 


LOGICAL ECONNIY 

DOUBLE PRECISION RWT^TOTWT, WEIGHT 


LOGICAL BYPASS , 

XIFSSR, f FTT , 

XSTCRYC, STORYX, 

LOGICAL XSAVE 
DIMENSION XCHKI60) 
COMMCN/CNTRS/ I, 

X I NEW , INOCH, 

X I RO, I RiP * 

X K, KBAR 


BZEROf 
IFWT* 
STORY I, 


DELETE, 
REPS* 
FIRST 


IFCHI, 
PREDCT, 
,ERRFXD 


IBC, 

I OLD , 

IS! 


IC, 

I OUT, 
ITC, 


I COL , 
IR* 

J , 


C******** ********** ***************************************************** 

c 

EQUIVALENCE I NAR AY ,C0RR ) | (S*BO)| ( SSQ ,REPVAR) , |ZCUT,SUMXX) 
DIMENSION N A-RAYt 18301,S<9) ,S$Q(9i 
DIMENSION Z0UTI1) 

C*** ***** ** ** ********* ************************************************** 
C ZERO OUT ALL DATA ARRAYS EACH NEW DATA SET 

100 «ZCUT=;9470 

DO 1C I J* 1 , I ZOUT 

101 SUMX X ( J ) =* OiO 
C 

C ******* **4 *********** ************************************************** 


c 

c 


READ IDENTIFICATION CARD AND OPTIONS CARD 

REACT INPUTS, 110) l, I DENT 
WRIT£!LIST,111) IDENT 
FIRST=.TRUE. 

ERRFXC=. FALSE. 

113 I F ( I I 120,120,115 
115 REACT INPUTS, 300) FMT 
WR I TE t L I ST , 30 1 ) FMT 
1=1-1 
GO TC 113 

120 REACT INPUTS, 1282 ) NOVAR , NODE P , NOTERM', NOOB 
WR I TE < L I ST , 1 2831 NOVAR , NODEP ,NOTERM , NOOB 

RE AC I IN PUT 5, 117) BZERffl ,.I F TT , I FWT , IFCH I , STORYC , STORYX , STORY I , 
X rECONMYfr ISTRAT 

WRITE (LI ST, 118) B ZERO *1 FTT, 1 FWT , IFCHI , STCRYC , STORYX , STORY I , 
X rECONMY, ISTRAT 


IFSSR 


IFSSR 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


GOPMON/ SMALL/ 

BYPASS < B2EftO» DELETE* 

FIRST:, IFCHI, 

IFSSR, 

14 

X 

I FTT , 

IFWT, 

INPUT, 

INPUTS* 

INTER, 

15 

X 

I STRATI 

JC0L , 

KONNO, 

LENGTH* 

LIST, 

16 

X 

NERRORi 

NODEP , 

NOOB, 


NOTERM* 

17 

X 

NOVAR, 

NPEEG, 

NRES , 

NTRANS, 

NWHERE , 

18 

X 

Pt 

PRiDCT | 

REPS, 

RWT , 


19 

X 

STORY I » 

STORYC, 

STORYX, 

TOTWTI 

WEIGHT, 

20 

X 

ERRFXC ♦ 

ECONMY 




21 


22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 
63 


53 


L 



non noon 


in min ii 


i hi ii i i 


(■I ■ 


£**** ******************************************************************* 64 

C THESE ARE INITIALIZATIONS MADE BEFORE EACH SET OF DATA 65 

C ICOL DETERMINES THE NUMBER OF VARIABLES REAO PER OBSERVATION 66 

C JCOL IS THE NUMBER OF TERMS IN THE TOTAL REGRESSION EQUATION 67 

C LENGTH IS THE NUMBER OF STORES NEEDED FOR THE MATRICES 68 

LENGTF= NOTE-RMMNOTERM+ll/2 69 

ICCL=NOVAR + NODEP 70 

JCOL * NQTERN! +NODEP 71 

NWHERE= NOTERM 72 

REWIND INTER 73 

DC 14C J*l,60 74 

IDCUT(J) = J 75 

140 NT ERR I J I* J 76 

DO 145 J*l,100 77 

NXCCD ( J ) * J 78 

NTRAM J )*0 79 

145 NCCNI2*JI=J 80 

C 81 

£**«**%* ******%«****** ** * * *** ******************* 4 **************** ******* 82 

IF ( BZEROI WRITEILISTI190I 83 

1F( .NOT.BZERC) W.R I TE (L I S T , 170 ) 84 

£********** ************************** ***** *************** *************** 85 

READ I INPUT5*282) NTRANS,KONNO 86 

I F ( NTRANS. EQ .0 ) GO TO 255 87 

220 READ ( INPLT5,230) (NTERM(K) ,K=1, JCOL) 88 

WRIT£(LIST,235) I NTE RM ( K ) , K = 1 , JCOL ) 89 

READ ( INPUT5.230I (NXCODI I ) ,NTRAN( I ) ,NCON ( 2* 1-1 ) ,NCON ( 2* I ) , I=1,NTR 90 

AANS I 91 

WRITE (LI ST, 240) <NXC0D< II.NTRANd) , NCON ( 2* I- 1 ) , NCON ( 2* I) , I = l , 92 

X NTRANS) 93 

IF(KONNO) 255,255,256 94 

250 REAO < INPUTS, 260) <CBN< I ), 1=1, KONNO) 95 

WR ITE (LIST, 262) ( CON t I ) ,1 = 1 , KONNO) 96 

******************************************************************* 97 

C 98 

255 REACI INPUTS, 257) DELETE, P 99 

1 F ( CELETE ) IFTT= .TRUE. 100 

C 101 

C*** ***************************************************** *************** 102 

IF THERE ARE REPLICATED POINTS READ IN THE NUMBER OF POINTS AND 103 

THE NUMBER OF REPLICATIONS. SINGLE OATA POINTS ARE DATA POINTS 104 

REPLICATED ONCE. OBSERVED OATA MUST BE ARRANGED IN THE ORDER 105 

IMPLIED HERE. 106 

READI INPUT5.257) REPS 107 

XSAV£=. FALSE! 108 

265 IF< .NOT. REPS) GO TO 290 109 

RE AC I INPUT5,282) I REP* (NARAY ( I ) , 1=1, IREP ) 110 

NPDEGnO HI 

I R E P= 1 112 

IC = NARAY ( 1 ) 113 

XS AVEt. TRUE. 11* 

DO 315 1*1, NODEP 115 

POCLED( I 1 = 0.0 116 

Sm*0.0 117 

315 SSC(l)=OLO 118 

119 

*********************************************************************** 120 

READ VARIABLE FORMAT FOR DATA 121 

290 READUNPUT5, 110) INPUT, FMT 122 

WR ITE (L I ST , 1 1 1 ) FMT 123 

310 TOTWT=O.ODC 124 

WE IGtrT= 110DC 125 

WRITE(LIST,30l) IDENT 126 

C 127 
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C* ************************ ******** ****************************** ******** 128 

C RE AC IN INPUT VARIABLES 129 

00 490 J* 1 , NQOB 130 

330 IFUNCT.IFWT) GO TO 350 131 

340 REAC (INPUT, FMT) ( X ( 1 ) , I *1 , ICOL ) 132 

GO TC 360 133 

350 RE AC ( IN PUT, FMT) < X C I 1*1 = 1, ICOL ) , WEIGHT 134 

360 CONTINUE 135 

‘IF! ECONMY ) WR ITE < L I ST , 38 1 ) J * ( X ( H , I - 1 , ICOL ) 136 

381 FORMAT ( 1H I 4 | 9G1 4* 6/ 1 5X , 9G14*6 ) ) 137 

IF(NTRANS.ECiO) GO TO 450 138 

I F ( ECCNMY I GO TO 390 139 

WRITE(LIST,370)WEIGH¥,.J 140 

WRITE ( LIST,3801(X( 1M = 1, ICOL) 141 

390 CALL TRANS 142 

420 CO 430 K=lf JCOL 143 

1 = NTERM (K 1 144 

X < K ) n CCN ( I) 145 

430 CONTINUE 146 

450 CONTINUE 147 

I F ( ECONMY X GO TO 4609 148 

WRITEfLIST, 460) J 149 

461 WRITE (LIST, 380) <X< 11*1-1, JCOL) 150 

4609 CONTINUE 151 

IF(IFCHI) WRITE( INTER) CXI I ), 1 = 1469), WEIGHT 152 

4F( .NCT.XSAVE) GO TO 4611 153 

DO 4ft 10 K=l, NOTERM 154 

4610 XOMK)=X(K) 155 

4611 CONTINUE 156 

C 157 

£ *********************************************************** ************ 15 Q 

C COMPUTE THE ERROR VARIANCE FROM REPLICATED DATA 159 

IF( .NOT. REPS) GO TO 480 160 

IGCTC =1 161 

IF(NARAYIIREP).GT.l) IG0T0=2 162 

TF(J*GE-*C) NRI T E ( 6 , 462 \ IREP 163 

DO 475 I * 1 , NGDEP 164 

4 F ( I- 1 ) 46 29) 462 9, 464 165 

4629 00 463 K*l, NOTERM 166 

€F(X(K) .NE.XCHKt K) ) €0 TO 2001 167 

463 CONTINUE 168 

464 CONTINUE 169 

K8AR*N0TERM+I 170 

5< I)*S( II+XIKBAR) 171 

SSC( I ) = SSQ( I >+:X< KB AR 1**2 172 

IF(J-IC) 475)465,465 173 

465 GO TC ( 468,466), IGOTffi 174 

466 ZEAN4 l)=S< I ) / FLO AT < N ARA Y t IREP) ) 175 

S5C(D = SSQII) - ZEAN( 1 )*S( I ) 176 

POOLED! I ) ^POOLED! I ) +SSQ ( 1 ) 177 

WRITE (6,467) I,SSQ( II, S( I) , ZEANI 1) 178 

468 €F( IiLT.NOCEP ) GO TO 469 179 

NPCEG=NPDEG+NARAY( IREP) *1 180 

4 REP* IREP+ 1 181 

(C = IC * N AR AY( IREP 1 182 

WR IT E ( L 1ST, 4671) 183 

469 S< I )*0.0 184 

SSCCIl-OkO 185 

XSAVE=.TRUE. 186 

475 CONTINUE 187 

C 188 
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noon non 


C**** ******** ******************************* ****** ******* *************** 189 

C CALCULATE SUMS, SUMS OF SQUARES AND SUMS OF CROSS PRODUCTS. 190 

480 CALL SUMUP 191 

490 CONTINUE 192 

C 490 CONTINUE IS THE END OF THE LOOP FOR READING DATA CARDS 193 

IFI.NCT.REPSI GO TO 496 194 

00 493 1*1, NODEP 195 

■REPVARI 11=5 POOLED! D/FLO AT(NPDEG) 196 

493 CONTINUE 197 

496 CONTINUE 198 

C 199 

C ********* * ********************************************* * *************** 200 

ALL DATA HAS BEEN READ IN AND THE XTRANSPOSEX AND XTRANSPOSEY 201 

MATRIX HAVE BEEN CALCULATED. 202 

NOW WRITE THE MATRICES 203 

CALL MF IX 204 

REWIND INTER 205 

GO TC 640 206 

207 

*********************************************************************** 208 

THIS CODING OELETES THE DATA FROM THE CORRELATION MATRIX 209 

CORRESPONDING TO THE INDEPENDENT TERM DELETED 210 

6500 CONTINUE 211 

4R= ICUT— 1 212 

IC= NCTERM - IOUT 213 

IF ( IC.EQ. 0) GO TO 6700 214 

INCCK= IOUTMR/2 215 

INEW * INOCH 216 

IOLC n INEW * IOUT 217 

IRC=C 218 

IBC=0 219 

ITC=C 220 

DO 6600 I^IOLD, LENGTH 221 

INEW * INEW+1 222 

IGLD* IOLD ♦ ! 223 

IF< ITC.GT.C) GO TO 6540 224 

•1RC= IRC * 1 225 

IF( IRC.GT. IRI GO TO 6530 226 

C0RR1 INEW ) = CDRfil IOED) 227 

GO TO 66C0 228 

6530 IBCMBC • 1 229 

ITC * IBC 230 

ICLD = ICLD + ,1 231 

I RC= 0 232 

6540 ITC * ITC -1 233 

CORR I INEW ) = CORRIIOfcD) 234 

6600 CONTINUE 235 

6700 LENGTH = LENGTH-NOTERM 236 

NCTERM= NOTERM -1 237 

JCOL* NOTERM+NODEP 238 

4FI8ZEROI WRITE ILIST46701 239 

4F( .NOT .BZERO) WR ITE I L I ST, 560 ) 240 

IF < ECCNMY ) GO TO 640 241 

CALL TR I ANG ( CORR ,N0TfRM»8, FMTTR I J 242 

C 243 

C*************** ************************************* ******************* 244 

C INVERT THE CORR COEF MATRIX AND COMPUTE REGRESSION COEFS 245 

C AND SUMS OF SQUARES HUE TO REGRESSION IN THE MATRIX INVERSION 246 

C ROUTINE 247 

640 CONTINUE 248 

CALL MAT IN V 249 

FIRST*. FALSE! 250 

C 251 


56 



C**+* **************************************************** *************** 252 

C WRITE(XTX) INVERSE. THIS MATRIX TIMES ERROR MEAN SQUARE (ERRMS) 253 

C I S TEE VARIANCE-COVARIANCE MATRIX OF REGRESSION COEFFICIENTS. 254 

TFIECONMT ) GO TO 970 255 

WRITEILIST »70O) 256 

CALL TRIANG(SUMXXI,N0TERM,8,FMTTRI ) 257 

c 258 

C *** ************************************************************** ****** 259 

C IF A VARIABLE HAS BEEN DELETED ADJUST COUNTERS AND RECOMPUTE THE 260 

C REGRESSION. IF NO VARIABLE HAS BEEN DELETED CONTROL VIILL PASS 261 

C FROM THE TTEST ROUTINE TO THE CHI-SQUARE OPTION. 262 

970 CONTINUE 263 

IF I .NOT .IFTT I GO TO 1020 264 

980 WRITE (LIST, 301) IDENT 265 

CALL TTEST ($1020) 266 

'IF(NCCEP-l) 985*990,985 267 

985 WRITE(LIST,986) NODEP 268 

N0CEP=1 269 

990 J*=JCCL-1 270 

00 995 K* I OUT *.J 271 

NTERMK )*NTERM(K+1) 272 

EEAMK) * ZEANIK + 1) 273 

SUMX1N) * SUMXIK+l) 274 

SUMX200 => SUMX2TK + 1) 275 

IDCUT(K) =s ICOUT ( K + 1 ) 276 

SUMXY(K,1)= SUMX Y ( K + 1 » 1 ) 277 

995 CONTINUE 278 

I F ( NGTERM . EQ « 1 ) GO TS 10C0 279 

GO TC 65C0 280 

1000 HRITEILIST ,1005) 281 

NOT ERM = C 282 

GC TC 1 C 3 5 283 

C 284 

q*** ***************************************************** *************** 285 

C 286 

1C2C «F(.6CT.IFCH ) GO TO 1C3E 287 

103C WRITEIL I ST , 2 € 1 ) ICEN1 288 

CALI CHISC 289 

C 290 

C * ** ****** $ * ********************************************* *************** 291 

1035 FEAU INPUTS ,1 17) PRECCT 292 

IF! .FCT.PRECCT) GO TC ICC 293 

CALL FRECIC 294 

1C4C GC TC ICC 295 

C 296 

C*44* ****** ********************************************** *************** 297 

2001 FR IT ML 1ST » 1506) 298 

STCF 299 

C********** ************************************************************* 300 

8CC1 FCPMATUHI 301 

8002 FCRMAT ( 1H2 ) 302 

11C fCRMAl ( 12 , 13A6) 303 

111 FCRMAT ( 1H 1 , 1 3A6 , A 2 ) 304 

117 FORMAT I 9L 1 , 1 1 ) 305 

118 FORMAT ( 1H 9Ll»Il) 306 

170 FORMAT ( 33H THERE IS SO BO TERM IN THE MOCEL) 307 

190 FORMAT ( 26H THERE IS A BO TO ESTIMATE > 308 

230 FORMAT ( 401 2) 309 

235 FORMAT! 11H NTERM ( K ) = / ( 1H 3014)) 310 

240 FORMAT! 25H THIE TRANSFORMATIONS ARE /( 1H 5(414, 5X) )) 311 

257 FORMAT! 1L1, F3..3) 312 

260 FORMAT! 5E15. 7) 313 

262 FORMAT! 19H THE CONSTANTS ARE /((llH 8G15i7))l 314 

282 FORMAT ( 2014 ) 315 

283 FORMAT! 1H 2014) 316 
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300 FORMAT! 13A6, 1A2) 317 

301 FORMAT ( 1H 13A6,A2) 318 

370 FORMAT! 1H0 *29H0B SERVED VAR I ABLE S j WEIGHT = G14.6,6X»15H OBSERVATION 319 

1 = ,15) 320 

380 FORMAT ! IH 9614.6) 321 

460 FORMAT ( 1H , 37HTE RMS OF THE EQUATION, OBSERVATION = ,15) 322 

462 FORMAT! 18HK** REPLICATE SET I 5 , 3X , 100 { 1W* ) ) 323 

4671 FORMAT ( 1H 125(1H*)) 324 

467 FORMAT! 14H 0EP. VAR J 16, 8H SSQ=G14.7,8H SUM^G14.7,8H M 325 

XEAN= G14J71 326 

500 FORMAT! 1M0 *2 1HT3 TAL OF THE WEIGHTS ,F8.0) 327 

510 FORMAT! 1H0*69HTQTAL WEIGHT TOO HIGH, EXCEEDS NUMBER OF OBSERVATION 328 

IS BY CNE PERCENT. > 329 

520 FORMAT! 1H0,69HT3TAL WEIGHT TOO LOW, LESS THAN 95 PERCENT OF NUMBER 330 

1 OF EBSEflVATIONS. ) 331 

540 FORMAT! 1H 8G14.7) 332 

560 FORMAT ( 21H2X TRANSPOSE X MATRIX > 333 

670 FORMAT! 25H2CGRRELATI (BN COEFFICIENTS ) 334 

700 FORMAT! 32H2 ( X TRANSPOSE X) INVERSE MATRIX ) 335 

986 FORMAT ( 39H THE NUMBER OF DEPENDENT VARIABLES WAS 13,83H IT IS BE 336 

XING 5ET TO CNE AND THE REJECTION OPTION EXERCISED CN DEPENDENT VAR 337 

XIABLf 1 ) 338 

1005 FORMAT! 39H THERE IS NO EVIDENCE OF A REGRESSION. / 339 

X 74H USE THE MEAN RESPONSE FOR THE BEST ESTIMATE OF THE DEPEND 340 

XENT VARIABLEiS). ) 341 

1282 FORMAT (314*15) 342 

1283 FORMAT ! 1H 314,151 343 

1306 FORMAT ( 4GH REPLICATE SETS ARE NOT GROUPED PROPERLY ) 344 

ENC 345 


$ IBFTC TRAhSX 


SUBRf LT INE TRANS 

*** # ******************* «*** *********************** *************** 

CCMMCN/B1G/S0MXX I 183e),SUMXX [( 1830), EIG! 1830) ,SUMXY (60,9) , 

X E ( 60, 9 I ,CCRR(1830) 

CC MMEN/M EC / B019)„ C0N<99), ERRMS ( 9 ) , 

X ICENTI 13)| IDOUT! 60 ) , NCON ( 200) i 
X NTERMI6C ) 4 NT RAN (100)4 NXCODilOO K P0CLE0(9), 


XREPV AR ( 9 ) * 


PESMS(9), 


SUMX! 138) 


1 

2 

3 

4 

5 

6 

7 

8 
9 


XSUMX2 1691* 

X ( S9 ) , 

ZEANI69), 

SUMY2! 18) 

10 

COMMCN/ SMALL / 

BYPASS j B ZE RC 

1, DELETE, FIRST, IFCHI, 

IFSSR, 

11 

X I FIT, 

IFWT, 

INPUT, 

I N PUT 5 fc 

INTER, 

12 

X I S T R A T 4 

JC0L , 

KCNNO , 

LENGTH* 

LIST, 

13 

X N ERROR * 

NO0EP, 

NCOB, 


NOTERM * 

14 

X NOVAR, 

NPBEG , 

NRES » 

NTRANS, 

NWHERE, 

15 

X P, 

PREDCT* 

REPS, 

RWTt 


16 

X STORY I % 

STORYC, 

STORYX, 

TOTWT , 

WEIGHT, 

17 

X ERRFXC* ECONMY 




18 

LOGICAL ECONMY 





19 

LOGICAL BYPASS 

, BZERO, 

DELETE 

, IFCHI, 


20 

XIFSSR, IFTT , 

IFWT, 

REPS, 

PREDCT 

, 

21 

XSTORYC, STORYX 

, STORYI 

, FIRST 

tERRFXD 


22 
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DOUBLE PRECISION 

RWT f TO TWT , 

WEIGHT 



COMMON/CNTRS/ 

It 

IBC, 

IC, 

ICOL, 

INEW, 

INOCH, 

I OLD * 

I OUT , 

IR, 

IRC, 

IREP, 

IS, 

ITC, 

J ? 

K, 

KBAR 





C ******* 4*********************************^ **************** ****** ****** 

C THIS SUBROUTINE PERFORMS TRANSFORMATIONS IF THIS OPTION IS 

C REGUESTED. 

C 

c 

C K TRANSFORMATION SET NUMBER. 

C NCON<2*K-l) CONSTANT NUMBER TO USE. 

C NCON ( 2*K ) DERIVED CONSTANT. 

C NTRAN(K) NUMBER OF TRANSFORMATION REQUESTED. 

C NXCOC(K) VARIABLE NUMBER 

C 

80 DO SCO K=1 , NTRAN S 
I=NCON ( 2*K-l) 

IF( I ) 100, ICO, 110 
100 C0NS=0. 0 
GO TC 120 
110 CONS=CON< I ) 

120 I =NXCOD ( K ) 

Y=X( I) 

MTRAN = NTRAN(K) 

IF(MTRAN.LE.O) M TRAN* 3 1 

1 AO GO TCt 150, 160, 170,180, 190,200,210,220,230,240,250,260,270,280,290, 
A300,310,32 0, 330, 340,350,360,370,380,390,400,410,420,430,440,450), 

B MTRAN 

150 CONS- Y+CONS 
GO TC 460 
160 CONS=Y*CONS 
GO TO 460 
170 CONS=CONS/Y 
GO TC 460 
180 CONS = EX P ( Y ) 

GO TC 460 
190 CONS=Y**CONS 
GO TC 460 
200 CCN S= ALOG { Y ) 

GO TC 460 

210 C0NS~ALOG10(Y ) 

GO TC 460 
220 CONS=S I N ( Y ) 

GO TC 460 
230 CCNS=COS ( Y ) 

GO TO 460 

240 C0NS=SIN(3.14159265*IC0NS*Y) ) 

GO TC 460 

250 CONS- CO S( 3. 14159265* f CONS*Y ) ) 

GO TO 460 
260 CONS= 1. 0/ Y 
GO TC 460 

270 CGNS- EXP (CONS /Y) 

GO TO 460 

280 CCNS-EXPI CONS/ ( Y *Y ) ) 

GO TC 460 
290 CON S= SORT ( Y ) 

GO TC 460 

300 CQNS=l*Q/SQRT ( Y 1 
GO TC 460 
310 CONS= CON $** Y 
GO TO 460 
320 CON S= 10 • 0**Y 


23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 
87 
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GO TC 460 

88 

330 CONS=SINH( Y ) 

89 

GO TO 460 

90 

340 CONS-COSHt Y ) 

91 

GO TC 460 

92 

350 CONS^U.O-CO$<Y) )/2.0 

93 

GO TO 460 

94 

360 CONS=AT AN( Y ) 

95 

GO TC 460 

96 

370 C0NS fc ATAN2( Y/CON S ) 

97 

GO TC 460 

98 

380 CON S = Y* Y 

99 

GO TC 460 

100 

390 CONS=Y*Y*Y 

101 

GO TO 460 

102 

400 CGNS-ARS IN ( $QRT ( Y ) ) 

103 

GO TC 460 

104 

410 C0NS=2.0*3.14159265*Y 

105 

GO TG 460 

106 

420 CONS=UO/<2.0*3. 14159265+Y) 

107 

GO TG 460 

108 

430 CONS= ERF ( Y ) 

109 

GO TO 460 

110 

440 CON S-GAMMA ( Y ) 

111 

GO TC 460 

112 

450 CONS^Y 

113 

460 I=NCCN(2*K) 

114 

I F ( I >470,470,480 

115 

470 CON ( K )=CONS 

116 

GO TC 500 

117 

480 CON ( I )=CONS 

118 

IFCI-60) 500,500,490 

119 

490 X( I ) =CONS 

120 

500 CONTINUE 

121 

RETURN 

122 

END 

123 
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$ I BF TC SUMUPX 


C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 

£**** ***************** ********************* ***** 4 ******** *************** 

SUBROUTINE SUMUP 


SUBROUTINE SUMUPS 
PURPOSE 

1 ) CALCULATE (X TRANSPOSE X) AND (X TRANSPOSE Y) MATRICES ONE 
OBSERVATION AT A TIME. 

2 ) COMPUTE TOTAL OF THE WEIGHTS 

** BOTH CALCULATIONS ARE IN DOUBLE PRECISION 


SUBROUTINES 

LOC 


NEEDED 


90 


50 

100 


15 


KB AR= J+NOTERM 

SUMXY ( K , J ) = SUMXY(KtJ) 

CONTINUE 

DO 50 1=1, K 


+ X(K)*X(KBAR)*WEIGHT 


CALL LOCI K , I , IR ) 

SUMXX(IR) = SUMXX(IR) + X ( I ) *X ( K )* WE I GHT 

CCNTINUE 

CONTINUE 

DO 15 J = 1 , NODEP 

KBAR=NOTERM + J 

SUMY 2 ( J ) = SUMY2( J )+X(KBAR )**2*WEIGHT 

TOTWT=TOTWT+WEIGHT 

RETURN 

END 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


COMMON/ BIG/SUMXX { 1830), EIG( 1830 ), SUMXY < 60, 9) ,CORR ( 1830 ) 

15 

DOUBLE PRECISION 

SUMXX, SUMXY 




16 

COMMON/ MED/ 

BO 19) , 


CON ( 99 ) , 

ERRMS ( 9 ) , 

17 

X I CENT < 13), ID0UTC60), NCON( 200) , 



18 

X NTERMI60), NT RAN ( 100 ) , NXCOD(IOO) 

, POOLED ( 9 ) , 


19 

XREPVARI 9) , 



RE SMS ( 9 ) , 

SUMX ( 69 ) , 

20 

XSUMX2 (69 ) , 

X ( 99 ) , 


ZE AN ( 69 ) , 

SUMY2 ( 9 ) 

21 

DOUBLE PRECISION 

SUMX , SUMY 2 




22 

COMMON/SMALL/ 

BYPASS, BZERO 

, DELETE, FIRST, IFCHI, 

IFSSR, 

23 

X I FT T , 

IFWT, 

INPUT, 

INPUT5, 

INTER, 

24 

X ISTRAT, 

JCOL, 

KONNO, 

LENGTH, 

LIST, 

25 

X NERROR, 

NOOEP, 

NOOB, 


NOTERM, 

26 

X NOVAR, 

NPDEG , 

NRES , 

NTRANS, 

NWHERE , 

27 

X P, 

PREDCT, 

REPS, 

RWT , 


28 

X STORY I, 

STORYC, 

STORYX, TOTWT , 

WEIGHT, 

29 

X ERRFXD , ECONMY 




30 

LOGICAL ECONMY 





31 

LOGICAL BYP A 

SS, BZERO, 


DELETE, IFCHI 

« 

32 

XIFSSR, IFTT 

* IFWT, 


REPS, PREDCT, 

33 

XSTORYC, STOR 

YX, STORYI 

t 

FIRST , ERRFXD 


34 

DOUBLE PRECISION 

RWT,TOTWT, WEIGHT 



35 

COMMON/ CNTRS/ 

I, 

IBC, 

IC, 

ICOL, 

36 

X INEW* 

INOCH, 

I OLD, 

I OUT , 

IR, 

37 

X IRC, 

IREP, 

IS, 

ITC, 

Jt 

38 

X K, 

KBAR 




39 

C 





40 

c ********* 4 *********** 

************************************************ 

41 

DO 110 1=1, JCOL 





42 

SUNX ( I ) = SUMX ( I) + 

X( I ) *WE I GHT 




43 

110 CONTINUE 





44 

DO ICO K=1 , NOTER 

M 




45 

DO 90 J=l, NODEP 





46 


47 

48 

49 

50 

51 

52 

53 

54 
5 5 

56 

57 

58 

59 

60 
61 
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SIBFTC MFIXXX 


SUBRCUTINE MFIX 

C THIS ROUTINE USES THE SUMXX MATRIX COMPUTE IN DOUBLE PRECISION 

C AND THE SUMX ARRAY COMPUTED IN DOUBLE PRECISION TO COMPUTE THE 

C SUMXX DEVIATIONS FORM OF X TRANSPOSE X IN DOUBLE PRECISION AND 

C THE RESULT IS TRUNCATED TO SINGLE* 

C THE MATRICES ARE ALSO PRINTED 

C 

£*********************************************************************** 


COMMON/ BIG/SUMXX < 1 8 30 ) * E IG ( 1 830 ) , SUMXY ( 60 , 9 > , CORR ( l 830 ) 
DOUBLE PRECISION SUMXX, SUMXY 


10 


COMMON/ MED/ 


BO (9) , 

CON ( 99 ) . 

ERRMS ( 9 ) , 

11 

X I CENT ( 1 3 ) , 

IDOUT (60) ,NCON< 200) , 



12 

X NTERMI 60 ) « 

NT RAN { 100) , 

NXCODI 100) , 

POOLED ( 9 ) , 


13 

XREPVARI9), 



RE SM S ( 9 ) , 

SUMX ( 69 ) , 

14 

XSUMX2 ( 69 ) , 


X( 30) , 

ZE AN < 69 ) , 

SUMY2 ( 9) 

15 

DOUBLE PRECISION 

SUMX . SUMX2 » SUM Y2 



16 

COMMON /FRMTS/ 

FMTl 13) .FMTTRI ( 14) 



17 

COMMON/SMALL/ 


BYPASS, 

BZERO, DELETE, 

FIRST, IFCHI, 

IFSSR, 

18 

X IFTT, 


IFWT, 

INPUT, 

INPUT5 , 

INTER, 

19 

X ISTRAT, 


JCOL, 

KONNO, 

LENGTH, 

LIST, 

20 

X NERROR, 


NODEP, 

NOOB, 


NOTERM, 

21 

X NOVAR, 


NPDEG , 

NRES , 

NTRANS , 

NWHERE , 

22 

X P, 


PREDCT, 

REPS, 

RWT, 


23 

X STORYI. 


STORYC 

, STORYX, 

TOTWT , 

WEIGHT, 

24 

X ERREXD, 

ECONMY 




25 


LOGICAL ECONMY 

DOUBLE PRECISION RWT , TOT WT , WE I GHT 
LOGICAL BYPASS, BZERO, 

XIFSSR, I FTT , IFWT, 

XSTORYC, STORYX, STORYI, 


DELETE , 
REPS, 


IFCHI, 

PREDCT, 


COMMON/CNTRS/ 
X INEW, 

X IRC, 

X K, 

EQUIVALENCE 


,ERRFXD 

IC, 

I OUT , 
ITC, 


ICOL, 
IRt 
J » 


FIRST 

It IBC, 

INOCH, I OLD , 

IREP, IS, 

KBAR 

( El G , SNGLX ) , I E IG ( 70) ,RIT£XY) , ( SNGLXX , SUMXX ) , ( SNGLXY , 

X SUMXY ) , ( SNGLX2 , SUMX2 ) , < SNGLY2 , SUMY2 ) 

DIMENSION SNGLX ( 69 ) , SNGL XY ( 60 , 9 ) , SNGLXX ( 1 8 30 ) , R I TE X Y ( 60 , 9 ) 

X , SNGL X2 ( 69 ) ,SNGLY2(9> 

C ******** ** *********** ******* * *** * *********************************** *** 

C 

DO 10 J = 1 , JCOL 
10 SNGLX ( J ) = SNGL ( S UMX ( J ) ) 

I F ( ECCNMY ) GO TO 500 
WRITE(LIST,530) 

WRITE (LIST, 540) 

WR I TE ( L 1 ST ,560) 

DO 20 1=1, LENGTH 
20 CORR ( l ) = S NGL ( SUMXX ( [ ) ) 

CALL TR I ANG ( CORR , NOTERM , 8 , FM TTR I ) 

DO 30 1=1, NOTERM 
DO 30 J = 1 , NODEP 

30 R1TEXY ( I , J ) = SNGL ( SUMXY ( I , J ) ) 

WRITE(L 1ST, 565) 

CALL RE CT( NOTERM, NODEP, 60, 9, RITEXY,FMTTR I ) 

C 

C* **************************************************** ****************** 

C COMPUTE ANO PRINT MEANS* COMPUTE AND PRINT THE ( X TRANSPOSE X) 

C MATRIX IN TERMS OF DEVIATIONS FROM MEAN. THE DEVIATIONS FORM 

C OF (X T X) IS THE VARIANCE-COVARIANCE MATRIX OF THE 

C INDEPENDENT VARIABLES. 


( SNGLX ( I ), 1 = 1, JCOL) 


26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
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500 CONTINUE 

RWT=1.0D0/TOTWT 
DO 570 1 = 1, JCOL 
570 ZEAN ( I ) =SUMX < I ) *RWT 
WR ITE(LIST,580) 

WRITE (L 1ST, 5 AO) < ZEAN II ) ,1=1 ,JCOL) 

IR = 0 

DO 600 J = 1 , NOTER M 
IR=IR ♦ J 

I F ( • NOT • BZERO ) GO TO 601 

$UMX2( J )=SUMXX< IR)-$UMX< J)**2 *RWT 

GO TC 600 

601 SUMX2IJ )=SUMXX( IR) 

600 CONTINUE 

602 CONTINUE 
I R= 1 

DO 620 J = 1 » NOTERM 

DO 618 K= 1 , NODEP 

I F ( B ZERO ) GO TO 617 

SNGLX Y ( J , K ) =SNGL ( SUMXYl J,K ) ) 

GO TO 618 

617 KBAR=NOTERM+K 

SNGLX Y( JfK) = SNGL < SUMXY ( J , K 1 - SUMX ( J ) * SUMX ( KBAR ) *RWT ) 

618 CONTINUE 

619 DO 620 K= 1 , J 

I F { • NOT * BZERO ) GO TO 6191 
SUVXX<IR)=SUMXX( IR)-SUMX<K)*SUMX{J )*RWT 
619 1 CORR ( IR ) =SNGL < SUMXX ( I R ) /DSQR T < SUMX2 ( J ) *$UMX2 I K ) ) ) 

6193 SNGLXXI IR)=SNGL< SUMXXUR > ) 

6194 I R= I R + l 

620 CONTINUE 


C*** ******* ********* ** ********** ** I **************************** ******* 

I F ( ECONMY ) GO TO 622 
I F ( • NOT • BZERO ) GO TO 621 
WRITEIL 1ST, 625) 

CALL TRIANGI SNGL XX * NOTERM, 8 » FMT TR I ) 

WR ITE(LIST ,630) 

CALL RECT( NOTERM, NODEP, 60,9, SNGLXY , FMTTR I ) 

621 WRITEIL 1ST, 670) 

CALL TRIANGI CORR , NOTERM , 8 , FMTTR I ) 

622 CONTINUE 

DO 640 J=l, NOTERM 
640 SNGLX2I J )= SNGL ( S UMX 2 < J) ) 

00 650 J = l, NODEP 

1 F ( B ZERO ) GO TO 645 
SNGLY2I J)=SNGL<SUMY2( J) ) 

GO TO 650 

645 K=NCTERM+J 

SNGLY2I J ) = SNGL ( S UMY 2 ( J ) - SUMX (K)**2*RWT ) 

650 CONTINUE 


C *** * * **** * ** ** * ** ** ** *** * **** ******** ** ******* * ************ ************ 



RETURN 

530 

FORMAT 

540 

FORMAT 

560 

FORMAT 

565 

FORMAT 

580 

FORMAT 

625 

FORMAT 

630 

FORMAT 

XMEAN 

670 

FORMAT 

END 


8G15.7) 


OF 


OF INDEP AND DEP VARIABLES ) 

)SE X MATRIX ) 

)SE Y MATRIX ) 

INDEP AND DEP VARIABLES ) 

)SE X MATRIX WHERE X IS DEVIATION FROM MEAN ) 
JSE Y MATRIX WHERE X AND Y ARE DEVIATIONS FROM 
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tIBFTC MATINX 


SUBROUTINE MATINV 

*********************************************************************** 

PURPOSE 

1 ) COMPUTE EIGENVALUES AND EIGENVECTORS OF (X-TRANSPOSE X) 

AND/OR CORRELATION MATRIX IF REQUESTED. 

STORYC =.T RUE . IF FOR CORRELATION 
STORYX = .TRUE IF FOR XTX 10 

2 ) I NVERT CORRELATION COEFFICIENT MATRIX BY EITHER BORDERING 11 

OR GAUSS ELIMINATION 12 

3 ICOMPUTE PRODUCT OF CORRELATION AND INVERTED CORRELATION 13 

MATRIX IF REQUESTED 14 

STORYI = .TRUE. IF THE PRODUCT IS TO BE PRINTED 15 

4 JCOMPUTE (X TRANSPOSE X) INVERSE FROM INVERTEO CORRELATION 16 

MATRIX 17 

5JC0MPUTE REGRESSION COEFFICIENTS 18 

6 ICOMPUTE OTHER REGRESSION STATISTICS 19 

20 

SUBROUTINES NEEDED 21 

BORD 22 

LOC 23 

EIGEN 24 

INVXTX 25 

RECT 26 

RSTATS 27 

TRIANG 28 

29 

REMARKS 30 

THE EIGENVALUES ARE COMPUTFD AS AN AID IN DETERMINING THE 31 

CONDITION OF THE SYSTEM OF EQUATIONS FOR THE REGRESSION 32 

COEFFICIENTS. EXAMINATION OF THEM AND THEIR ASSOCIATED 33 

EIGENVECTORS MAY SHOW THAT CERTAIN SETS OF INDEPENDENT 34 

VARIABLES ARE HIGHLY CORRELATED AND NOT EASILY LIABLE TO 35 

INDEPENDENT STUDY. 36 

37 

SUBROUTINE MATINV 38 

39 

********************************************************************** 40 

COMMON/ BIG/SUMXX ( 18 30 ) , SUMXX I ( 1 830 ) , E l G < 18 30 ) , SUMXY < 60 ,9 ) , 41 

X B<60,9) ,CORR< 1830) 42 

COMMON/MEO/ 80<9), C0N<99), ERRMS (9) , 43 

X IDENTJ 13),ID0UT(60),NC0N(200>, 44 

X N TERM 1 60 ) , NTRAN(IOO), NXCOD(IOO), P00LED(9), 45 

XREPVAR(9), RESMS ( 9) * SUMX(138), 46 

XSUMX2 ( 69 ) , X ( 99 ) , ZEAN (69) , SUMY2(18) 47 

CCMMCN /FRMTS/ FMT ( 1 3 ) , FMTTR I ( 14 ) 48 

COMMON/SMALL/ BYPASS , BZERO, DELETE , FIRST, IFCHI, IFSSR, 49 

X IFTT, I FWT , INPUT, INPUT5, INTER, 50 
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X 

ISTRAT, 

JCOL , 

KONNO, 

LENGTH, 

LIST, 

51 

X 

NERRGR, 

NODEP , 

NOOB, 


NOTERM, 

52 

X 

NOVAR, 

NPDEG , 

NRES , 

NTRANS, 

NWHERE , 

53 

X 

P, 

PREDCT, 

REPS, 

RWT , 


54 

X 

STORYI , 

STORYC, 

STORYX, 

TOTWT , 

WEIGHT, 

55 

X 

ERRFXD, 

ECONMY 




56 


LOGICAL 

LOGICAL 

XIFSSR, 

XSTORYC, 


ECONMY 

BYPASS, 
I FTT , 
ST3RYX, 


BZERO, 
IFWT, 
STORYI , 


DELETE, 
REPS, 
FIRST 


IFCH I , 
PREDCT , 
,ERRF XD 


RWTtTOTWT, WEIGHT 
It IBC, 

INOCH, IOLD, 

IREP, IS, 

KBAR 

DIMENSION A { 1 ) , C ( 1) , XTX ( 3 ) , CMAT ( 3 ) , HOL ( 3 ) 
EQUIVALENCE <SUMXX,A), (SUMXXI,C) 

DATA (XTXII),I*I,3) /6HX TRAN, 6HSP0SE , 
DATA ( CMAT (I), 1=1,31/ 6HC0RREL , 6HATI0N 


DOUBLE PRECISION 
COMMON/ CNTR S / 

X INEW, 

X IRC, 

X K, 


1C, 
IOUT , 
ITC, 


6HX 

tlH 


ICOL, 
IRt 
J t 


O*********** ** ******* ******************** ****** ************************ 

IORDER= NOTERM 
I F ( NCTERM- 1 ) 10, 10,12 

10 SUMXX I ( 1 ) = 1 .0/SUMX2 I 1 ) 

GO TO 350 
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62 

63 

64 
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69 

70 

71 

72 

73 

74 
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76 

TRANSFER CORR TO A FOR INVERSION 77 

AT THIS POINT THE INFORMATION IN A MAY BE USED FOR ANY DESIRED 78 

CALCULATIONS EIGENVALUES , RANK, ETC* 79 

JUST PUT CORR INTO A BEFORE PROCEEDING TO REMAINDER OF ROUTINE 80 

81 

12 BYPASS=. FALSE* 82 

I F ( • NCT • STORYX) GO TO 30 83 

15 DO 14 1-1,3 84 

14 HOL ( I)=XTX< I ) 85 

16 CALL E I GEN ( A , SUMXX I , I ORDER , 0 ) 86 

WRITE! LI ST, 17)(H0L( I), l = 1,3) 87 

J~0 88 

00 18 1-1,1 ORDER 89 

J=J+I 90 

18 A( I )=A< J) 91 

WRITEIL 1ST, 19) I A< I ) , 1=1 , IORDER ) 92 

WR I TE ( L I ST * 20 ) 93 

CALL RECTI IORDER , I ORDER, I ORDER, IORDER , SUMXX I , FMTTR I ) 94 

30 DO 35 1=1, LENGTH 95 

35 A I I )=CORR ( I ) 96 

I F ( B YPA SS ) GO TO 49 97 

IF ( • NCT • STCRYC ) GO TO 49 98 

DO 36 1=1,3 99 

36 HOL I I ) = CM AT ( I ) 100 

BYPASS= . TRUE * 101 

GO TC 16 102 

C 103 

q************ ***************************************************** ****** 104 

C NO SUBMODELS TO ANALYZE SO INVERT A DIRECTLY BY GAUSS 105 

49 IF(IFSSR) GO TO 50 106 
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CALL INVXTX( A, NOTERM, D, 1.0) 107 

GO TO 60 108 

109 

********************* ************************************************** no 

C SUBMODELS HAVE BEEN REQUESTED SO WE USE BORDERING 111 

50 I0RDER=0 112 

55 I ORDER 3 IORDER +1 113 

CALL BORDI IORDER, A) 114 

60 IF ( .NOT .STORY I ) GO TO 200 115 

116 

*********************************************************************** 117 

WRITE INVERSE OF CORRELATION COEFFICIENT MATRIX 118 

WRITEIL 1ST, 65) 119 

CALL TRIANGIA, IORDER, 8, FMTTRI) 120 

COMPUTE A TIMES A INVERSE AND WRITE IT 121 

WRITEIL 1ST, 70) 122 

ITC= 0 123 

DO 150 IC=1, IORDER 124 

DO 150 I R= 1 * IC 125 

I TC= ITC+1 126 

C( ITC ) = 0.0 127 

DO 130 1=1, IORDER 128 

CALL LOCI IR, I, IRC) 129 

CALL LOCI I , I C , I B C ) 130 

C ( I TC ) = C ( ITC) + A I IRC) *CORR (IRC) 131 

130 CONTINUE 132 

150 CONTINUE 133 

134 

CALL TRIANGIC, IORDER, 8, FMTTRJJ 135 

136 

200 CONTINUE 137 

138 

********************* ************************************************** 13g 

COMPUTE SUMXXI FROM CORR INVERSE. SUMXXI TIMES THE ERROR MEAN 140 

SQUARE IS THE VARIANCE-COVARIANCE MATRIX OF THE REGRESSION 141 

COEFFICIENTS. 142 

I R=0 143 

DO 340 1=1, IORDER 144 

DO 340 J=1,I 145 

IR= IR+1 146 

SUMXXI! IR) = A( IR J/SQRTI SUMX2I I )*SUMX2 I J) ) 147 

340 CONTINUE 148 

149 

********************* ************************************************** 150 

COMPUTE COEFFICIENTS ANO PRINT THEM 151 

152 

350 DO 370 J = 1 , NODEP 153 

DO 370 K=l, IORDER . 154 

B(K,J)=0.0 155 

DO 370 L=l , IORDER 156 

CALL LOCI L , K , IR ) 157 

B I K , J ) = B 1 K , J) + SUMXX I I IR )*SUMXYIL»J) 158 

370 CONTINUE 159 

C 160 

WRITEILIST-,380) IDENT 161 

WRITEILIST ,382) 162 
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IF ( .NOT.BZERO) GO TO 400 163 

DO 390 J=l • NODEP 164 

SUM=0.0 165 

KBAR* NOTERM + J 166 

DO 385 K = 1 » IORDER 167 

SUM = SUM + B!K, J)*ZEAN!K) 168 

385 CONTINUE 169 

BO ( J ) * ZEANIKBAR) -SUM 170 

390 CONTINUE 171 

WRITE(LIST,395> 172 

WRITE!LIST,397) I BO (K ) , K* 1, NODEP ) 173 

400 WR I TE ! L l ST, 410) 174 

DO 430 J=1 , IORDER 175 

WR IT E (L 1ST » 432) IDOUT ( J ) , < B! J ,K ) ,K=l , NODEP ) 176 

430 CONTINUE 177 

178 

********* ******* ***** ******************** ****** * ********* ************** 179 

COMPUTE REGRESSION STATISTICS IN RSTATS 180 

181 

CALL RSTATS! IORDER) 182 

183 

*********************************************************************** 184 

IF IORDER IS LESS THAN NOTERM WE HAVE USED THE BORDERING OPTION 185 

AND MUST GO BACK TO FINISH. 186 

187 

I F ( I CRDER-NOTERM ) 55,500,500 188 

500 STORYC=. FALSE. 189 

STORYX=. FALSE. 190 

RETURN 191 

17 FORMAT! 34H2THE FOLLOWING ARE EIGENVALUES OF 2A6,A1, 7H MATRIX) 192 

19 FORMAT! 1H 8G16.7) 193 

20 FORMAT! 132H2THIS IS THE MODAL MATRIX OR MATRIX OF EIGENVECTORS. El 194 

1GENVECT0RS ARE WRITTEN IN COLUMNS LEFT TO RIGHT IN SAME ORDER AS 195 

2E IGENVALUES ) 196 

65 FORMAT! 14H2C0RR INVERSE ) 197 

70 FORMAT! 118H2AS A PARTIAL CHECK ON INVERSION ACCURACY THE (CORR)*!C 198 

XORR INVERSE) MATRIX FOLLOWS. IT SHOULD BE THE IDENTITY MATRIX. ) 199 

380 FORMAT! 1H1,13A6, 1A2) 200 

382 FORMAT! 61H EACH COLUMN CONTAINS THE COEFFIC I ENTS .FOR ONE DEPENDEN 201 

XT TERM ) 202 

395 FORMAT! 20H CONSTANT TERM ! BO ) ) 203 

397 F0RMATI4X*9G14. 6) 204 

410 FORMAT! 36H REGRESSION COEFFICIENTS !B1,...,BK) » 205 

432 FORMAT! 1H I3,9G14.6) 206 

END 207 
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11BFTC RSTATX 


SUBRCUT I NE RSTATS 
PURPOSE 

DCOMPUTE AND PRINT THE ANALYSIS OF VARIANCE TABLES ON 
REGRESSION AND LACK-OF-FIT IF APPROPRIATE. 

2 )C OMPUTE AND PRINT R-SQUARED AND STANDARD ERROR OF 
ESTIMATE 

3 JCOMPUTE AND PRINT SUMS OF SQUARES DUE TO EACH VARIABLE 
IF IT WERE LAST TO ENTER REGRESSION 

4 JCOMPUTE AND PRINT THE STANDARD DEVIATIONS OF EACH 
REGRESSION COEFFICIENT. 

SUBROUTINE RSTAT S (I ORDER ) 
c * *********** ****************************** 

COMMON/ BIG/SUMXX ( 1830) , SUMXX I ( 1 830 ) , El G ( 1830 ) * SUMXY < 60 , 9 ) , 


X ERRFXD* ECONMY 



LOGICAL 

ECONMY 



LOGICAL 

SATRTD 



LOGICAL 

BYPASS, 

BZERO, 

DELETE 

XIFSSR, 

I FTT , 

IFWT, 

REPS, 

XSTORYC, 

STOR YX, 

STOR Y I , 

FIRST 


RWT « TOTWT , 

WEIGHT 
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I» 

IBC, 

IC, 

ICOL, 

38 

INOCH, 

I OLD, 

IOUT , 

IR, 

39 

IREP, 

IS, 

ITC, 

J, 

40 

KBAR 




41 


, I FCH I v 

PREDCT, 

,ERRFXD 

DOUBLE PRECISION 
COMMON/ C NT RS/ 

X I NEW , 

X IRC, 

X K, 

C 

£******************************************************** **** *********** 
C PROGRAMMING NOT E*** ******** 

C SOME OF THESE EQUIVALENCES ARE USED TO COMMUNICATE WITH 

C SUBROUTINE PREDCT. BE CAREFUL ABOUT CHANGES INVOLVING THE 

C ARRAY EIG . 

DIMENSION SSQREG ( 9 ) , SSQRES(9) f REGMS ( 9 ) , 

X XLOF ( 9 ) , XLOFMS ( 9 ) , FRATI0<9), RSQD(9), R(9) t 

X S SQL ST ( 9) , DEVB<60,9) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


X 

B(60,9) ,CORR( 1830) 




17 


COMMON/MED/ 

BO ( 9 ) , 

CON ( 99 ) , 

ERRMS < 9 ) , 

18 

X 

IDENTi 13) , 

IDO UT ( 60 ) , NCON ( 200 ) , 



19 

X 

NTERM ( 60 ) , 

NTRAN( 100) , 

NXCOD( 100) > 

POOLED ( 9 ) , 


20 

XREPVAR ( 9 ) » 


RE SMS ( 9 ) , 

SUMX( 138) , 

21 

XSUMX2 ( 69 ) » 

X ( 99 ) , 

ZE AN ( 69 ) , 

SUMY2< 18) 

22 


COMMON /FRMTS/ FMT ( 13 ) , FMTTR 1(14 ) 



23 

COMMON/ SMALL/ 

BYPASS, 

BZERO, DELETE, 

FIRST, I FCH I , 

IFSSR, 

24 

X 

I FTT , 

IFWT, 

INPUT, 

INPUT5 , 

INTER, 

25 

X 

ISTRAT, 

JCOL, 

KONNO, 

LENGTH, 

LIST, 

26 

X 

NERRCR , 

NODEP , 

NOOB, 


NOTERM, 

27 

X 

NOVAR, 

NPDEG, 

NRES , 

NTRANS , 

NWHERE , 

28 

X 

Pt 

PREDCT, 

REPS, 

RWT, 


29 

X 

STORY I , 

STORYC 

, STORYX, 

TOTWT, 

WEIGHT, 

30 


31 

32 

33 

34 

35 

36 


42 

43 

44 

45 

46 

47 

48 

49 

50 
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C 


EQUIVALENCE 

X 

X 

X 

X 


( E I G ( 1) ,SSQREG), 
i EIG ( 19) tREGMS) , 
(EIG(46),XL0FMS) , 
( EIG ( 64 ) *RSQD ) , 
(EIG (82) tSSQLST) , 


(EIG( 10) * SSQRES ) , 
( E I G ( 37 ) « XLOF ) , 

(E IG ( 55 ) , FR ATI 0 ) « 
( E I G ( 73 ) , R ) , 
(EIG(9l) » DEVB) 


C ***************************** ******* ******* 4c * ************************** 


C COMPUTE DEGREES OF FREEDOM AND RECIPROCALS 

NREG = I ORDER 
NTOT = I F IX ( TOTWT )- 1 
IF( .NCT.BZERO) NT0T= NTOT+1 
NRES= NTOT-NREG 
NL0F= NRES - NPDEG 
RNREG= l.O/FLOAT (NREG) 

IF(NRES.EQ.O) GO TO 980 
RNRES= l.O/FLOAT (NRES) 

SATRTC=. FALSE. 

IF(NLCF.EQ.O) GO TO 90 
RNL0F=1 .0/ FLOAT ( NLOF ) 

GO TO 100 
90 SATRTC= .TRUE. 

100 continOe 

NXTERM= IORDER 
RNOOB=R WT 
C 

C*** ********* ********* ************************************************** 

COMPUTE RESIDUAL SUM OF SQUARES, RESIDUAL VARIANCE, VARIANCE 
FROM REPLICATIONS IF APPROPRIATE, AND THE F-RATIO OF MEAN SQUARE 
LACK-OF-FIT AND MEAN SQUARE RESIDUALS. 

DO 210 J= 1 , NODEP 
SSQREGl J ) =0. 0 
DO 2C0 1*1 , NXTERM 

SSQREG( J )=SSQREG ( J ) + B ( I , J ) *$UMXY ( I , J ) 

200 CONTINUE 

SSQRES( J ) =SUMY2( J ) -SSQREG( J ) 

REGMS ( J ) = SSQREG ( J ) * RNREG 
RESMS ( J ) = SSQRES (J )*RNRES 
RSQD( J)*SSQREG( J ) /SUMY2 ( J ) 

R ( J ) =SQRT ( RSQD( J ) ) 

I F ( ( .NOT .REPS).OR.SATRTD) GO TO 210 
XLOF(J)=SSQRES(J ) -POOLED ( J ) 

XLGFMS ( J )= XLOF ( J ) *RNLOF 
FRATIO( J)*XLOFMS (J)/REPVAR( J ) 

210 CONTINUE 


*********************************************************************** 
DETERMINE WHICH ESTIMATE OF SIGMA SQUARED SHOULD BE USED IN 
HYPOTHESIS TESTS. PUT THE PROPER ONE IN ERRMS AND SET ERRFXD 
TO TRUE IF THE PRESENT VALUE IS TO BE USED FOR ALL FOLLOWING 
TESTS AND T-STA T I ST ICS . 

I QUT= I STR AT 
I F ( ERRFXD ) GO TO 250 
IF( ISTRAT.NE.3) GO TO 214 
211 DO 213 J=l, NODEP 
213 ERRMS ( J )= RESMS< J) 

NERRGR = NRES 


51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 
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I OUT = 3 107 

GO TO 250 108 

214 IF( ISTRAT.NE.l) GU TO 218 109 

IF(. NOT. REPS) GO TO 211 110 

DO 215 J=1 »NODEP 111 

215 ERRMS I J ) = REPVARIJ) 112 

NERROR= NPDEG 113 

ERRFXD= .TRUE. 114 

I0UT=1 115 

GO TC 250 116 

218 IFIFIRST. AND. I IDRDER.EQ. NOTERM) ) GO TO 220 117 

GO TO 211 118 

220 ERRFXD= .TRUE. 119 

DO 222 J = 1 1 NODEP 120 

222 ERRMS ( J ) = RESMS(J) 121 

NERRCR= NRES 122 

I 0UT=2 123 

124 

*********************************************************************** 125 

WRITE ANOVA TABLES 126 

250 I S=2 127 

IF(REPS) I S=4 128 

DO 500 J=l, NODEP 129 

IF I ERRM S I J ) .EQ.O.O) ERRM S ( J ) = 1 . 0E-30 130 

WR I T E ( L I ST * 1001 ) IS,J 131 

WRITEIL 1ST, 1002) 132 

WRITEIL 1ST, 1003) SSQREG(J), NREG, REGMS(J) 133 

WRITE(LIST»1004) SSQRES(J), NRES, RESMS(J) 134 

WRITE(LIST,1005) 135 

WRITEIL 1ST, 1006) SUMY2 I J ) , NTOT 136 

WRITEIL 1ST, 1007) 137 

WRITEIL 1ST, 1500) RSODIJ), RIJ) 138 

ST D= SORT I RE SMS I J ) ) 139 

WRITEILIST, 1600) STD 140 

WRITEIL 1ST , 1700) I OUT i ERRMS I J ) , NERROR 141 

F=REGMS(J)/ERRMS(J) 142 

WRITEILIST, 1750JF, NREG, NERROR 143 

IF! I .NOT.REPS l.OR.SATRTD ) GO TO 500 144 

WRITEILIST, 2001) 145 

WRITEILIST, 1002) 146 

WRITEIL 1ST, 2005) XLOFIJ), NLOF, XLOFMSIJ) 147 

WRITEILIST, 2006) POOLEDIJ), NPDEG, REPVARIJ) 148 

WRITEI L 1ST, 1004 ) SSQRE S I J ) , NRE S ,RE SMS I J ) 149 

WRITEILIST, 1005) 150 

WRITEIL 1ST, 2008) FRATIOIJ) 151 

WRITEIL 1ST, 1007) 152 

500 CONTINUE 153 

154 

****************** ******************************************* ********** 155 

COMPUTE CONTRIBUTION OF EACH INDEPENDENT VARIABLE TO REG SUM 156 

OF SQUARES AS IF IT WERE LAST TO ENTER 157 

WRITEILIST, 370) 158 

I R= 0 159 

DO 8635 K= 1 , NXTE RM 160 

I R= IR+K 161 

DO 8632 J=l, NODEP 162 
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8632 SSQL ST ( J ) = 8 ( K, J ) **2/SUMXX I ( IR ) 163 

WRITE(LIST,380) I DOUT ( K ) , ( SSQLST ( J ) , J= l , NODEP > 164 

8635 CONTINUE 165 

C 166 

C * ** * ****** ** *** ****** *********************************************** *** 167 

C COMPUTE STANDARD DEVIATION OF REGRESSION COEFFICIENTS 168 

WRITE(L IST.375) 169 

IF! .NOT.BZERO) GO TO 959 170 

DO 910 J = l ,NXTER M 171 

R ( J ) =0. 0 172 

DO 910 1 = 1 , NXTER M 173 

CALL LOC (IfJ,IR) 174 

R( J)=RI J)+ZEAN( I )*SUMXXI( IR) 175 

910 CONTINUE 176 

XX T = 0 . 0 177 

DO 920 J=1,NXTERM 178 

920 XXT=XXT+ZEANI J)*R( J ) 179 

DO 930 K=l, NODEP 180 

930 DEVB( 1»K)=SQRT(ERRMS(K)*( RNOOB+XXT ) ) 181 

K=0 182 

WRITEU 1ST, 380) K ,( DEVB ( 1 , J ) ,J= 1 .NODEP ) 183 

959 I R=0 184 

DO 970 J=1,NXTERM 185 

I R= IR+J 186 

DO 960 K=l, NODEP 187 

DEVB(J.K) =SQRTI ERRMS(K)*SUMXXI ( IR) ) 188 

960 CONTINUE 189 

WRITE! L 1ST, 380) I DOUT ( J ),( DEVB ( J , KR ), KR = 1 .NODEP ) 190 

970 CONTINUE 191 

C 192 

C ******* «** ** * ** ****** ******** **** ************************************** 193 

C FORMATS 194 

1001 FORMAT! II, 42H ANOVA OF REGRESSION ON DEPENDENT VARIABLE 15) 195 

1002 FORMAT! 1H 79I1H*)/79H SOURCE SUMS OF SQUARES DEG 196 

XREES OF FREEDOM MEAN SQUARES /1H 79!1H-)) 197 

1003 FORMAT I 1 7H REGRESSION G20.8, 5X , 1 10, 5X , G20. 8 ) 198 

1004 FORMAT ( 17H RESIDUAL G20-8, 5X , I 10 , 5X , G20. 8 ) 199 

1005 FORMAT! 1H 79 ( 1H- ) ) 200 

1006 FORMAT! 17H TOTAL G20.8, 5X.I10) 201 

1007 FORMAT! 1H 79!1H*)) 202 

2001 FORMAT! IX/ 1X/22H ANOVA OF LACK OF FIT » 203 

2005 FORMAT! 17H LACK OF FIT G20.8, 5X , I 10 , 5X , G20 . 8 ) 204 

2006 FORMAT! 17H REPLICATION G20.8, 5X,I 10, 5X, G20. 8) 205 

2008 F0RMAT128H F = MS ! LOF ) /MS! REPS ) = F10.3 ) 206 

370 FORMAT! 74H1 SUMS OF SQUARES DUE TO EACH VARIABLE IF IT WERE LAST T 207 

XO ENTER REGRESSION ) 208 

375 FORMAT! U5H2 STANDARD DEVIATION OF REGRESSION COEFFICIENTS (DERIVE 209 

XD FRCM DIAGONAL ELEMENTS OF (X TRANSPOSE X) INVERSE MATRIX )) 210 

380 FORMAT! 1H I3.9G14.6) 211 

1500 FORMAT ( 40H R SQUARED = SSQ(REG) / SSQ(TOT) = F8.6, 212 

X 5 X , 4HR = F7.6) 213 

1600 FORMAT! 34H STANDARD ERROR OF ESTIMATE G14.6) 214 

1700 FORMAT! 24H USING POOLING STRATEGY I2.25H THE ERROR MEAN SQUARE = 215 

X G14.7, 26H WITH DEGREES OF FREEDOM = 16) 216 

1750 FORMAT! 5X, 19HF=M S ( REG ) /M S( ERR ) = F6.2.5X, 13HC0MP ARE TO F(I2,1H,I3,1 217 

XH)) 218 

RETURN 219 

980 WRITEIL 1ST, 981) 220 

981 FORMAT! 41H ZERO RESIDUAL DEGREES OF FREEDOM. STOP. ) 221 

STOP 222 

END 223 
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$ I BFTC TTEST 


C ********************* ************************************************** 

C 

C SUBROUTINE TTEST 

C 

C PURPOSE 

C COMPUTE THE T-STATISTICS FOR EACH REGRESSION TERM ANO 

C ITS TWO TAILED SIGNIFICANCE LEVEL. THEN DETERMINE THE 

C TERM WITH THE LEAST SIGNIFICANCE AND RETURN THIS 

C INFORMATION TO RAPIER. 

C 

C* *************************************************************** ******* 
SUBROUTINE TTEST I * ) 

C 

C*** ****************** *4**4*****$**************$************************ 
COMMON/ B I G/SUMXX ( L 830 ) , SUMXX I ( 1 830 ) , E I G I 18 30 ) , SUMX Y ( 60 , 9 > , 

X B ( 60 y 9 ) , CORR ( 1830) 

COMMON/ MED/ B0(9), C0N(99) f ERRMS 1 9 ) f 

X I DENT ( 13 ) , IDOUTI 60),NC0NI 200) , 

X NTERMI60), NTRAN(IOO), NXC0DI100), P00LED(9) f 

XREPV AR ( 9 ) , RESMSI9), SUMXI138), 

XSUMX 2 ( 69 ) , X(99), ZEAN(69), SUMY2I18) 

COMMON/SMALL/ BY P AS S , B ZE RO , DELE TE , FIRST, IFCHI, IFSSR, 

X I FTT , IFWT, INPUT, INPUTS , INTER, 

X ISTRAT, JCOL , KONNO, LENGTH, LIST, 

X N ERROR , NODEP, NOOB , NOTERM, 

X NOVAR, NPDEG , NRES, NTRANS, NWHERE , 

X P, PREDCT , REPS, RWT , 

X STORYl, STORYC, STORYX, TOTWT , WEIGHT, 

X ERRFXD , ECONMY 

LOGICAL ECONMY 

LOGICAL BYPASS, BZERO, DELETE, IFCHI, 

XIFSSR, I FTT , IFWT, REPS, PREDCT, 

X STORYC , STORYX, STORYI, FIRST , ERRFXD 

DOUBLE PRECISION WEIGHT, RWT 

COMMON/C NT RS/ I, IBC, IC, ICOL, 

X INEW, INOCH, IOLO, IOUT, IR, 

X IRC, IREP, IS, ITC, J, 

X K, KBAR 

C 

Q*** ******* *********** ************************************************** 

LOGICAL MAKENU, NOZERO 
DIMENSION T(35, 13),PLEVEL( 13) 

DIMENSION DEVB ( 60,9) , PR0B(60,9) ,TT<60,9) 

EQUIVALENCE ( E I G ( 9 1 ) , DE VB , TT ) , (EIGI650) ,PROB) 

EQUIVALENCE (P,PWANT) 

C 

C *********************************************************** ************ 

DATA (PLEVELIJJ ) , J J-I, 13 ) / 0. I 0 , 0. 20 , 0 . 30 , 0 . AO , 0. 50 , 0 . 60 , 0 . 70 , 
10.80,0.90,0.95,0.98,0.99,0.999 / 

DATA (T(lfJJ) , J J= 1 , 1 3 ) /O. 15 8, 0.325, 0.5 10,0. 727, 1.000,1.376, 


1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

1A 

15 

16 

17 

18 

19 

20 
21 
22 
23 

2 A 

25 

26 

27 

28 

29 

30 

31 

32 

33 

3 A 

35 

36 

37 

38 

39 
AO 
A1 
A2 
A3 
AA 
A5 
A6 
A7 
A8 
A9 
50 
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1.963,3.078,6.314,12.706,31.821,63.657,636.619 /, 51 

(T(2,JJ) ,JJ=1,13) /0. 142, 0.289, 0.445, 0.617, 0.816, 1.061, 52 

1.386,1.886,2.920,4.3027,6.965,9.925,31.598 /, 53 

(T(3,JJ) *JJ=1»13) /0. 137, 0.277, 0.424, 0.584, 0.765, 0.978, 54 

1.250,1.638,2.353,3.1825,4.541,5.841,12.924 /, 55 

(T(4,JJ) ,JJ=1,13) /0. 134, 0.271, 0.414, 0.569, 0.741, 0.941, 56 

1.190,1.533,2.132,2.7764,3.747,4.604,8.610 /, 57 

(T(5,JJ) , JJ=1,13) /0. 132, 0.267, 0.408, 0.559, 0.727, 0.920, 58 

1.156,1.476,2.015,2.5706,3.365,4.032,6.869 /, 59 

(T(6,JJ),JJ=1,13) /0. 131, 0.265, 0.404, 0.553, 0.718, 0.906, 60 

1.134,1.440,1.943,2.4469,3.143,3.707,5.959 /, 61 

(T(7,JJ) * JJ= l , 1 3 ) /0. 130, 0.263, 0.402, 0.549, 0.711, 0.896, 62 

1.119,1.415,1.895,2.3646,2.998,3.499,5.408 /, 63 

(T(8,JJ) ,JJ=1,13) /0. 130, 0.262, 0.399, 0.546, 0.706, 0.889, 64 

1.108,1.397,1.860,2.3060,2.896,3.355,5.041 /, 65 

(T(9,JJ) , JJ= 1, 1 3 ) /0. 129, 0.261, 0.398, 0.543, 0.703, 0.883, 66 

1.100,1.383,1.833,2.2622,2.821,3.250,4.781 /, 67 

(T( 10, JJ) ,JJ = 1,13) /0. 129, 0.260, 0.397, 0.542, 0.700, 0.879, 68 

1.093,1.372,1.812,2.2281,2.764,3.169,4.587 / 69 

DATA IT( 11, JJ) ,JJ=1,13) /O. 129, 0.260, 0.396, 0.540, 0.697, 0.876, 70 

1.088,1.363,1.796,2.2010,2.718,3.106,4.437 /, 71 

ITl 12, JJ! , JJ=1,13) /O. 128, 0.259, 0.395, 0.539, 0.695, 0.873, 72 

1.083,1.356,1.782,2.1788,2.681,3.055,4.318 /, 73 

( T< 13, JJ ) , JJ=1, 13) /O. 128, 0.259, 0.394, 0.538, 0.694, 0.870, 74 

1.079,1.350,1.771,2.1604,2.650,3.012,4.221 /, 75 

(T( 14, JJ) , JJ=1, 13 ) /O. 128, 0.258, 0.393, 0.537, 0.692, 0.868, 76 

1.076,1.345,1.761,2.1448,2.624,2.977,4.140 /, 77 

(T( 15, JJ) ,JJ=1,13) /O. 128, 0.258, 0.393, 0.536, 0.691 ,0.866, 78 

1.074,1.341,1.753,2.1315,2.602,2.947,4.073 /, 79 

(T( 16, JJ) ,JJ=1»13) /O. 128, 0.258, 0.392, 0.535, 0.690, 0.865, 80 

1.071,1.377,1.746,2.1199,2.583,2.921,4.015 /, 81 

(T(17,JJ),JJ=1,13) /O. 128, 0.257, 0.392, 0.534, 0.689, 0.863, 82 

1.069,1.333,1.740,2.1098,2.567,2.898,3.965 /, 83 

<T( 18, JJ) ,JJ=l,13) /O. 127, 0.257, 0.392, 0.534, 0.688, 0.862, 84 

1.067,1.330,1.734,2.1009,2.552,2.878,3.922 /, 85 

( T( 19, JJ) , JJ=1, 13) /O. 127, 0.257, 0.391, 0.533, 0.688, 0.861, 86 

1.066,1.328,1.729,2.0930,2.539,2.861,3.883 /, 87 

<T(20,JJ) ,JJ=1,13) /O. 127, 0.257, 0.391, 0.533, 0.687, 0.860, 88 

1.064,1.325,1.725,2.0860,2.528,2.845,3.850 / 89 

DATA (T(21,JJ),JJ=1»13) /O. 127 , 0. 25 7 , 0 . 39 1 , 0 . 5 32 , 0 . 686 , 0. 85 9 , 90 

1.063,1.323,1.721,2.0796,2.518,2.831,3.819 /, 91 

(T(22,JJ) ,JJ=1,13) /O. 127, 0.256, 0.390, 0.532, 0.686, 0.858, 92 

1.061,1.321,1.717,2.0739,2.508,2.819,3.792 /, 93 

(T(23,JJ) ,JJ=1,13) /O. 127, 0.256, 0.390, 0.532, 0.685, 0.858, 94 

1.060,1.319,1.714,2.0687,2.500,2.807,3.767 /, 95 

(T(24»JJ) ,JJ=l,13) /O. 127, 0.256, 0.390,0. 531 ,0.685,0.857, 96 

1.059,1.318,1.711,2.0639,2.492,2.797,3.745 /, 97 

(T( 25, JJ) , JJ=1,13> /O. 127, 0.256, 0.390, 0.531, 0.684,0. 856, 98 

1.058,1.316,1.708,2.0595,2.485,2.787,3.725 /, 99 

(T(26,JJ) ,JJ=1,13) /O. 127, 0.256, 0.390, 0.531, 0.684, 0.856, 100 

1.058,1.315,1.706,2.0555,2.479,2.779,3.707 /, 101 

(T<27, JJ) ,JJ=1,13> /O. 127, 0.256, 0.389, 0.531, 0.684, 0.855, 102 

1.057,1.314,1.703,2.0518,2.473,2.771,3.690 /, 103 

(T(28,JJ) ,JJ=l,13) /O. 127, 0.256, 0.389, 0.530, 0.683, 0.855, 104 

1.056,1.313,1.701,2.0484,2.467,2.763,3.674 /, 105 

(T(29»JJ) ,JJ=1*13) /O. 127, 0.256, 0.389,0. 530, 0.683, 0.854, 106 
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H 1.055, 1.31 1,1.699,2.0452,2.462,2.756,3.659 /, 107 

I l T( 30, JJ ) ,JJ=1,13) /0. 127, 0.256, 0.389, 0.5 30, 0.683, 0.854, 108 

J 1.055,1.310,1.697,2.0423,2.457,2.750,3.646 / 109 

DATA (T(31»JJ) ,JJ=i,13) /O . 126, 0 . 255 , 0. 388 , 0. 529 , 0. 681 , 0 . 85 1 , 110 

1 1.050, 1.303, 1.684,2.0211,2.423,2.704,3.551 /, ill 

2 <T(32, JJ) ,JJ=1,13) /O. 126,0.254,0. 387,0. 527,0.679,0.848, 112 

3 1.046,1.296,1.671,2.0003,2.390,2.660,3.460 /, 113 

4 <TC 33, JJ ) , JJ=1, 13) /O. 126, 0.2 54, 0.386, 0.526, 0.677, 0.845, 114 

5 1.041, 1.289,1.658, 1.9799,2.358,2.617,3.373 /, 115 

6 ( T ( 34 , J J ) ,JJ=1,13) /O. 126, 0.253, 0.385, 0.524,0.674,0.842, 116 

7 1.036,1.282,1.645,1.9600,2.326,2.576,3.291 / 117 

118 

T ( I I , J J ) IS THE T-STATISTIC AT THE TABULATED OEGREES OF FREEDOM 119 

(II) AND AT THE TABULATED PROBABILITY LEVELS (JJ). 120 

I I=DEGREES OF FREEDOM, EXCEPT FOR 121 


C 


11=31 IS FOR 40 DEGREES 


122 

C 


11=32 IS FOR 60 


123 

c 


11=33 IS FOR 120 


124 

c 


11=34 IS FOR INFINITY 


125 

c 




126 

c 


JJ PROBABILITY LEVEL * 

JJ PROBABILITY LEVEL 

127 

c 


1 0.10 * 

8 0.80 

128 

c 


2 0.20 * 

9 0.90 

129 

c 


3 0.30 * 

10 0.95 

130 

c 


4 0.40 * 

11 0.98 

131 

c 


5 0.50 * 

12 0.99 

132 

c 


6 0.60 * 

13 0.999 

133 

c 


7 0.70 * 


134 

c 




135 

c**** 

******************************************************************* 

136 

c 


CALCULATE T STATISTICS 


137 

c 




138 


220 

WRITE (LIST, 230) 


139 


230 

FORMAT! IH0.23HCALCULATED T STATISTICS /75H THE T STATISTICS CAN BE 

140 



I USED TO TEST THE NET REGRESSION 

COEFFICIENTS B( I). ) 

141 



DO 260 J=l, NOTERM 


142 



DO 240 K = l, NODEP 


143 



TT (J,K)=A8S( B(J, K)/DEVB( J»K ) ) 


144 


240 CONTINUE 

WRITE (LIST, 250) ( TT ( J , K ) , K= 1 , NODEP ) 

250 FORMAT ( 1H 9G14.6 ) 

260 CONTINUE 
C 

£***************************************** ****************************** 
NDEG = NERROR 
C 

C **************************************************************** ******* 


C 

C 


SEARCH THE TABLE OF TABULATED DEGREES OF FREEDOM 

MAKENU=. FALSE. 

IF (NDEG -30 )290,2 90 ,300 
290 I I = NCEG 
GO TC <*00 

300 I F ( NCEG-40 ) 3 10, 320, 330 
310 F I NV = 1 • 0/40.0 
FM1INV- 1.0/30.0 


145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 
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no on 


MAKENU= . TRUE • 163 

320 11=31 164 

GO TO 400 165 

330 IF(NCEG-60)340,350,360 166 

340 FINV=1. 0/60.0 167 

FM1INV=1. 0/40.0 168 

MAKENU= .TRUE . 169 

350 11=32 170 

GO TO 400 171 

360 IF(NCEG-120)370, 380,390 172 

370 FINV=1. 0/120.0 173 

FM1INV=1. 0/60.0 174 

MAKENU= .TRUE . 175 

380 11=33 176 

GO TO 400 177 

390 11=34 178 

F I NV=0. 0 179 

FM1INV=1. 0/120.0 180 

MAKENU= .TRUE. 181 

182 

183 

400 WR1TEIL I ST , 4 10 > 184 

410 FORMAT t 1 04H UNDER NULL HYPOTHESIS THE INTERVAL <-T,T) WHERE T IS G 185 

XIVEN ABOVE, HAS APPROX PROBABILITY LISTED BELOW. /42H MINUS S 186 

XIGN INDICATES PROB EXCEEDS .999. > 187 

IF ( .NOT.MAKENU) GO TO 430 188 

FNDEG=NDEG 189 

DO 420 J J= 1 » 13 190 

T(35,JJ)=T( I I,JJ ) + ( ( l.O/FNDEG - F I NV ) / ( FM 1 1 NV-F I NV ) ) * ( T ( I I - 1 , J J ) 191 

1 -T ( II, JJ) ) 192 

420 CONTINUE 193 

11=35 194 

430 DO 560 J=l, NOTERM 195 

DO 540 K= 1 , NO DEP 196 

DO 440 J J= 1 , 1 3 197 

IFITt II ,JJ)-TT(J,K> >440,450,460 19 8 

440 CONTINUE 199 

PROB(J,K)=-0.999 200 

GO TO 540 201 

450 PROBt J,K)=PLEVEL(JJ) 202 

GO TC 540 203 

460 IFIJJ.LE.9) GO TO 470 204 

J J 1= J J-2 205 

JJ2= J J- 1 206 

JJ3= J J 207 

GO TO 490 208 

470 IFIJJ.LE.4 )G0 TO 480 209 

JJ1=JJ-1 210 

J J 2= J J 211 

J J3= J J+ 1 212 

GO TO 490 213 

480 JJ1=JJ 214 

JJ2= JJ+ 1 215 

J J3= J J+2 216 

217 

PERFORM A THREE-POINT LAGRANGE INTERPOLATION 218 
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C 219 

490 X- ALCG ( T T ( J , K )) 220 

X 1 = ALOG (T(IlfJJl)) 221 

X2 = AL0G ( T < I l , JJ2 ) ) 222 

X3 = ALCG ( T ( 1 1 » JJ3 ) ) 223 

IF(TTU f K).LE.i.O> GO TO 500 224 

Y1=AL0G( i.0-PLEVEL< JJ1) ) 225 

Y2= A LOG ( 1 • O-PLEV EL ( J J2 ) ) 226 

Y3= A LOG ( 1*0~PLEVEL< JJ3) ) 227 

GO TO 510 228 

500 Y1=AL0G(PLEVEL< JJ1) ) 229 

Y2=ALGG(PLEVEL( JJ2) ) 230 

Y3=AL0G(PLEVEL(JJ3) ) 231 

510 PROB ( J * K ) - ( (X-X2)*(X-X3)*Yl)/< ( X 1-X2 ) * < X1-X3 ) ) + ( ( X- X i ) * ( X-X3 ) 232 

1 *Y2)/( <X2-X1)*<X2-X3> ) + < < X-Xl )* < X-X2 ) *Y3 > / < < X3-X1 ) * < X3-X2 ) > 233 

I F ( TT ( J * K ) — l • 0 ) 520,520,530 234 

520 PROB< J,K)=EXP(PROB( J,K) ) 235 

GO TO 540 236 

530 PROB( J»K)=1.0-EXP(PR0B< J,K>) 237 

540 CONTINUE 238 

q*** ****************** ************************************************** 239 

C WRITE THE PROBABILITIES ( 1.0— ALPHA ) 240 

WRITEU 1ST, 5 50) IDOUT ( J ) , < PROB ( J ,K ) , K=1 , NODE P ) 241 

550 FORMAT ( 1H I 3 , 9( 8 X , F 6 . 3 ) ) 242 

560 CONTINUE 243 

244 

*********************************************************************** 245 

LIST THE DESIRED VALUE OF PROBABILITY ( PWANT ) 246 

247 

570 PERCEN=PWANT*100.0 248 

WRITE(LIST ,580) PERCEN 249 

580 FORMAT! 1H0,36HTHE DESIRED VALUE OF PROBABILITY IS ,F5.1, 8H PERCEN 250 

IT ) 251 

252 

DELETE THE TERM WITH THE LOWEST COMPUTED PROBABILITY IF THAT 253 

PROBABILITY IS LESS THAN THAT DESIRED (PWANT) 254 

255 

I F ( • NOT • DELETE ) GO TO 660 256 

I OUT=0 257 

590 AM I N=PW ANT 258 

DO 620 J= 1 , NOTERM 259 

I F ( ABS ( PROB ( J , 1 ) )-PW ANT ) 600,620,620 260 

600 IF(ABS(PROB< J,l) ) -AM IN ) 6 10, 620 , 620 261 

610 AMIN=ABS(PROB(J, 1) ) 262 

I OUT= J 263 

620 CONTINUE 264 

I F ( I OUT ) 660,660,630 265 

630 WRITEtL 1ST, 650) IDOUTUOUT) 266 

650 FORM AT ( 1H 10X,11HTHE TERM X(,I2,18H) IS BEING DELETED ) 267 

GO TO 670 268 

C ALL VARIABLES REMAINING HAVE BEEN CONCLUDED SIGNIFICANT 269 

660 RETURN 1 270 

670 RETURN 271 

END 272 
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SIBFTC PR EC IX 


SUBROUTINE PREOIC 
PURPOSE 

DREAD INPUT LEVELS OF INDEPENDENT VARIABLES AND COMPUTE 

A PREDICTED RESPONSE FROM THF ESTIMATED REGRESSION EQUATION. 
2 ) C OMPUT E VARIANCE AND STADARD DEVIATION OF THE PREDICTED 
MEAN VALUE AND A SINGLE FURTHER OBSERVATION. 

SUBROUTINES NEEDED 
TRANS 
LOC 

REMARKS 

VALUES FOR DEPENDENT VARIABLES ARE NOT NECCESARY FOR THE 
PREDICTING OF VALUES. HOWEVERt A DUMMY VALUE MAY NEED TO 
BE SUPPLIED IF A ZERO (BLANK) INPUT VALUE WILL CAUSE AN 
IMPOSSIBLE OPERATION TO BE ATTEMPTED DURING THE 
TRANSFORMATIONS. 

**** * ** * * * ** * * * ** * ***** ************************** ************* ********* 


SUBROUTINE PREDIC 

COMMON/ B I G/SUMXX ( 1830 ) , SUMXX I < 1 830 ) , E I G ( 18 30 ) , SUMXY ( 60 , 9 ) , 

X B ( 60 , 9 ) , CORR ( 1830) 

COMMON/MED/ B0(9), C0N(99), ERRMS ( 9) * 

X I CENT { 13) , ID0UT(60) ,NCON( 200) , 

X NTERM ( 60 ) , NT RAN ( 100 ) , NXCOD(IOO), P00LED19), 


XREPVAR( 9 ) t 
XSUMX2 ( 69 ) , 

COMMON /FRMTS/ 
COMMON/ SMALL/ 

X IFTT, 

X ISTRAT, 

X NERROR, 

X NOVAR, 

X P, 

X S TOR Y I , 

X ERRFXD , 


X( 99) , 

FMT ( 13) tFMTTRI ( 14) 
BYPASS, BZERO, DELETE, 


RESMS( 9) ,. 
ZE AN (69 ) , 


IFWT, 
JCOL , 
NODEP , 
NPDEG , 
PREDCT, 
STORYC, 
ECONMY 


INPUT, 

KONNO, 

NOOB, 

NRES, 

REPS, 

STORYX, 


FIRST, IFCHI, 
INPUT5 , 
LENGTH, 

NTRANS, 

RWT , 

TOTWT , 


SUMX ( 138) , 
SUMY2( 18) 

IFSSR, 
INTER, 
LIST, 
NOTERM, 
NWHERE , 

WEIGHT, 


LOGICAL 

LOGICAL 

XIFSSR, 

XSTORYC, 


ECONMY 

BYPASS, 
IFTT , 
STORYX, 


BZERO, 
IFWT, 
STORY I , 


DELETE , 
REPS, 
FIRST 


IFCHI , 
PREDCT, 
, ERRF XD 


WEIGHT, RWT, TOTWT 
I, IBC, 

INOCH, IOLD, 

IREP, IS, 

KBAR 


DOUBLE PRECISION 
COMMON/CNTRS/ 

X INEW, 

X IRC, 

X K, 

DIMENSION YC ALC ( 9 ) , V(60), 

X VARP ( 9 ) , SEEP( 9) 

EQUIVALENCE <YCALC< 1), SUMXX ( 1) ) , 

X ( VARMCl ) ,SUMXX(71) ), (SEEM(l) 
X , ( SE EP ( 1 ) , SUMXX ( 98 ) ) 

EQUIVALENCE ( RNOOB , RWT ) 


IC, 
IOUT , 
ITC, 


VA RM ( 9 ) i 


ICOL. 

IR, 

J t 

SEEM ( 9 ) , 


( V( 1) ,SUMXX( 10) ) , 

, SUMXX (80 ) ) , ( VARP ( 1 ), SUMXX (89) ) 


C* ** * ***** * ** ** * ** **** *** * *** * **** *** ***** ******* * * ************* ******** 

C 

IF (NOTERM. EQ.O) RETURN 
WRITE<6,3) 

RE AD ( 5 , 5 ) NPRED 
C 

DO 500 KK*1, NPRED 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
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c 62 

105 READ ( 5» FMT ) ( XU » , 1= l , ICOL ) 63 

WRITEI6, 110) (XI I ), 1=1, ICOL ) 64 

125 CALL TRANS 65 

DO 130 K= 1 , JCOL 66 

I = NT ERM ( K ) 67 

X ( K 1 = CON ( I ) 68 

130 CONTINUE 69 

HR I TE ( 6 » 135 ) IXI I ), 1=1, NOTERM) 70 

71 

COMPUTE PREDICTED RESPONSE 72 

140 DO 150 K= 1 , NODEP 73 

YCALCIK) = BOCK) 74 

IF( .NGT.BZERO) YCALC!K)=0.0 75 

DO 150 J=l, NOTERM 76 

YCALCIK ) = YCALCIK) + B(J,K)*X(J) 77 

150 CONTINUE 78 

79 

COMPUTE VARIANCE AND STANDARD DEVIATION OF REGRESSION LINE 80 

AND VARIANCE AND STANDARD DEVIATION OF PREDICTED VALUE 81 

AT THE POINT XO 82 

83 

DO 250 K=l, NOTERM 84 

V I K ) =0. 0 85 

DO 250 J=l, NOTERM 86 

CALL LOCI J,K, IR) 87 

V ( K ) =V I K ) + I XI J )-ZEANI J) )*SUMXXI I IR) 88 

250 CONTINUE 89 

XXT=0. 0 90 

DO 275 K=l, NOTERM 91 

XXT = XXT + IXIK )-ZEANIK) ) * V I K ) 92 

275 CONTINUE 93 

XRNOCB = RNOOB 94 

I F I . NOT . BZ ERO ) XRN0OB = 0.0 95 

DO 300 K=l, NODEP 96 

VARM I K ) = ERRMSIK )*( XRNOOB + XXT) 97 

SEEM(K)=SQRTIVARM(K) ) 98 

VARP ( K ) = ERRMSIK )+VARM(K) 99 

SEEP(K)=SQRTIVARP(K) ) 100 

300 CONTINUE 101 

HR ITEI6, 310)1 YCALC IK ) ,K=1, NODEP ) 102 

WRITE 1 6, 320) IVARMIK) ,K=1, NODEP) 103 

WRITE 1 6, 320) I SEE M I K ) , K= 1 , NODEP ) 104 

WRITE(6»320) I VA RP I K ) , K= 1 , NODEP ) 105 

WRITEI6,3?0) (SEEPIK) ,K=1, NODEP) 106 

C 107 

C 108 

500 CONTINUE 109 

RETURN 110 

3 FORMAT! 54H1F0R EACH SET OF INDEP VARIABLES THERE IS COMPUTED... / 111 

X 20H PREDICTED RESPONSE / 112 

X 29H VARIANCE OF REGRESSION LINE / 113 

X 34H STANDARD DEVIATION OF REGRESSION / 114 

X 29H VARIANCE OF PREDICTED VALUE / 115 

X 39H STANDARD DEVIATION OF PREDICTED VALUE ) 116 

5 FORMAT! 14) 117 

110 FORMAT! 39HKINPUT DATA FOR THIS PREDICTED RESPONSE /IlH 9G14.6M 118 

135 FORMAT! 56HK INDEPENDENT TERMS ACCORDING TO FINAL REGRESSION MODEL 119 

X /IlH 9G14.6) ) 120 

310 FORM A T I 55HKPRED ICTED RESPONSE FOR ABOVE INDEP VARIABLES 121 

X /<1H 9G14.6)) 122 

320 FORMAT! 1H 9G14.6) 123 

END 124 
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SIBFTC CHISQX LIST 


C*** + ***************** **** * * ** ********** ********** **** *** * ** ** **** * ***** 

C 

C SUBROUTINE CHISQ 

C 

C PURPOSE 

c 

C DCOMPUTE PREDICTED VALUE DF DEPENDENT VARIABLES AND RESIDUALS 

C AT INPUT DATA POINTS. 

C 2) COMPUTE STANDARDIZED RESIDUALS. 

C 3 ) COMPUTE SKEWNESS AND KURTOSIS OF SAMPLE DISTRIBUTION OF 

C RESIDUALS* ALSO USE THE SAMPLE DISTRIBUTION TO COMPUTE 

C THE CHI-SQUARE STATISTIC. 

C 4) PRINT HISTOGRAMS OF THE DISTRIBUTION OF RESIDUALS. 

C 

C SUBROUTINES NEEDED 

C HIST 

C 

£***** ************* ***************************************************** 

C 

SUBROUTINE CHISQ 

COMMON/ B I G/SUMXX i 1 830 ) , SUMXX I ( 1 830 ) * E I G ( 18 30 ) , SUMXY ( 60 , 9 ) , 

X B (60,9) , CQRR ( 1830) 

COMMON/MED/ B0<9), C0N(99>, ERRMS ( 9 ) , 

X IDENT( 13) * I DO UT ( 60 ) , NCON ( 200 ) , 

X NTERM ( 60 ) , NT RAN ( 100 ) * NXCODUOO), P00LED(9>, 

XREPV AR ( 9 ) 9 RE SM S ( 9 ) , SUMX<138), 

X SUMX 2 ( 69 ) » X ( 99 ) , ZEAN(69) f SUMY2 ( 18 ) 

COMMON/SMALL/ BYPAS S , BZE RO , DELE T E , FIRST, IFCHI, IFSSR, 

X IFTT, IFWT, INPUT, INPUT 5 , INTER, 

X ISTRAT, JCOL , KONNO, LENGTH, LIST, 

X NERROR, NODEP , NOOB, NOTERM, 

X NOVAR, NPDEG , NRES , NTRANS, NWHERE , 

X P, PREDCT, REPS, RWT, 

X STORY I , STORYC, STORYX, TOTWT , WEIGHT, 

X ERRFXD , ECONMY 

LOGICAL ECONMY 

DOUBLE PRECISION RWT , TD T WT , WE I GHT 

LOGICAL BYPASS, BZERO, DELETE, IFCHI, 

XIFSSR, I FT T , IFWT, REPS, PREDCT, 

XSTORYC, STORYX, STORYI, FIRST , ERRFXD 

COMMON/ C NT RS/ I, IBC, IC, ICOL, 

X INEW, INOCH, I OLD , IOUT, IR, 

X IRC, IREP, IS, ITC, J, 

X K, KBAR 

C 

C* *************************** ******* *************************** ********* 

INTEGER CELLS, PLUS1 

DIMENSION BOUND ( 45 ) , CELLBD (21), CH I ( 9 ) , 

X OBS ( 20 , 9 ) , RCT (212) , RELKUR ( 9 ) , RELSKW ( 9 ) , 

X STDERR ( 9 ) * VAR ( 9 ) , YCALC(9), YD I FR ( 9 ) , 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
2? 

23 

24 

25 

26 
27 
2B 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 
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X Z(9) 



51 

EQUIVALENCE ( CEL LBD , SUMXX ( l ) ) 

, (CHI ,SUMXX(22) ) , 

52 

X (OBS,SUMXX( 31) ) » 

( RCT , SUMXX ( 221 ) ) , 

53 

X (RELKUR, SUMXX(434 

)), (RELSKW 

t SUMXX (443) ) , 

54 

X ( STD ERR,SUMXX(452 

)), (VAR, RESMS), 

55 

X (YCALC, SUMXX (461 ) 

), (YDIFR, 

SUMXX (470 ) ) , 

56 

X ( Z, S UMXX ( 479 ) ) 



57 

ODATA BOUND/. 674 4890 7 

,.43072720 

,.96741604 

58 

I, .31863932 ,.67448907 

,1.15033859 

,.25334708 

59 

2, .52439979 ,.84162001 

,1.28153777 

,.2 1042836 

60 

3, .43072721 , .674 48922 

, .96741746 

,1.38298075 

61 

4,. 18001235 ,.36610621 

,.56594856 

, .79163735 

62 

5,1.06756653,1.46521688 

, .15731067 

,.31863932 

63 

6, .48877614 ,.67448930 

,.88714436 

,1.15034184 

64 

7, 1.53411831, .13971028 

, .28221612 

,.43072722 

65 

8, .58945544 ,.76470731 

, .96741836 

, 1.22062834 

66 

9,1.59323335, .12566134 

,.25334701 

,.38532026 

67 

A, .52440000 ,.67448801 

, .84161868 

,1.03642921 

68 

B, 1.28154233, 1.64490172 / 



69 


• 


70 

JCOL=NOTERM+ NODEP 



71 

NUVAR=N0TERM+1 



72 

BYPASS* .FALSE. 



73 

KOUNT* 0 



74 




75 

t ************************************************************ ********** 

76 

DETERMINE IF SAMPLE SIZE IS LARGE ENOUGH TO 

PERMIT CHI-SQUARE 

77 

CALCULATION. IF SO, DETERMINE 

NUMBER OE CELLS AND CELL BOUNDARIES 

78 

I F ( NERROR-30 ) 110,120,120 



79 

110 BYPASS 55 • TRUE . 



80 

GO TC 125 



81 

120 CELLS=N00B/5 



82 

CELLS =MIN0(C ELLS, 20) 



83 

1= MCD( CELLS ,2) 



84 

IF(I.NE.O) CELLS=CELLS + l 



85 

FCELLS= FLOAT (CELLS ) 



86 

PLUS 1= CELLS + 1 



87 

MINUS 1 = CELLS -1 



88 

NDEGCH = CELLS-3 



89 

IR = CELLS/2-1 



90 

IC=IR*( IR-1 )/2 



91 

I S= I R+2 



92 

00 122 J ”1 » I R 



93 

IC=IC+1 



94 

I BC= I S- J 



95 

I RC= I S+ J 



96 

CELLED ( I BC ) =— BOUND ( IC) 



97 

CELLBDI I RC ) = BOUND ( IC) 



98 

122 CONTINUE 



99 

CELLBD ( 1 ) = — 1 • 0E+ 37 



100 

CELL BC ( PLUS 1 ) =1.0E37 



101 

CELL E0( IS ) = 0.0 



1 02 

DO 124 K=l, NODEP 



103 

CHI (K)=0.0 



104 

DO 124 1=1, CELLS 



105 

OBS ( I , K ) =0 .0 



106 
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no on on 


124 CONTINUE 


107 

C 108 

C* *************************************************** ******************* 109 

INITIALIZE SKEWNESS AND KURTOSIS ARRAYS. COMPUTE STANDARD ERROR 110 

OF ESTIMATE. Ill 

125 DO 130 K=1 .NODEP 112 

RELKUR(K)=0.0 113 

RELSKW ( K ) =0. 0 114 

STDERRI K )= SORT! ERRMSIK >) 115 

130 CONTINUE 116 

WRITEIL 1ST, 135) 117 

135 FORMAT! 51H FOR EACH DEPENDENT TERM AND OBSERVATION IS PRINTED 118 

X /31H OBSERVED RESPONSE (Y OBSERVED) 119 

X /29H CALCULATED RESPONSE (Y CALC) 120 

X /28H RESIDUAL (YOBS- YC ALC=YDI F ) 121 

X /28H STANDARDIZED RESIDUAL ( Z ) ) 122 

123 

*********************************************************************** l2 4 

DO 430 J= 1 » NOOB 125 

READ! INTER ) (X(I >.1 = 1,69 ), WEIGHT 126 

DO 142 1=1 .NOTERM 127 

K= IDOUT(I) 128 

X! I )= X (K) 129 

142 CONTINUE 130 

KB AR=NWHERE 131 

DO 143 1=1, NODEP 132 

IC= NOTERM* I 133 

KB AR=KB AR+ 1 134 

X ( I C ) = X(KBAR) 135 

143 CONTINUE 136 

137 

*********************************************** 

DO 160 K= 1, NODEP 139 

YC ALC ( K ) = BO ( K ) 140 

IF! .NCT.BZERO) YCALC!K)= 0.0 141 

KB AR= K+NOTERM 142 

00 150 1=1, NOTERM 143 

YCALC(K) = YCALC(K) + B(I,K)*X!I) 144 

150 CONTINUE 145 

ACTDEV= X(KBAR)- YC ACC ( K ) 146 

YD I FR ( K ) = ACTDEV 147 

Z(K)=ACTDEV/STDERR(K) 148 

ACTDE3=ACTDEV**3 149 

RELSKW (K)= RELSKW(K)*ACTDE3 150 

RELKUR(K) = RELKUR(K)+ACTDE3*ACTDEV 151 

160 CONTINUE 152 

WRITEIL 1ST, 180) ( X ( K ) , K =NUVAR , JCOL ) 153 

WRITEIL 1ST, 190) ( YC ALC ( K ), K= 1 , NODEP ) 154 

WRITEILIST ,200) ( YD I FR ( K ) , K= 1 , NODEP ) 155 

WRITEIL 1ST, 210) ( Z ( K ), K= 1 , NODEP ) 156 

180 FORMAT! 12HKY OBSERVED ,9G13.4) 157 

190 FORMAT! 12H Y CALC .9G13.4) 158 

200 FORMAT! 12H Y DIF .9G13.4) 159 

210 FORMAT! 12H STUDENTIZED ,9G13.4> 160 

IF ( BYPASS ) GO TO 410 161 

C 162 
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C*** ********* ********************* ************** 163 

DO 250 K= 1 » NODEP 164 

DO 230 1 = 1, PLUS1 165 

IF ( 2 IKJ-CELLBDi I ) ) 220,220,230 166 

220 OBS! I-1,K)=0BS( I-1.KJ+ 1.0 167 

GO TO 250 168 

230 CONTINUE 169 

250 CONTINUE 170 

C 171 

410 KOUNT = KOUNT +1 172 

IFIKOUNT.LT. 10) GO TO 430 173 

WRITE(LIST,270) IDENT 174 

270 FORMAT! 1H113A6.A2) 175 

K0UNT=0 176 

430 CONTINUE 177 

C 178 

C********** ************************************************* ************ 179 

C PRINT SKEWNESS AND KURTOSIS 180 

DO 440 K=l, NODEP 181 

RELSKW!K)=RELSKW<K)**2/(FL0AT!N00B)**2*ERRMS(K)**3) 182 

RELKURIK )=RELKUR ( K ) / I FL OAT ( NOOB ) *ERRMS ( K ) **2 ) 183 

440 CONTINUE 184 

WRITE!! 1ST, 450) IOENT, (RELSKW!K ) ,K=I ,NOOEP) 185 

450 FORMAT! 1H1 1 3A6 // 1 OX , 30HSKE WNESS < SHOULD BE NEAR ZERO) // 186 

X 12X , 9F 12 .4) 187 

WR I TE ( L I ST , 460 ) < RELKUR < K ), K = 1 , NODEP ) 188 

460 FORMAT! 10X, 31HXURT0SIS < SHOULD BE NEAR THREE) // 1 2X , 9F 12 . 4 ) 189 

IF!. NOT. BYPASS) GO TO 480 190 

WRITEfLIST ,470) 191 

470 FORMAT! 74HKCH I-SQUARE IS NOT COMPUTED FOR LESS THAN 30 DEGREES OF 192 

XFREEDOM FOR ERROR. ) 193 

RETURN 194 

C 195 

C**** **************************************************** *************** 196 

C COMPUTE CHI-SOU ARED AND PRINT HISTOGRAMS OF RESIDUALS 197 

480 DO 580 K=l, NODEP 198 

DO 570 1=1, CELLS 199 

CH I ( K )= CHIIK) +OBS ( I , K ) **Z 200 

570 CONTINUE 201 

CHI! K)=FCELLS*CHl (K) /FLOAT (NOOB (-FLOAT! NOOB) 202 

580 CONTINUE 203 

WRITE!! 1ST, 590) NDEGCH, ( CH I ( K ) , K= 1 , NODEP ) 204 

590 FORMAT! 55HKTHE CHI-SQUARED VALUES ARE LISTED BELOW. COMPARE WITH 205 

X I 10, 20H DEGREES OF FREEDOM / 1H 9G14.6) 206 

C 207 

RELFRQ = TOTWT/FCELLS 208 

DO 650 K=l, NODEP 209 

WRITEfLIST, 620) K,K, RELFRQ 210 

620 FORMAT! 1H1, 15, 18H DEPENDENT TERM ,15, 59H IF THE DISTRIBUTION 211 

XWERE NORMAL EACH CELL WOULD CONTAIN F6.2, 8H COUNTS ) 212 

DO 640 1=1, CELLS 213 

RCT ( I ) = OBS ( I , K ) 214 

640 CONTINUE 215 

CALL HISTIK, RCT, CELLS) 216 

650 CONTINUE 217 

RETURN 218 

END 219 
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$1 BFTC RECTXX 


SUBROUTINE RECT ( I ROW , J JCOL , I MAX , JMAX , A , FM T ) I 

DIMENSION A! IMAX.JMAX) ,FMT( 14) ,X0UT(8) 2 

DATA J8 / 8/ 3 

LOGICAL OUT 4 

OUT =. FALSE. 5 

JT I MES=0 6 

JCOL= J JCOL 7 

5 JNXT= JCOL— J8 8 

IFIJNXT) 10,20,30 9 

10 JP= JCOL 10 

GO TO AO 11 

20 JP= J 8 12 

GO TO 40 13 

30 JCOL= JNXT 14 

JP = J8 15 

GO TO 50 16 

40 OUT= .TRUE. 17 

50 DO 100 1 = 1, I ROW 18 

DO 60 J=1,JP 19 

JJ=JTIMES +J 20 

60 XOUT ( J ) = A( I,JJ ) 21 

WRITE ( 6, FMT ) I , ( XOUT ( K ) , K= 1 , JP ) 22 

100 CONTINUE 23 

I F ( OUT ) RETURN 24 

WRITEI6, 110) 25 

110 FORMAT! 1H /1H ) 26 

JTIMES= JTIMES +JP 27 

GO TO 5 28 

END 29 


SIBFTC LOCXXX 

SUBROUTINE LOCI I , J, IR ) 1 

IX= I 2 

JX= J 3 

20 IF(IX-JX) 22,24,24 4 

22 I RX= IX + (JX*JX-JX)/2 5 

GO TO 36 6 

24 I RX= JX + (IX*IX - IXI/2 7 

36 I R= IRX 8 

RETURN 9 

END 10 
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n n ° ° no o oo on onnonooooooono 


$ I BFTC BORCXX 


SUBRCUT I NE BORD 
PURPOSE 

TO COMPLETE THE INVERSION OF A SYMMETRIC POSITIVE DEFINITE 
MATRIX A OF ORDER N GIVEN THAT THE UPPER LEFT SUB- 
MATRIX OF ORDER N-l HAS ALREADY BEEN INVERTED. 

SUBROUTINES NEEDEO 
LOC 

REMARKS 

ONLY THE UPPER TRIANGULAR PART OF A IS STORED AS A 
VECTOR IN THE ORDER A C 1 1 1 ) » A C 1 f 2 1 t A ( 2 , 2 > * A ( l ♦ 3 ) f . . .ETC 
SUBROUTINE BORD < I ORDER , A ) 1 

3 

10 

DIMENSION BETA( 60) f A( 1 ) 11 

12 

13 

ALPHA= 0.0 14 

NM 1= IORDER-1 15 

IF(NMl) 100,100,200 16 

100 A( 1) = 1.0/A! 1) 17 

GO TG 600 18 

200 M=NMl*( NMi + 1 )/2 19 

LEN = M + I ORDER 20 

21 

DO 400 1 = 1 1 NM 1 22 

BETA ( I) = 0.0 23 

M I = M ♦ I 24 

DO 350 J = 1 , NM 1 25 

CALL LOC( I , J, II ) 26 

MJ= M+J 27 

BET A ( I ) = BETAUJ-A! in*A(MJ) 28 

350 CONTINUE 29 

ALPH A= ALPHA + A (MI )*BETA( I ) 30 

400 CONTINUE 31 

32 

33 

ALPHA = ALPHA + A ( LEN ) 34 

RALPH A= 1.0/ALPHA 35 

A ( LEN ) = RALPHA 36 

37 

DO 500 1=1, NM1 38 

DO 500 J=1,I 39 

CALL LOC ( I ,J, II ) 40 

A ( I I ) = A ( 1 1 ) + BETA! I)*BETA< J)*RALPHA 41 

500 CONTINUE 42 

43 

DO 550 J= 1 , NM 1 44 

M J = M4J 45 

A ( M J ) = BETA! J)*RALPHA 46 

550 CONTINUE 47 

48 

49 

600 CONTINUE 50 

RETURN 51 

END 52 
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APPENDIX B 


BORROWED ROUTINES 


Some of the routines used in the program were taken from the literature. Both 
BSTVXTX and TRIANG are by Webb and Galley (ref. 9), and EIGEN and HIST are from the 
IBM programmer's manual (ref. 10). 

Listing of INVXTX and TRIANG are given here, as follows: 


41BFTC INVXXX 


SUBROUTINE 1NVXTXIA, NN, D, FACT! 

C 

C ASSUMES THE MATRIX A IS SYMMETRIC ANO POSITIVE DEFINITE » AND ONLY 

C THE UPPER TRIANGLE IS STORED AS A ONE-DIMENSIONAL ARRAY IN THE 

C ORDER All, 11. All, 21* A(2,2), All, 3), A(2,3), A(3,3), A ( N , N ) • 

C NN IS THE ORDER N OF THE INPUT MATRIX A, 

C D IS ION EXIT) THE DETERMINANT OF A, DIVIDED BY FACTOR**NN. 

C 

DIMENSION A(l) 

D * 1,000 

N = NN 
ITR1 = 0 
FACTOR = FACT 
DO l A 5 K = l,N 
C 

ITR1 = ITR1+K-1 
KP1 * K + l 
KM1 = K— 1 
KK = ITR1+K 

C CONTINUED PRODUCT OF PIVOTS 

D = D*A(KK)/FACTOR 
PV = 1 • 000/ A ( KK ) 

C 

I TR2 * 0 

IF ( K— 1 I 150,80, 50 
C 

C REDUCE TOP PART OF TRIANGLE, LEFT OF PIVOTAL COLUMN 

50 DO 60 J = 1 * KM 1 
ITR2 = ITR2+J-1 
KJ = ITR1+J 
F = A ( K J ) *PV 
DO 60 1 = 1, J 
IJ = ITR2+I 
IK = ITR1 ♦ I 

60 A( I J) = A( I J) + A ( IK ) *F 
C 

IF (K-N ) 70,120, 150 
C 

C REDUCE REST OF TRIANGLE, RIGHT OF PIVOTAL COLUMN 

70 I TR2 = ITR1 
80 DO 110 J=KP 1 , N 
ITR3 = ITR1 
I TR2 = ITR2+J-1 
KJ = ITR2+K 
F = A ( K J ) *PV 
DO 100 1=1, J 
IF (l-K) 90,100,95 
90 IJ * ITR2+I 
IK = IT R l ♦ I 
A( IJ) = A( IJ) - A( IK )*F 
GO TO 100 
95 IJ = ITR2 + I 

I TR3 = I TR 3 + 1-1 
IK = IT R3 + K 
A( IJ) = A( IJ) - A( IK )*F 
100 CONTINUE 
110 CONTINUE 
C 

C 01VIDE PIVOTAL ROW-COLUMN BY PIVOT, INCLUDING APPROPRIATE SIGNS 

120 I TR2 * I TR 1 


2 

3 

4 

5 

6 

7 
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19 

20 
21 
22 

23 
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27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 
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56 
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60 
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i i min ii i is 


DO 140 1=1, N 61 

IF (I-K) 125,130,135 62 

125 IK = I TR 1+ I 63 

At IK) = —A ( I K )*P V 64 

GO TO 140 65 

C (REPLACE PIVOT BY RECIPROCAL) 66 

130 A ( KK ) = PV 67 

GO TO 140 68 

135 I TR2 = ITR2+I-1 69 

K I = IT R2+K 70 

A ( K I ) = A ( K I )*P V 71 

140 CONTINUE 72 

C 73 

145 CONTINUE 74 

C 75 

150 RETURN 76 

END 77 


SIBFTC TRIANX 


SUBROUTINE TR IAN G ( A , NN , NKOL , FORMAT ) l 

DIMENSION FORMAT ( 1) 2 

DIMENSION A( 1) 3 

1 FORMAT (1H1) 4 

3 FORMAT! 1H /1H /1H ) 5 

N = NN 6 

NCOL = NKOL 7 

KLUMPS = N/NCOL 8 

C 9 

KEEPTR =0 10 

K 1 = 1 11 

K2 = NCOL - 1 12 

K3 = NCOL 13 

IF (KLUMPS .EQ. 0) GO TO 120 14 

C 15 

DO 90 KLUMP=1, KLUMPS 16 

ITR 1 = KEEPTR 17 

I = -1 18 

ILO = (KLUMP-1)*NC0L + ITR1 + 1 19 

DO 30 K=Kl,K2 20 

1 = 1+1 21 

I TR 1 = I TR 1 + K - 1 22 

ILO = ILO + K - 1 23 

I H I = ILO + I 24 

30 WRITE (6, FORMAT) K, (A(J), J=IL0,IHI) 25 

KEEPTR = I TR 1 + K2 26 

00 60 K=K3,N 27 

ITR1 = I TR 1 + K - 1 28 

ILO = ILO + K - 1 29 

I H I = ILO + NCOL - 1 30 

60 WRITE (6, FORMAT) K, (A(J), J=ILO,IHl) 31 

K 1 = K1 + NCOL 32 

K2 = K2 + NCOL 33 

K3 = K3 + NCOL 34 

90 WR I TE ( 6 , 3 ) 35 

C 36 

120 I TR 1 = KEEPTR 37 

IF ( K 1 ,GT. N) GO TO 180 38 

1 = -1 39 

ILO = KLUMPS+NCOL + ITR1 + 1 40 

DO 150 K=K 1 » N 41 

1=1+1 42 

ITRl = I TR 1 + K — 1 43 

ILO = ILO + K - 1 44 

I H I = ILO + I 45 

150 WRITE (6, FORMAT) K, (A(J), J=ILO,IHI) 46 

C 47 

180 RETURN 48 

END 49 
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